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CHAPTER I 


INTRODUCTION 

As the exploration of space progresses from the near earth 
environment to the moon and the planets of the solar system, study 
of the atmospheric entry of objects in excess of escape velocity 
(about 11 km/sec) becomes necessary. In addition to stuides con- 
cerning manmade objects, there is considerable interest in the entry 
of meteoroids into the earth's atmosphere at velocities from 20 to 
70 km/sec. At these large speeds radiant energy transfer is an 
important factor governing the behavior of the hot shock layer gas 
enveloping the object. 

Consequently, a number of investigators have addressed them- 
selves to the problem of the radiating shock layer. The first 
analyses assumed that the flow processes were uneffected by (or 
uncoupled from) the transfer of energy by radiation ( see for example 
refs. 1 and 2). Thermodynamic and flow properties were calculated 
neglecting radiation. The radiant energy flux was then calculated 
from measured or theoretically determined optical properties for 
these conditions. While this approach provides acceptable engineering 
estimates at speeds less than escape velocity, it is not sufficient to 
describe the effects of radiation at higher speeds. The next step was 
to take into account the loss of energy from the shock layer due to 
radiation. This cooling of the shock layer tends to reduce the 
emergent radiant energy flux. This reduction is often termed 
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"radiation decay." Radiation decay was studied by a number of authors 
(see, for example, refs. 3 - 8 ). All of the cited works, with the 
exception of reference 8, used the transparent approximation*, which 
neglects absorption within the shock layer, and assumed that the 
radiation cooling of the shock layer gas was small and did not influ- 
ence the mass transport. The process of absorption by a gray** gas 
was studied in reference 8. However, the flow model used in that 
investigation only roughly approximates the flow in the stagnation 
region of a shock layer. Consequently, the analysis was unable to 
describe details of conditions in the shock layer or to provide reli- 
able quantitative results. 

Perhaps the most ambitious analysis to appear to date is the 
work of Howe and Viegas (ref. 9)* They obtained numerical solutions 
to the integrodifferential system of equations governing the flow in 
the stagnation region including the effects of radiation decay, 
absorption by a gray gas, viscosity, and surface mass injection. An 
indication of the complexity of this numerical approach is the reported 
computation time for a single example of 5 hours on the IBM 7090 
electronic digital computer. 


*So called because the shock layer gas is considered to be 
transparent to its own radiation. 

A gray gas is one for which the optical properties are 
independent of the photon energy or wavelength. 
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All of the works discussed above are restricted to velocities 
less than about 20 km/ sec, although the work of Howe and Viegas was 
so restricted simply because they did not choose to make calculations 
for higher velocities. Fay, Mbffatt, and Probstein (ref. 10) under- 
took an analysis of meteoroid entry, in the speed range of 20 to 
70 km/sec. Since they were interested only in obtaining upper bound 
estimates of radiant heating, they ignored radiation decay and absorp- 
tion (except that they did not allow the radiant energy flux to 
exceed the black-body limit) both of which can be quite important 
at these speeds. 

While the existing studies (which include many works in addition 
to those cited) have contributed much to the qualitative and quanti- 
tative understanding of the physical processes taking place in 
radiating shock layers, a great amount of work remains. For example, 
parametric studies of absorption in a realistic shock layer flow are 
lacking, the effects of surface reflectivity have been generally 
ignored, and there have been no reported attempts (at least in the 
knowledge of this investigator) to study shock layer gases with non- 
gray optical properties.* 

The investigation reported herein was undertaken to provide a 
parametric study of the influence of radiation on blunt objects 
large and small, travelling at speeds up to 70 km/ sec. The approach 


*Lick (ref. 11) and Greif (ref. 12) considered nongray optical 
properties in their studies of combined radiation and conduction. 
Their results indicated that nongray effects can be significant. 
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has been to seek simple approximate solutions where available in the 
hope that they would lead to a better understanding of the physical 
processes involved. The parameters to be studied include the radia- 
tion cooling parameter e, which characterizes the relative importance 
of radiation as an energy transport mechanism compared to convection, 
the Bouguer number, which indicates the importance of absorption in 
radiant transport, the surface reflectivity, indicative of the ability 
of the surface of the object to accept the incident radiant energy 
flux, and the spectral variation of the absorption coefficient. (There 
is no single quantity or even group of quantities which characterizes 
the important effect of spectral variation on the flow. ) Definitions 
of these parameters and their role in influencing the flow will be 
discussed in greater detail in subsequent chapters. 

In order to facilitate this investigation without sacrificing 
physical significance, analysis is limited to the stagnation region 
and the following conditions are assumed to apply: (l) the shock 

layer gas is in local thermodynamic and chemical equilibrium, (2) the 
body geometry is axisymmetric, (3) there is no mass addition to the 
flow from the body surface, (k) the thicknesses of the shock and the 
viscous boundary layer are small in comparison to the shock standoff 
distance, and (5) absorption in the free stream ahead of the object 
is negligible. 

In this investigation, solutions will be obtained for four 
limiting cases of the radiation cooling parameter and the Bouguer 
number. The first of these, which is presented in chapter III, is 
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a small perturbation expansion in the radiation cooling parameter e 
valid when the influence of radiation is small. The second solution, 
presented in chapter IV, holds when the shock layer is optically thin. 
This solution is presented as a small perturbation expansion in the 
Bouguer number. A solution valid when the shock layer is optically 
thick (Bouguer number » l) and the final solution, which is 
restricted to the case when radiation is the principal mode of energy 
transport within the shock layer, are presented in chapters V and VI, 
respectively. The first and second solutions have been formulated to 
include the effects of nongray radiation. The third and fourth 
solutions are restricted to the gray case. In each of the four 
limiting cases, it is possible to approximate the governing integro- 
differential system of equations by a purely differential system 
which leads to a singular perturbation problem. 

The results obtained by means of the various approximations are 
combined in chapter VII to give the radiant heat transfer rate and 
an estimate of the effect of radiation on the convective heating rate 
at the stagnation points of blunt objects traversing a gray model 
earth atmosphere. The effects of the nongray character of air on 


these results is discussed. 
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CHAPTER II. STAGNATION MODEL FOR 
A RADIATING SHOCK LAYER 

A. Fundamental Equations of Radiation 
Gas Dynamics 

Prior to setting up a particular flow model for the problem at 
hand, it is desirable to examine briefly the fundamental equations 
of radiation gas dynamics. An excellent discussion of these equations 
has been presented by Goulard in the volume "High Temperature Aspects 
of Hypersonic Flow" (ref. Ij), and the reader is referred to this 
vork for a more detailed exposition. 

In the first chapter, it was indicated that the studies of this 
paper are limited to the steady flow of gases in local thermodynamic 
and chemical equilibrium. In addition, the effects of radiation 
pressure and radiation energy density are ignored. These effects 
are important only when the radiant energy flux is extremely large 
as it is deep in the interior of a stellar atmosphere. Finally, the 
presence of external forces, such as gravity and electromagnetic 
forces, are neglected. With these restrictions in mind, the conser- 
vation equations for a radiating gas can be written 



(Continuity) 


(2.l)* 


The double subscript notation is employed. 


6 



7 


PU j u i, j = p ,i + T ij,j (Momentum) (2.2) 

pu i h t,i - T ij) -j -«°,i -%,i (a»w) (2-3) 

where the quantity h^ is the total specific enthalpy of the gas 

h t = h + i u i u i (2.u) 

The static specific enthalpy h includes the chemical energy of the 
gas in terms of the heats of formation of the various gaseous species* 
An expression relating the thermodynamic variables is needed to 
complete the set of equations* A convenient form is 

h = h(p,p) (2*5) 

The molecular transfer processes are represented by the classical 
expressions 

t « ■ + -N 6 u \,k < 2 - 6) 

4 " ’ k eff T 'i (2-7) 

The quantity k is an effective coefficient of heat conduction 
which includes the effects of energy transport by molecular collisions 
and by the diffusion of reacting species. These two processes cam be 
lumped together like this only when the conditions of local thermo- 
dynamic and chemical equilibrium hold (see ref. 14). 



8 


The radiant energy flux vector q ± R is defined in terms of the 

radiation intensity J, . 

A 

_ 00 

q i R s f o ax ’ 4i s f Un h h * < 2 - 8 > 

and is the rate of flow of radiant energy per unit area across an 

element of area whose normal points in the ith direction. Hie 

quantity l ^ is the direction cosine between the direction of a 

single beam of intensity J. and the ith direction. J, can be 

A A 

determined from the conservation equation of radiation transfer 

dT= " pP a[ J A (2*9) 

where 0^ is the mass extinction coefficient. It is composed of 

the mass absorption coefficient and the mass scattering 

coefficient ol. 

A 

\ \ 

The ratio of mass emission coefficient to the mass extinction 

coefficient 3^ is often called the source function . 

For nonscattering media in a state of local thermodynamic 
equilibrium = oj , the source function reduces to the Planck 


function 
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2hc 2 fhc/WcT 
'A= x 5 1 


( 2 . 11 ) 


provided that the mass absorption coefficient includes the 

effects of induced emission. Here h and k are the Planck and 
Ppl~h 7.TTifl.n n constants, respectively, and c is the speed of light. 
Throughout the remainder of this paper, it is assumed that the gas 
in the shock layer is nonscattering. This assumption is reasonable 
as the number of large solid particles which might scatter radiation 
is expected to be negligible in the shock layer. A few such particles 
might exist in the cooler regions of the boundary layer adjacent to 
the body surface as a result of "spalling" of this surface. However, 
their presence could be accounted for, if necessary, by changing the 
effective reflectivity of the body surface. 

At the extremely high shock layer temperatures for which the 
gas is multiply ionized and free electrons are plentifully, Thomson 
scattering can become important. For example, Kivel and Mayer 
(ref. 15) show that scattering cannot be neglected when the tempera- 
ture reaches 550,000° K at densities less than about 0.01 of the sea 
level value. 

For the nonscattering case the intensity of radiation at a 
point M in the direction i* follows from a formal integration 
of equation (2.9) 




*A1 though the terms "optical thickness" and "optical path length," 
long established in astrophysics! literature, seem to imply a dimension 
of length the quantity ^ is dimensionless and is indicative of 

the number of photon mean free path lengths in the physical distance 
between M and P. 
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The quantity (Q) represents the contribution to the 

intensity at point M from the boundary and, in general, includes 
emission from the surface, reflection from the boundary of radiation 
originating from within the region, and transmission through the 
boundary of radiation originating from without the region. 

The integral term represents the summation of the contributions 
from all points P along the beam reduced by the attenuating factor 
exp(-T^ which acounts for absorption by the intervening matter. 

The divergence of the radiation flux vector can be found with 
the aid of solution (2.12) with the result 

( (s) 1 

_ , poo n exnl -t' ' 

- i, H)m - L K)m J t K b a)p dV(p)aA 

S6r 

( 2 . 13 ) 

-/»« s k ^ (q > c ° s 9 ds(e)dA 

The integrations over the volume V of the gas and the area A of 
the bounding surface include only those portions of the volume and 
surface which are visible to an observer stationed at point M. 

B. Stagnation Flow Model 

The study of three-dimensional flow of compressible gas in the 
vicinity of the forward face of a blunt body cannot be reduced via 
transformation to the study of an equivalent one -dimensional system 
as can be done in the incompressible case. However, available 
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numerical solutions (see for example refs. 15-18) indicate that for 
all practical purposes a reduction from a three -dimensioned, to a 
nearly equivalent one -dimensional problem can be carried out in the 
stagnation region. The reason that this simplification can be 
applied is that the flow behind a strong bow shock is nearly incom- 
pressible in the stagnation region. Also the various thermodynamic 
properties are nearly independent of the lateral or radial coordinate. 
While the same arguments apply in the stagnation region of a 
radiating shock layer, it is not possible to postulate the existence 
(even approximately) of a one -dimensional solution solely on this 
basis. Some additional assumption is required regarding the effect 
of the far -field on the radiant heat flux and its divergence. This 
effect, of course, cannot be obtained a priori as it depends on the 
solution to the entire flow field. Fortunately, the shock layer is 
thin and only a small portion of the radiant energy emitted by gas 
removed from the stagnation region actu all y passes through the 
stagnation region. If absorption is small, only a small portion of 
this is absorbed in the stagnation region. If , on the other hand, 
absorption is large, the beam is greatly attenuated when it reaches 
the stagnation region leaving only a small portion of the energy 
which started the journey to be absorbed in the stagnation region. 

The divergence of the radiant flux is influenced only by the amount 
of energy absorbed and emitted. Consequently, the far-field effect 
on the divergence of the radiant flux is a result of that small 
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portion of radiant energy originating in the far-field and absorbed 
in the stagnation region. In the transparent and optically thick 
limits , this effect of the far field vanishes. 

It would appear from the above discussion that a stagnation 
model for a radiating gas can be postulated as long as the assumptions 
concerning the far-field are not grossly unrealistic. In what follows, 
a particular stagnation model will be formulated and an estimate of 
the inaccuracy resulting from the assumption concerning the far field 
will be obtained. 

A schematic of the flow in the stagnation region of a blunt 
body is shown below 
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At very high speeds, the ratio of the shock standoff distance A to 
the shock radius R is very much smaller than one (a typical value 

S 

is 0.05)* Under these conditions, the geometry of the stagnation 
region closely resembles a plane parallel gas slab. In addition, 
the enthalpy in the shock layer varies slowly with respect to r/A 
so that the stagnation region may be approximately represented by 
a gas slab in which the thermodynamics as well as the geometry is 
one-dimensional . 

As a result of the above considerations, the model described 
below has been chosen to represent the flow of a radiating gas in 
the stagnation region of a blunt object. The model consists of an 
axially symmetric flow impinging upon an infinite flat plate normal 
to the stream direction. At a plane which is parallel to the plate 
and a distance A in front of it, the gas is suddenly raised to a 
total specific enthalpy of ^ tT. The plate is held at a constant 
temperature T^. A sketch illustrating the geometry of the flow 
model is shown below. 



Axis of 
symmetry 


Shock 


Body 
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The general equations of motion (eqs. (2.1) and (2.5)) when 
specialized to the asisynmetric geometry become 
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(2.14) 

(2.15) 

(2.16) 


(2.17) 


where and q z are the r- and z -components, respectively, of 

the heat flux vector which includes, conduction, diffusion of reacting 
species, and radiation. The stress components are given by the 
expressions 


T 


rr 


T 


zz 











(2.18) 

(2.19) 
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( 2 . 20 ) 
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dw 
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( 2 . 21 ) 
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The equation of state is 


h s h(p,p) (2.22) 

In order to completely specify the problem, a consistent set 
of boundary conditions must be provided. The kinematical conditions 
on the velocity are 


v(r,0) = 0 (2.23) 

p(r,A)w(r,A) = -p W (2.2U) 

The first of these conditions restricts the analysis to one 
for which there is no injection from the surface of the object. When 
gas injection is important, it is necessary to replace the zero on 
the right-hand side of equation (2.2p) with w^, the normal velocity 
of the gas at the wall. The second condition was obtained from 
continuity across the normal shock at the stagnation point. A third 
kinematical condition is introduced here in order to relate the 
standoff distance A to the variation of the tangential velocity 
along the surface z = A. This variation in velocity is taken to be 
equal to that behind the near normal portion of a spherical shock, 
that is 

u(tA) = cos 3 ** W m ~ ( 2.25) 

where 3 is the local inclination of the shock from the free 
stream direction and R g is the radius of the spherical shock. 
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The dynamical "no- slip" condition at the surface is 


u(r,0) = 0 (2.26) 

The conditions on the enthalpy and pressure are 

h^(r,A)=iw? (2.27) 

^ C- 

(r,0) = h w (2.28) 


p(r,A) 


* P. ^ (1 - X) 1 - 


(2.29) 


where X » j p(0,A) is the density ratio across the normal shock. 

Condition (2.27) comes from the conservation of energy across a strong 

normal shock and does not take into account absorption in the free 

stream of radiant energy emitted by the shock layer. Condition (2.28) 

restricts the analysis to those conditions at which a temperature 

"jump" or discontinuity is not present at the body surface. Such a 

"jump" can occur only when the molecular mean free path in the gas 

is not negligible in comparison to the characteristic length of the 

domain (in this case, the thickness of the ther mal boundary layer) . 

Condition (2.29) is the pressure distribution behind the near normal 

portion of a spherical shock of radius R . 

s 

In addition to the boundary conditions listed above, boundary 
conditions on the radiant energy flux must be specified. These 
conditions are: 
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(1) The boundary at z = A (which corresponds to a bow 
shock) is transparent. 

(2) There is no radiant energy transfer from the free stream 
to the shock layer. 

(5) The boundary at z = 0 (which corresponds to the body 
surface) is cold and reflects diffusely and independently of wave- 
length a fraction r of the incident radiation. 

V 

The statement (contained in condition ( 3 )) that the body surface 
is cold means that emission from the body surface has a negligible 
influence on the gas in the shock layer. When the hot (temperatures 
in excess of 10,000° K) shock layer is optically thin emission from 
the relatively cool (temperatures less than 4,000° K) body surface 
may be comparable to emission from the shock layer gas. However, 
because the shock layer is optically thin very little of the radiant 
energy emitted at the body surface will be absorbed by the shock 
layer gas. On the other hand, when absorption in the shock layer 
is important the shock layer gas emission will approach the black- 
body value corresponding to the high shock layer temperature. Since 
black -body radiation is proportional to the fourth power of temperature 
the gas emission from an optically thick layer will greatly exceed 
the emission from the body surface. Thus whenever the body surface 
temperatures are small compared to the 6hock layer gas temperatures 
the influence of emission from the body surface on the shock layer gas 
is unimportant. 
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Prior to assuming that the shock layer is one-dimensional it 
is necessary to specify whether the body surface reflects diffusely, 
specularly, or in some combination of the two. However, in a one- 
dimensional system this specification is superfluous because the 
difference in effect of the two types of reflectivity vanishes . 

Since the surface reflectivity of most solid materials at high 
temperatures varies little with the wavelength, the assumption that 
the surface reflectivity is independent of wavelength provides a 
simplification in the analysis without sacrificing physical 
significance. 

It can be seen from the definition of the total enthalpy 



and the boundary conditions (2.2*0 and (2.27) that the magnitude of 

O 

the kinetic energy in the shock layer is order X compared with 
the static specific enthalpy. For a strong shock, which is the only 
case of interest here, 0.05 is a typical value for X, the density 
ratio across the shock. As a consequence of the above the kinetic 
energy terms will be neglected in the subsequent analysis. The 
viscous dissipation terms (the last three terms on the right-hand 
side of equation (2.17)) will also be neglected because only kinetic 
energy is dissipated through the action of the viscous forces. 

It is desired that the solutions to the one -dimensional model 
represent, as closely as possible, the phenomena in the stagnation 



region of a blunt body. For simplicity, the blunt body geometry, 
flow field, and thermodynamic properties are considered to be 
axially symmetric about the stagnation streamline. Expanding the 
solutions for the real blunt body problem in terms of the radial 
coordinate r and arguing on physical grounds that w(r,z), p(r,z), 
p(r,z), and h(r,z) are even functions of r while u(r,z) is 
odd, gives 

w = (z) + 0(r 2 ) 

u = ru^ (z) + 0(r :> ) 

P = P^ (z) + 0(r 2 ) (2.30) 

p = p^ (z) + 0(r 2 ) 
h . h^°^ (Z) + 0(r 2 ) 


In addition, the heat flux components will have the form 


= q* 0) « + 0(r 2 ) 


(2.31) 


Qj. = ( z ) + ©(r 5 ) 

p 

Neglecting terms of order r and higher restricts the solutions 
to the vicinity of the stagnation point. Since stagnation region 
solutions are desired, it will be assumed that the solutions in the 
plane parallel model have the functional forms of equations (2.30) 
and (2.31) truncated after the linear term in r. For these assumed 



forms, the continuity equation (2.14) requires 


pu = rg'(z) 
pv = -2g(z) 


(2.J2) 


That portion of the heat flux due to conduction and diffusion 
of reaction species is proportional to the enthalpy gradient, that is 


_ c dh c dh 

q ~ -r— , a ~ ~ — 

z dz' t dr 


From conservation of energy across the near normal portion of a strong 
spherical shock 


h ~ 



1 



Thus 


Hh ^ c ~i^ r#- 


Comparing terms that appear in the energy equation one finds that 

2 


V + V ~ (a 

or r 


« 1 


Thus, the terms containing can be neglected in the formulation 
of the stagnation flow model. 

If the shock layer is optically thick, that is the photon mean 
free path is very small compared to the shock standoff distance, the 
radiation flux terms take on the same form as the conduction 
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* R/ R / 

terms and / Sr + /r may be neglected. On the other hand, 

if the gas is not optically thick, this simple order of magnitude 
analysis no longer suffices because the divergence of the radiant 
energy flux depends not only on local conditions, but on conditions 
throughout all of the shock heated gas which can be seen by an 
observer located at the point in question. 

Calculations were made of the divergence of the radiant flux 
for a gray isothermal gas in a shock layer formed by two concentric 
spherical surfaces with a standoff distance to shock radius ratio 
of 0.05 • A sketch showing the volume of gas which contributes to 
the radiant flux at a point on the stagnation streamline is shown 
below . 



According to the Rosseland or diffusion approximation. 
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The largest difference between this result and the divergence 
of the radiant flux for a plane -parallel layer occurred adjacent to 
the wall for an optical thickness of about 0.13* The difference 
amounted to 2.2 percent of the value for the plane -parallel layer. 

A second set of calculations was made to determine the effect 
of a nonuni form enthalpy distribution, in the lateral direction, 
on the magnitude of the divergence of the radiant flux. The enthalpy 
distribution was given by 


h(r) 


h(0) 



for r < V^ R S 


0, for r > -\/2 R s 


(2.33) 


V. 

This expression approximately corresponds to the enthalpy distribu- 
tion in the shock layer about a spherical body. The absorption 

•ft 

coefficient was assumed to vary as the third power of the enthalpy 

i 

and the shock standoff distance to shock radius ratio, Aj R g , was 
chosen to be 0.05. A comparison of calculations for a plane -parallel 
layer in which the enthalpy was assumed to vary according to 
equation (2.33) and of calculations for a plane -parallel layer in 
which the lateral enthalpy distribution was uniform (i.e., h(r) = h(o)) 
indicated that the largest difference in the magnitude of the 


This variation is consistent with the correlations of the optical 
properties of air to be discussed in a subsequent section of this 
chapter. 
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divergence of the radiant flux occurred for a shock layer optical 
thickness of about 0.1. This difference amounted to 2.3 percent of 
the value for the uniform distribution. These results are inde- 
pendent of the value of h(0). 

Since the errors in the divergence of the radiation flux for 
the one -dimensional shock layer due to the separate effects of 
geometry and nonuniforra lateral distribution of enthalpy are small, 
their combined effect should be given approximately by the sum of 
the separate effects. That is, the maximum error due to the combined 
effects of geometry and nonuniform lateral distribution is probably 
not much greater than 5 percent for A^R g = 0.05. This, of course, 
does not imply that the final results for the enthalpy (for example) 
would be in error by 5 percent but only that one term in the energy 
equation is in error by 5 percent. In any event, the results of the 
calculations mentioned above are considered to give sufficient 
justification for choosing the plane -parallel layer as a model for 
the stagnation region of a blunt body. 

The expression for the radiant energy flux is more seriously 
affected by the plane -parallel layer assumption than is the divergence. 
For example, Koh (ref. 19) has computed the radiant flux at the vail 
for an isothermal shock layer formed by two concentric spheres. For 
a shock standoff distance to body nose radius ratio of 0.05 and a 
vanishingly small value of optical thickness the result is about 
17 percent less than for a plane- parallel isothermal, layer of equal 
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optical thickness. This difference decreases with increasing optical 
thickness . Koh also computed the effect of nonuniform lateral en- 
thalpy distribution using an assumed enthalpy distribution similar 
to that given by equation (2.33). He found that the flux at the 
wall for the nonuniform distribution was about 1.5 percent less than 
that for an isothermal layer for a shock standoff distance to shock 
radius ratio = 0.05 and a vanishingly small optical thickness. 

As expected, the difference decreases as the optical thickness 
increases. 

It is apparent from Koh's results, that an accurate estimate 
to the rate of radiant heat transfer to the stagnation point cannot 
be obtained through the use of the plane -parallel layer approximation 
unless some correction factor, which takes into account the actual 
geometry of the shock layer, is employed. However, because this 
investigation is concerned with obtaining a general understanding 
of the problem of radiating shock layers rather than specific numerical 
results, such a correction factor will not be used herein. 

At this point, it is convenient to introduce the variable 
transformation 

H = f Pdz ( 2 . 3*0 

u o 

The new variable i) is often called the Dorodnitsyn variable. 

Under this transformation, the normal and tangential velocity 
components become 
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w 


-2f(q) 


dz 

dr] 


(2.35a) 


The two momentum equations 


u = rf ’ (tj) (2.35b) 

(eas. (2.13) and (2.16)) take the form 


(p f"(n)] + 2f(n) f”(n) - [i'Wj = ± 5F (2 ' 56) 


and 


Phf"(n) - 2p n^f'(i}) - ^ p’(n) f(q)j 


+ p 


' - 1 p'(n) f(n) 


- 2f( n ) f * ( tj ) 


(2.37) 


+ | p'(n)[f(n)] = + |pp’(n) 


An order of magnitude analysis of equation (2.37) indicates 
that p'(q) is order X or Re"'*', whichever is larger. Since both 
X and Re ^ are very small compared to unity equation (2-37) will 
be replaced by the simple approximate expression 

P’(ti) = 0 (2.38) 


Thus, the pressure is a function of r only. In particular, the 
strong shock relations for the near normal portion of a spherical 


shock give 
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To first-order in r 


* 


1 dp 
r dr 


-2P_ 



(1 - X) -L 

R 


so that equation (2.56) becomes 


[puf + 2f (n) r»(,) -[f(i0 2 - - M ” f ' X) 

pR s 

(2.40) 

Under the foregoing assumptions and the coordinate transformation, 
the energy equation (eq. (2.17)) becomes 

“2f ( T) ) h * ( r) ) + q'(q) = 0 (2.4l) 

The boundary conditions are 



f(0) = 0 

(2.42) 


f’(o) = 0 

(2.43) 


f (\) - \ »» a . 

(2.44) 

* 

r x W 

f 'W - 1 

(2.45) 


h(0) = h 

V 

(2.46) 


h(n /v ') = W 2 
l 'A J 2 «> 

(2.47) 
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where 


1 A = / Pdz 

^ J o 


(2.43) 


The heat flux terra q(rj) in the energy equation (2. Hi) is 
composed of a combined conduction and diffusion terra 


q C (T!) = -pkT'(n) - - h ’ ( tj ) 


and a radiation terra 


(2.49) 


r,\ 

i 




i B-. 
A 


(h) E L> ( T X • \) dt X 


(2o0) 


“ 2r w E J ( T a) j o A B A (**) E 2W dt X> d ^ 


This radiation terra is representative of the case of a plane -parallel 

geometry with a transparent wall (shock) and a cold wall, which 

reflects diffusely and independently of wavelength a fraction r of 

the incident radiation, separated by a nonscattering, nongray gas. 

The variable t. is called the "optical path length" and is defined 
A 

by the expression 
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r z 

T = 1 

* Jo 

Expression (2.p0) was specialized from the more general expression 
of Goulard (ref. l). Goulard derived the expression for the radiant 
flux in a plane parallel geometry with arbitrary reflecting, absorbing, 
and emitting walls separated by a nonscattering, nongray gas. His 
expression was restricted to isotropically emitting and diffusely 
reflecting surfaces and hence, so is equation (2.50). However, this 
restriction is of little consequence in this problem because emission 
from the wall will be neglected (the wall is cold) and there is no 
difference in effect between specular and diffuse reflection in the 
one -dimensional case. 

The first term in equation (2.p0) represents the radiant energy 
flux passing through the plane r, = const, and which originated in 

A 

the region between this plane and the shock at t- = t\ . This 

A 

radiant flux has been attenuated by partial absorption in the inter- 
vening gas. The second term represents the radiant energy flux 
passing through the plane = const, and which originated in the 

region between this plane and the wall at = 0. This flux has 

A 

also been attenuated by partial absorption in the intervening gas. 

The last term represents the radiant flux passing through the plane 
= const, and which was reflected from the wall and attenuated 
by the intervening absorbing gas. 


pK A 


dz = 


I V * 


(2.51) 
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Substituting the expressions for the energy flux equations (2.49) 
and (2.90) into the energy e uation (2.41) gives 


2f ( tj) h 1 (t)) + 




+ I 


H- 0 


(2.52) 


where the divergence of the radiant flux is represented by the 
integral terra 


b (t) 


*f 0 ^ V v) e i(J t a ( ti) ‘ 

r v v>'> E 2 ^(T)JdXx 


(2.53) 


The final step in the derivation is to reduce the equation to 
nondimensional form. For this purpose, the following set of non- 
dimensional quantities is introduced 


i ■ P s A A i, f(i) = I p.vW, h(n) . I v L 2 h(n) 


2 00 


If ■ «* - (p, p' 1 - p^w 

I[3]= 2aT e k ^ = K p, T x (t|) = k p ? A (i) 


(2.54) 


oT 


s /-X 2 af l - X X A A ^ 

A r s y 


B x<^ = -T- a - 6 [ : 


2 
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u c 

s P« P W A. 

= — £ > Re = ~ ° ” , Pe = Pr Re 

s K C s 6 S S 


( 2 . 5 ^) 


4cT 


r = 


l 3' *P = p s K P V € = lc p 


00 00 


The subscripts » and s indicate conditions evaluated in the 
undisturbed free stream and immediately behind the shock, respec- 
tively. The quantity A^ is the shock standoff distance for 
the nonradiating (or adiabatic) shock layer. The property variations 
represented by and 3 ^ are functions of h only as 

the gas is in local thermodynamic and chemical equilibrium and the 
pressure has been assumed constant throughout the stagnation region. 
The quantities Pr g , Re g , and Pe g are the Prandtl, Reynolds, and 
Peclet numbers, respectively, based on conditions immediately behind 
the shock. The parameters T, ly and e are the inverse of the 
Boltzmann number, the Bouguer number, and the radiation cooling 
parameter, respectively. These parameters are fundamental to the 
study of radiation gas dynamics and have been discussed by a number 
of investigators (see, for example, refs. 21 and 22). 

Substituting the above listed nondimens ional quantities into 
equations (2.1)0), (2.42) through (2.1)7), (2.52), and (2-53) yields 
the nondimens ional system governing the flow in the stagnation 
region of a blunt body traveling at hypersonic speeds. 
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f(l) h'(Ti) + Pe 


-1 


^ 1 (h) h(n) 


+ el 


GO- 


(2.55) 


2Re 


3 2 (h) f"(n) 


+ 2f(r,) f ’* ( tj ) - 


f(n) 


2 + a 2 3 5 (h) = 0 (2.56) 


f(o) = 0 


(2.57) 


f’(O) = 0 


(2.58) 


f h* - 1 


(2 09) 


r « ■ ®) ■ 


V 2 ^ 1 - *) 


(2.60) 


h(0) = h = | V/ 2 h 
' w 2 00 w 


(2.61) 


^ =■ 


(2.62) 


where is the value of the nondimens ional Dorodnitsyn variable 

at the shock. This quantity is determined from the expression 


A 


A ’-'O 


A h(q)dq 


(2.63) 


The integral term IpqJ is given by the expression. 
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lGD= -2Kp(n) B(n) 




w 



In two chapters (V and VI) of this paper, it will be convenient 
to express the energy equation in terms of the optical path length as 
independent variable. In both cases the optical properties of the gas 
will be assumed to be independent of wavelength. In this event, the 
energy equation (less the thermal conductivity term) becomes 


f(7) h>(7) + el(7]= 0 ( 2 . 65 ) 

where 



— iR] = - 2 B( t ) 

k(^) 


T 



(2.66) 




+ 2r 


w 


E, 


0* 7 ) / 0 Ta 5(T) % 
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Throughout the remainder of this paper, the bars over the non- 
dimensional variables will be dropped. This should not lead to any 
confusion because only the nondiraensional form of the governing 
equations will be employed. 

C. The Divergence of the Radiant Flux 

The nondimensional form of the divergence of the radiant flux 
is 

* -2e"p(n) B( tj ) 

+ f Q ^ {n) lf 0 A I T A^ - T A (V)|JdV 

n -a < 2 * 6 7) 

+ 2r w ^(^P T A (| 1 ) ) f Q A “x^') E 2 (%> T A (q')JdT)'ldA 

The first term on the right -hand -side of this expression is the local 
emission term which represents the rate at which energy is emitted per 
unit volume of gas at the location tj. The integration over all wave- 
lengths A has been performed for this term with the aid of the 
definition of the Planck mean mass absorption coefficient (see below). 
The second and third terms represent the rate at which radiant energy 
is absorbed per unit volume at the location tj. 

It is the presence of the second and third terms which so greatly 
complicate the radiation problem. These terms are integral expressions. 
In addition, their presence makes it impossible to define a wavelength 
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averaged absorption coefficient by which the wavelength dependence 
might be eliminated. The importance of these terms is indicated by 
the magnitude of the Bouguer number kp which is the ratio of the 
shock standoff distance for a nonradiating shock layer to the photon 
mean free path evaluated at conditions immediately behind the shock. 

The radiation cooling parameter e is a ratio of the rate of 
energy loss per unit area by radiation from both sides of a non- 
aosorbing isothermal layer of gas of thickness A. to the rate at 
which kinetic energy enters the shock layer per unit area of shock 
surface. Alternatively, the parameter e may be interpreted as the 
ratio of the radiationless standoff distance to the decay length where 
the decay length is the length required by an element of gas to lose 
all the energy it possessed upon emerging from the normal shock if 
it loses this energy by radiating (without reabsorbing) at a con- 
stant rate. This parameter modifies the entire radiation term 
thus, acts as a measure of the relative efficiency of radiation 
compared to convection as energy transport mechanisms within the 
shock layer. ' In addition, the surface reflectivity and the 

wavelength dependence of the absorption coefficient influence the 
character of the radiation terms and will be considered as parameters 
in this study. 

Most investigators who have studied problems in which a term 
similar to ijjJ appears have assumed that the gas and its sur- 
roundings are gray, that is the optical properties are independent 
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of wavelength. This allows the integration over frequency to be 
performed analytically. Accurate results can be achieved in the 
two extreme cases of optically thin « lj and optically thick 

» lj gases. When the gas is optically thin at all wavelengths, 
the gray absorption coefficient is correctly given by the Planck 
mean mass absorption coefficient 


K p ■ J K x ®x 

oT ° o 


dA 


( 2 . 68 ) 


Where k.^ is the monochromatic mass absorption coefficient and 
the weighting function is the Planck black-body function. 

When the gas is optically thick at all wavelengths, the gray absorp- 
tion coefficient, in the interior of the gas, is correctly given 
by the Rosseland mean mass absorption coefficient 



-i 4B * dX 

x 5T tt 


(2.69) 


Near a radiation boundary or in regions of rapid (with respect to 
the optical path length) variations in thermodynamic properties the 
Rosseland mean is not valid. At intermediate values of optical depth, 
no single mean absorption coefficient, which depends only on local 
thermodynamic conditions can be defined. In fact, as has been pointed 
out by Krook (ref. 23) it would be necessary to define an infinite 
number of such mean coefficients. This, of course, does not preclude 
the possibility of defining approximate mean coefficients under 


these conditions. 
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Stone (ref. 24) introduced a model in which the monochromatic 
absorption coefficient was a step function of frequency with the size 
of the steps independent of the geometry or thermodynamics of the 
system. By means of this method, the integral over all wavelengths 
is reduced to a finite series. Carrier and Averrett (ref. 25) con- 
sidered an absorption coefficient with only two steps, one of which 
was very much larger than the other. Both of the papers noted above 
were concerned with Milne's problem of a stellar atmosphere in 
radiative equilibrium. Lick (ref. 11) and later Grief (ref. 12) 
studied the problem of one -dimensional energy transfer between two 
walls separated by a radiating and conducting gas. A picket fence 
model, which is a specialization of the step function model, for the 
absorption coefficient was used. Krook (ref. 26) derived expressions 
by means of the P-L-K perturbation procedure for a slightly nongray 
gas. The solution represents a perturbation to the gray gas solution. 
Rhyming (ref. 27) considered wave propagation in a simple dissociating 
flow of a radiating gas where the absorption coefficient was given as 
a Gaussian function of the frequency. 

However, even with the above simple models for the absorption 
coefficient, the term retains an integral character and the 

solution to the set of equation is still very difficult to obtain. 
Numerical procedures are extremely tedious. For example, it was 
pointed out in reference 9 that the time to obtain solutions on the 
IBM 7090 to a similar (though not identic!) set of equations with the 
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gray gas assumption ranged from 20 minutes to 5 hours. As a result 
of this difficulty, several approximate analytical methods have been 
derived in order to reduce this term to purely differential form. 

One such technique is the Milne -Eddington approximation (ref. 28), 
the derivation of which has been based on physical considerations, 
but which may also be thought of as a substitute kerned approximation 
(ref. 29)* The integral terms can then be eliminated by means of a 
double differentiation (for a gray gas only). Of course, this 
increases the order of the differential equation by two. This 
technique has been used by a number of authors in the study of the 
dynamics of radiating gases (see for example, refs. 30 and 3l)** 

Barbier (ref. 28) introduced the method of expanding the source 
function in a Taylor series about the zero of the argument of the 
exponential integral kemal, E^|i\(q) - t-^( Q 1 )|1 • Because the kemal 
function has a logarithmic singularity at the zero of its argument, 
the integral over the first few terms of the series should provide a 
good approximation. The resulting integrals can then be evaluated 
analytically and the equation becomes purely differential in character. 
Yoshikawa and Chapman (ref. 8), Thomas (ref. 33 ) > and Viskanta (ref. 3*0 
all used the method of Barbier to different degrees of approximation. 


Traugott (ref. 32) has introduced a "method of moments" in order 
to reduce the integral term to differential form. This method may be 
taken to any degree of approximation desired (not without a considerable 
sacrifice in simplicity however). The first approximation is identical 
to the Milne -Eddington approximation. Traugott' s higher approximations 
can also be obtained by a substitute kernal method similar to that of 
Krook (ref. 29). 
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When more than the constant term in the Taylor series is retained, 
it may become necessary to introduce additional boundary conditions. 
In fact, it may not be enough to merely specif^ a new condition, it 
may be necessary to modify the existing conditions as well so that 
in the limit as the parameter H = ki^crT^ (for example) tends to 
zero, the solution will approach the proper pure radiation solution.* 
Apparently, this point was overlooked by Viskanta. In his paper, 
Viskanta blamed the failure of his pure radiation solutions for a 
finite optical thickness to exhibit a temperature jump at the wall 
on a premature truncation of the Taylor series expansion of the source 
function B(t) . Actually, this failure was a result of improperly 
specifying the boundary conditions. 

The diffusion approximation for optically thick gases has been 
used extensively in astrophysics and gas dynamics. Probstein 
(ref. 35) has shown how to extend the usefulness of this approxi- 
mation to gas layers of finite optical thickness by means of radiation 
slip boundary conditions. It is not at all clear, however, that these 
slip conditions can be used in the problem of this paper because of 
the presence of the convection term pwdh/dz. 


*This parameter, which appears in the literature concerning 
energy transport by radiation and conduction, represents the relative 
importance of conduction compared to radiation. When N tends to 
zero radiation is the dominant mode of energy transport. 
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The optically thin approximation of hot cases, in which absorp- 
tion is neglected in comparison to emission has also been used exten- 
sively in gas dynamics. As Thomas (ref. 55) has pointed out, this 
approximation is not valid in those portions of the gas which are 
considerably cooler than the remainder of the gas. 

In this paper, the integral term lQ)J will be reduced to alge- 
braic or differential form through the use of various approximations 
similar to those described above. The manner in which this is to be 
accomplished will depend on the order of magnitude of the parameters 
e and k p and will be discussed in detail in the next four chapters. 
Whenever possible, the gas will be treated as nongray. 

D. The Inviscid Shock Layer 

As was pointed out in chapter I, the studies of this paper will 
be concerned only with those cases for which the thicknesses of the 
wall boundary layers due to the presence of viscosity and thermal 
conductivity are very much less than the shock standoff distance. For 
a nonradiating gas , the shock layer can be separated into an outer 
inviscid and nonheat conducting region and an inner viscous and heat 
conducting region or boundary layer. Considerable simplication will 
result if a similar separation can be achieved in the case of a 
radiating gas. As will be shown, such a separation can be obtained 
when the boundary is either optically thin or optically thick. Only 
the former situation will be considered herein. The method of 



separation follows the procedures delineated by Van Dyke (ref. 36). 
Mathematical details are presented in appendix A. 

It is shown in the appendix that the significant parameter which 
determines the extent of the boundary layer is the inverse square root 
of the Peclet number, Pe The zero-order in Pe*"^"” system of 

equations which governs the flow in the inviscid region is 


f 0 h) + gI ( £t] = 0 

(2. TO) 

f o ( '> - f 


2 + a 2 3t 5 (h) . 0 

(2.71) 

L 



r o<°) ■ 

. 0 

(2.72) 

f 

0 

■w- 

= 1 

(2.73) 


rA. A 



f ‘ a —I 

-i. 

a 

(2.7 J 0 

o^ xl 

, R J 

i 

X 

h 

0 

N- 

= 1 

(2.73) 


The dependent variables f^q) a nd h Q (q) are the asymptotic values 

of f( r) ) , the nondimens ional stream function, and h(q), the non- 

-1 /2 

dimensional enthalpy, respectively as Pe approaches zero. 
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In the boundary layer , the zero -order system of equations is 


% 

' 3 l (i o )i o (l) j' + 6 0 (E)iod) + 0 

(2.76) 

► 

Ws, i„ e"( S)”|' + SB„( e)B*(£) -f E AU)1 a + - 0 

L J ' L' j ' v y 

(2-77) 


g o (0) = 0 

(2.73) 


g^(o) = 0 

(2.79) 


lira g;U) = f^(0) 
1 00 

( 2 . 80 ) 


lira i (0) = h 
t o N w 

(2.31) 


yu = h o (o) 

( 2 . 82 ) 


The independent variable i is the "stretched" boundary layer coor- 
dinate defined by the relation 

I -Pe 1/? tj (2.33) 

The dependent variables g( l) and i(|) are defined by the expressions 

i(t) = h(t)) (2.34) 

g'(0 = f’(n) (2.35) 

in the boundary layer as Pe - ^^ approaches zero. 
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The integrals which appear in the second of these expressions are 
definite integrals. Consequently, the system of equations governing 
the flow in the boundary layer is a purely differential, system. 

It must be realized that expressions (2.86) and (2.87) are 
restricted to the case of an optically thin boundary layer. It is 
only in this case , and the case for which the boundary layer is 
optically thick, that a complete separation between the inviscid 
region and the boundary layer can be achieved. At intermediate 
values of optical depth, the integral term I q Qi 1 is a function 
of the enthalpy distribution in the boundary layer in addition to 
being a function of the enthalpy distribution in the inviscid region 
so that the equations in the inviscid region and the boundary layer 
are coupled. The influence on the inviscid region of radiation from 
an optically thick boundary layer cannot be neglected. However, most 
of this radiation originates at the outer edge of the boundary layer. 
The boundary layer solution in this region is constrained by matching 
conditions to approach asymptotically the value of the inviscid 
solution at the wall. Hence, the radiation contribution to the 
inviscid region from the boundary layer can be obtained from the 
inviscid solution at the wall, leaving the inviscid solution uncoupled 
from the boundary layer solution. 

This restriction to an optically thin boundary layer is not so 
severe as it might first appear. This is because the optical thick- 
ness of a boundary layer in which the absorption coefficient is the 
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same order of magnitude as it is for shock heated air will not exceed 
about 0.1 at any altitude and velocity (up to km/sec) for a shock 
radius of 1 meter or less. In fact, the optical thickness of the 
boundary layer ’..'ill be less than 0.1 at that point of the trajectory 
of a Martian or Lunar return vehicle with a shock radius of about 
1 meter for which heating is a maximum even if the absorption coeffi- 
cient in the boundary layer is 100 times that of shock heated air. 

That this should be the case is not so difficult to see when it is 
realized that both the optical path length and the boundary layer 
thickness decrease rapidly with decreasing altitude. Thus, at low 
altitudes where the optica] path length is small and the shock layer 
may be optically thick, the boundary layer thickness is very small. 

For larger objects, the boundary layer need not be optically thin at 
the lower altitudes because the boundary layer thickness depends on 
the size of the object while the optical path length does not. 

These conclusions regarding the optical thickness of the boundary 
layer generally concur with the observations of Fay, Moffatt, and 
Probstein (ref. 10). Henceforth, the discussions of this paper will 
be limited to the case of an optically thin boundary layer and 
radiation from this boundary layer will be considered to have no effect 
on the solution in the inviscid region of the shock layer. 

If the inviscid system of equations ( 2.70 ) through (2.75) is 
solved for the nonradiating case (e = 0) along with condition (2.6y) 
one finds that the ratio of the shock standoff distance to the shock 


radius is given by the expression 
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x 

1 + V^vl - X ) 


(2.83) 


Hayes (ref. 37) obtained the same result when the shock and body 
surfaces near the stagnation point are concentric spheres. When 
the shock and bodv surfaces are not concentric (i.e.. R ± R.. + A) 

” * ' S ' 1M 

condition (2.83) is still approximately true over a wide range of 
body shapes (see, for example, refs. 38 and 39)- With this result 


a . (2.0. 

1 + Y2X(1- X) 

This value for a, the constant appearing in the momentum equa- 
tion (2.7l), will be used throughout the remainder of this 
investigation. 

E. Thermodynamic and Optical Property 
Correlations 

In order to achieve meaningful results , an attempt has been 
made in this paper to use simple yet physically reasonable approxi- 
mations to the thermodynamic and optical properties of high tempera- 
ture gases. In particular, coorelation formulas were derived from 
the existing store of information about equilibrium air. The thermo- 
dynamic properties were obtained from reference 40, for temperature 
up to 100,000° K and pressures from 10 to 10 2 times atmospheric. 

The optical properties were obtained from a variety of sources which 


will be noted later. 
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It was noted from the data of reference 40, that both the density 
and temperature could be approximately represented by functions 
separable in the variables pressure and enthalpy. More specifically 
in the form /PqJ* f ^/^o) ' Plots of the factions f ^h^RT^ for 


the density and temperature at various pressure levels are presented 
in figures 2.1 and 2.2, respectively. It is apparent from these plots 
that the density and temperature can adequately be represented by the 
expressions 



(2.90) 


(2.91) 


A number of investigators (see, for example, refs. 41-45) have 
calculated the radiant properties of equilibrium air for temperatures 
up to 25,000° K and for densities from 10 ^ to 10^ Amagats. Because 
of the extremely complex nature of these calculations, the many 
physical processes which produce radiation, and the uncertain knowl- 
edge of cross sections and transition probabilities the scatter among 
the various calculations is often quite large. Some of the results 
for the Planck mean mass absorption coefficients are presented in 
figure 2 . 5 . 

A correlation formula can be obtained from figure 2.5 by 
approximating the curves of Log^ pKp versus Log^T with straight 
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(°J,H/^V ,T fx, T 


Figure 2.1.- Density correlation. 
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Figure 2.2.- Temperature correlation. 




Figure 2.3* - Planck mean mass absorption coefficient correlation. 




Figure 2.3*- Continued 
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lines. The resulting formula is 


pK P 


= 7-94 X 10 


-26 


y 


' 3 ' 25 6.0-0. 5 Log ~ -1 

T 10P o , cm 


(2.92) 


ITTVs A r« -5 r. -v-» /-> 4 . A /»' ,1 n V.1 -.*■ n ^ «n, ^ 4, ^ . 4 „ iv A ^ 1 1 f 

1UJ.O -to iwu o i>ax uxtuxaixj ^uiivciucao xuiiii iui uoc xii uic tanviuatiuili 

of this paper. It is much more desirable to express the Planck mean 
mass absorption coefficient K p in terms of the pressure and enthalpy. 
This was done by cross -plotting the logarithm of the absorption 
coefficient data shown in figure 2.3 against the logatithm of the 
temperature at constant pressure. Straight lines were then fitted 
to the resulting curves. Finally, the correlation formulas (2.90) and 
(2.91) for density and temperature were used to obtain the formula 


10 


Kp = 1.39 x ia*~[J2. 


^.34-0.44Log 10 ^ 


(“o 


3 • 55 -0 . 24Log 


10Pr 


Z h 


, cm /gm 


(2.95) 


This formula is valid for temperatures up to 20,000° K at the higher 
pressures ^"p/p q = 10 ^ to lO^ j and to somewhat lower temperatures 
at the lower pressures ^ for example, when p/p Q = 10 the maximum 
temperature at which the formula is valid is 15,000° k| . 

In addition to the Planck mean, it is necessary to know the 
spectral variation of the mass absorption coefficient. Results of 
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some typical calculations of the monochromatic mass absorption 
coefficient plotted as functions of wavelength for constant tempera- 
ture and pressure are presented in figure 2.4. No attempt was made 
to correlate these data.* 

At temperatures above about 20,000° K, the information about the 
radiative properties is not so comprehensive. Most of what exists 
consists of Planck and or Rosseland mean absorption coefficients for 
continuum radiation. Line radiation is neglected. At these high 
temperatures , the radiation consists of spectral lines of the various 
ions which may be appreciable Stark -broadened at high electron 
densities , and a continuum due to free bound and free-free transitions 
of electrons in collisons with the ions. Since the integrated line 
emission is proportional to the ion density while the continuum 
emission is proportional to the product of ion and electron densities, 
the ratio of the latter to the former increases with increasing 
density. Thus, at the higher density levels, the continuum calculations 
may be adequate. 

* 

It was noted from the results presented in references 45 and 46 
that the functional form of the monochromatic mass absorption coeffi- 
cient is approximately 

tc^(h) = \ Kj(h) gj(A) (2.94) 


where the subscript refers to the 1th radiating species. 






Figure 2.4.- Continued 





results of calculations of the Planck mean volume absorption coefficient for 
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CHAPTER III. THE SMALL PERTURBATION SOLUTION 


A. The Conventional Method 

As was pointed out in chapter I, there is a flight regime of 
considerable importance in which the radiation cooling parameter e 
is very much less than unity. In this regime, the energy transferred 
by radiation is small compared to the influx of kinetic energy across 
the bow shock, and it would be reasonable to expect the flow prop- 
erties to be only slightly perturbed from the radiationless case. 
Lunev and Murzinov (ref. 4) and Goulard (ref. 5) took advantage of 
this and developed what amounted to first order perturbation solutions 
of the temperature distribution in the inviscid region of an 
transparent, gray gas layer. In both these papers, simplifying 
assumptions concerning the gas properties and flow model have been 
included. 

In this section, the perturbation solutions will be generalized 
to include nongray gases with arbitrary thermodynamic and optical 
properties. These solutions will not be limited to shock layers of 
small optical thickness. Also, the solutions will be extended to 
second order. As will be shown, the second-order solutions can be 
quite important when the absorption coefficient varies rapidly with 
temperature. 

The integrodifferential system which governs the flow in the 
inviscid region of the shock layer is 
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f(n) h'(T,) + elQ]- 0 

2f(n) f"( n ) - [fh)] + a 2 h(n) = o 


( 5 * 1 ) 


( 5 . 2 ) 


f(o) = 0 


( 5 - 5 ) 


f (^)= 1 

f (\)= 1 ( 5 - 5 ) 


f 



s/a'i _ « 

X W " VsxU -X) 


( 3 . 6 ) 


Here f(n) and. h(rj) are the nondimens ional stream function and 
enthalpy, respectively. The quantity tj^ is the value of the 
Dorodnitsyn coordinate at the location of the shock. The constant 
a can be expressed in terms of X, the density ratio across the 
shock, through expression (2.89). When the radiation cooling 
parameter € is very small, the integral term in equation (5«l) 
becomes of only secondary significance throughout most of the domain 
of the problem. Neglecting the integral term ifjl”) reduces the 
problem to one in which radiation does not play a part. If, as 
expected when e is small, the presence of radiation only slightly 


With the obvious exception of the region ^ 88 0 where f(rj) 85 0. 
The difficulties presented by this exception will be discussed later. 
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influences the solution one can, to reasonable accuracy, evaluate 
I [jl^j usin S the radiationless solution for h so that equation (3.1) 
becomes purely differential. Thus, when the small perturbation 
procedure (which roughly proceeds in the manner outlined above) is 
applied to this problem, the integrodifferential system is simplified 
to a purely differential system. In addition, as a result of the 
nature of the lowest order solution for the enthalpy distribution, 
the two differential equations become uncoupled and can be solved 
independently. Hence, it becomes possible to obtain analytic solutions 
to any order of approximation to the flow in the inviscid region of 
the shock layer. Details of the derivation of these solutions are 
presented in appendix B. 

The zero-order, or radiationless, solution is simlpy 

h 0 ( T ‘) " 1 (^* 7 ) 

f Q (n) = (1 - a)t] 2 + at} (5.8) 

The first-order solution, which represents the effect of radiation 
assuming that the emissive power of the gas is independent of 
temperature, is 


M 1 * * f 


(l - &)x^ + 


( 3 .; 3 ) 


ax 
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\ (x) 

ftl " a)x + aj 


dx 


+ 


2 

n 



2(1 - a)x + a 
x 2 [(l - a)x + 0 


(x) dx\ 


( 5 . 10 ) 


Here x is a dummy variable of integration. The quantity I Q QfJ is 
given by the formula 



( 5 . 11 ) 


The notation has been simplified somewhat in this expression by 

omitting the arguement h in the terms k and B, and by 

o A A 

introducing the quantities 


k x = ‘V ^ 


( 3 - 12 ) 



Also 

^(l) = - | a 2 h^Tj) 


( 3 . 15 ) 


(3-lM 


The second-order solution takes into account the change in 
gas properties with changes in enthalpy. This solution is 
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h 2 (n) - ^ ijg + f l fl ^ . fffigjl 

\ 1®(*) J 


dx (3.15) 


f 2 (o) » - (l - a) 2 t] 2 


pT * 2 (x) 


1 r- r 'l 

^[ 2 (l-a) n + g J — 

J o [(l - a)x + a] 


dx (3»lo) 


1 Tj 2 f' 1 ^(1 - a)x + a 

2 ^ tj x 2 jjl - a)x + aj 


2 2 


4> 0 (x) dx 


f'here 


J3 - V’’’ j7 [** {‘bCM 1 - iD 


+ f 1 • r 0?> ) E 2[ k A l]} - 2 "a \ 


> dA 


• pi 

k vK E A h l<*> ®i|>X 1 11 - x l> (5-^T) 


+ S> { e iJ?a (1 • tO 7 h l (x) 4:1 

1 

+ (* - "o^V] P + \ E 2 [V>] ,' dX 
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Here the arguement h is omitted in the terms A. and B, and the 

o A A 

quantity r.. is defined by the expression 
X A 


r n = 2r 

h ” 


• r 


4 t i U) dx 


+ \ \ E 2 (^ I 0 h l (x)d!C + \ \ \ E 2 (>) 


Also 


O.lo) 


® 2 (n) = -f^n) fj( n ) + | ff^i) 


| &2 hg(n) 


( 5 - 19 ) 


The quantities , tj. , and are given by the formulas 

o x 2 


% 


= 1 


( 5 - 20 ) 


! ^1 * x (x) 

^1 ^ [(1 - a)x + a] 2 


( 5 - 21 ) 



(1 - a)r^ + 


1 

2 




^(x) 

- a)x + a]^ 


dx 


( 5 - 22 ) 


It can be seen upon inspection of relation (5.17) that a large 
value of the rate of change of the Planck mean absorption coefficient 
with enthalpy will lead to large values of I^JjQ . Thus, it is clear 
that at shock temperatures of less than about 30,000° K, for which the 
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absorption coefficient does vary rapidly with enthalpy, the second- 
order solutions can become more important to the overall solution 
than their order in e might at first indicate. 

B. The P - L - K Method 

As can be seen from an inspection of the expressions (3.9) and 
(3.15) the first order solution for the enthalpy distribution has a 
logarithmic singularity at the point rj = 0 and the second-order 
solution has a singularity of greater strength at this point. As a 
consequence, the assumed expansion diverges as the origin is approached 
and the small perturbation solution is not uniformly valid throughout 
the domain of the problem. This divergence can lead to serious errors 
in the calculation of the radiant heat flux to the wall because those 
regions close to the wall, in which the largest errors occur, are 
given the most weight in the calculation. This is particularly true 
for shock layers which are not optically thin. Additional difficulties 
are encountered when attempting to specify the proper outer boundary 
conditions for the viscous boundary layer equations. In classical 
boundary layer theory, the outer boundary conditions are obtained 
from the values of the outer (or inviscid) solution at the wall 
(t = 0 in this problem). Because of the divergence of the outer 
solution, no finite value exists at q = 0. 
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In this section, the Poincare-Lighthill-Kuo perturbation of 
coordinate procedure (ref. 47 ) is used to obtain a solution which 
is uniformly valid over the domain of the problem. The details of 
the application of this method to the problem of this paper are 
presented in appendix B. This method utilizes a coordinate trans- 
formation in the form of a perturbation expansion of the coordinate 
to remove the singularity (which caused the divergence of the con- 
ventional solution ) from rj = 0 to a small negative value of r) 
which lies outside the domain of the problem. The P-L-K expansions 
are 

T) = X + €T1* (x) + ... (3«23) 

h(Tj;e) = h Q (x) + el^ (x) + ... (3*24) 

f(a; e ) = f* (x) + €f* (ti) + ... (3.25) 

where x is the coordinate in the transformed plane, and the 

superscript * has been used to differentiate between the coeffi- 
cients in the P-L-K expansion and the coefficients in the conventional 
expansion. Pritulo (ref. 48) has derived a general relation between 
the P-L-K and conventional coefficients. Adapted to this problem, 
the relationships become 

Variously called the P-L-K method, the P-L method, Lighthill's 
technique, the method of strained coordinates, and the method of 
perturbation of coordinates. 
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* 

h 

0 

(x) 

- h o w 

( 3 . 26 ) 

* 

h i 

(x) 

N, 

X 

II 

(3.27) 

* 

f 

0 

(x) 

= f o 

( 3 - 28 ) 

f* (*> = f x 

(x) 

+ n* (x) f^ (x) 

(3-29) 

\ (x) = 

: - h 2 

, (x)/h£ (x) 

(3.30) 


The second -order terra (x) introduces the effects of variable 
thermodynamic and optical properties, so it is apparent that these 
effects are contained in the first-order P-L- K solution. 

A comparison of the P-L-K and conventional perturbation solutions 
for the enthalpy distribution for a constant density, transparent 
shock layer is presented in figure 3*1. The divergent character of 
the conventional solutions is apparent. Also shown on this figure 
is the exact analytic solution which can be obtained in this simple 
case. The formula for this exact solution is 



where y (the exponent in the correlation formula k^B = h^) was 
taken to be 6 and the constant a (which appears in the momentum 





Exact solution (eq. (3.31)) 
Zero-order perturbation solution 
First-order perturbation solution 
Second- order perturbation solution 


00 


O 


P-L-K solution 



72 


eq. (3*2)) was taken to be 0.5. The good agreement between the 
P-L-K solution and the exact solution indicates that the accuracy 
of the P-L-K solution is probably second-order in the radiation 
cooling parameter € throughout the domain except in the immediate 
neighborhood of the wall. It is clear that quantities such as the 
radiant heat flux at the wall, which depend upon an integration 
over the enthalpy distribution, will be considerably more accurate 
if the P-L-K solution rather than the conventional perturbation 
solution is used. 

It should be noticed that the P-L-K solution does not lead to 
zero enthalpy at the wall as the exact transparent solution does. 

The reason for this disparity can be found in the fact that the 
coordinate stretching displaces the boundary with regard to both the 
energy and momentum equations but not by a uniform amount. Thus, a 
physical interpretation of the first order P-L-K solution is that the 
normal velocity of the flow at the boundary for the energy equation 
is not quite zero, and a particle approaching this boundary will 
reach it in a finite time before losing all its energy by radiation. 

It can be shown that since the expected error in the Dorodnitsyn 
coordinate tj in terms of the stretched coordinate x is order e 
and since the gradients in h^ (x) are very large in the vicinity 
of the wall, the difference between the P-L-K and exact solutions at 
the wall lie within expected limits. Convergence to the correct 
solution should be attained with the addition of higher order terms 
to the expansion of h and tj. 
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C. The Method of Matched Asymptotic Expansions 


Van Dyke (ref. 36) has pointed out that the method of matched 
asymptotic expansions is applicable whenever the P-L-K method can be 
used. Thus, it would be interesting to formulate the solution when 
radiation is a snail perturbation using the method of* matched asyrup— 
totic expansions. Use of this method implies that the domain of the 
problem can be divided into at least two regions in which the governing 
equations take on different asymptotic forms. There must also be some 
overlap between adjacent regions so that a smooth transition between 
solutions valid in these adjacent regions can be affected. In the 
problem of this chapter, the regions are the "outer" region in which 
the conventional perturbation solutions are valid and the "inner" 
region in the vicinity of the wall at q = 0. The equations which 
describe the conditions in the outer region are simply the system (3-1) 
to (3*6). In order to obtain the "boundary layer" form of these 
equations, it is necessary to stretch the coordinate tj in the 
vicinity of the wall. This stretching takes the nonlinear form 


F(0=[f(Ti)] e (3-32) 

F’U) = f f (n) 


where I is the stretched boundary layer coordinate and F( i) is 
the velocity function written in terms of t* It follows from 
above that 


at f F(l) 
a n = £ KrO 


(3.33) 
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and the energy and momentum equations, respectively, take the forms 


V 

fU)h’(i) + i[t) =0 

(3-34) 


2F(4)F"U) -[f'U)] + a 2 H(4) = 0 

(3-35) 

4k 

vhere 



in* 

II 

(3.36) 


and 



~r--> , , 




(3.37) 


One boundary condition is available, that is 

F(0)=0 ( 3 . 38 ) 

The remaining two constants of integration can be obtained by 
matching the inner and outer solutions according to the matching 
principle put forth in reference 36 . 

The boundary layer system is seen to be quite complex. The 
energy and momentum equations remain coupled so that it is necessary 
to obtain a simultaneous solution to the two equations. Thus, as 
is often the case when the P-L-K method can be applied, its 
application is much simpler them the method of matched asymptotic 
expansions. 


* 
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D. Results and Discussion 

The formulas derived in the preceeding sections of this chapter 
were programed for numerical computation on the IBM 7094 electronic 
digital computer. The value of X, the density ratio across the 
normal shock, was fixed at a constant value of 0.06 for the calcula- 
tions reported on in this and subsequent chapters. This choice is 
justified because X varies but little with altitude and velocity 
and the effects of this variation on the stagnation solutions are 
slight. The value X = 0.06 is typical for hypervelocity flight in 
the atmosphere of the earth. 

The numerical calculations indicate that the enthalpy is a 
double valued function of the Dorodnitsyn coordinate i] in the 
vicinity of the shock for large values of the Bouguer number. 

An examination of the governing equations failed to show the presence 
of any singularities which might adversely influence the solution In 
this region when kp is large and e small. On the other hand, 
the results of numerical calculations with varying mesh size seemed 
to rule out the possibility that the doubled valued behavior can 
be attributed solely to numerical inaccuracies. Consequently, it is 
suspected that the difficulty results from truncation of the 
perturbation expansion and that inclusion of higher order terms would 
either eliminate the problem or increase the value of kp at which 
it first appears. For truncation after the second order term the 
conditions for validity of the solution is ek„ < 1. 
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Gray gas results 

Shock layer enthalpy distributions for a gray gas with differing 
values of the radiation cooling parameter €, the Bouguer number kp, 
the variation with enthalpy of the Planck mean mass absorption coeffi- 
cient k , and the reflectivity of the body surface r are presented 

ir W 

in figures 3*2 to 3-^- While the gray gas assumption may not be 
realistic for most gases of interest, its use is felt to be justified 
in the study of the above listed parameters for two reasons. First, 
the highly complex and varied spectral structure of absorption 
coefficients makes a general parametric study of nongray gases imprac- 
tical. Second, experience with nongray calculations indicates that 
the qualitative dependence of the gray results on the various parame- 
ters will carry over to most nongray cases. 

The decrease in enthalpy level with increasing e is illustrated 
in figures 3*2a to 3 -2c. These results indicate that the loss of 
energy from the shock layer by radiation (i.e., radiation cooling) can 
produce a noticeable drop in enthalpy for values of e as small as 
0.01. The dependence of the enthalpy distribution on the Bouguer 
number (hence, optical thickness) is also shown in these figures. 

As expected an increase in the Bouguer number (or optical thickness) 
inhibits shock layer cooling and leads to higher values of enthalpy 


near the wall. 









1 (b) kp = 0.3. 


Figure 3-3*- Continued. 
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Figure 3*^»- Effect of surface reflectivity on the shock layer 

enthalpy distribution. 
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The variation of the enthalpy distribution with Kp (the 
enthalpy variation of the Planck mean mass absorption coefficient) 
for several values of the Bouguer number kp is shown in figures 5* 3a 
to 3* 3c. These effects are most noticeable for optically thin shock 
layers ^kp « l.oj and tend to vanish as the optical thickness 
increases. In a transparent layer, the rate of emission of radiant 
energy is proportional to the Planck mean mass absorption coefficient 
Kp. Thus, gases with small values of Kp (which mean larger values 
of Kp when the nondimensional enthalpy is less than l) will be 
cooled more than gases with large values of Kp. As the optical 
thickness increases smaller fep still implies greater emission rates 
but it also means greater absorption and more radiant energy available 
for absorption. The process of absorption tends to counteract the 
differences in emission rates due to differences in Kp. Finally, 
when radiation equilibrium is reached (this state Is achieved in the 
interior of optically thick regions) the energy of the particle is 
independent of its optical properties. Of course, in those regions 
optically close to the shock and the wall the amount of radiant energy 
available for absorption is not so great as in the interior of the 
shock layer and particles in these regions cannot approach the state 
of radiation equilibrium (except in a region optically close to a 
perfectly reflecting surface). Thus, the enthalpy distribution 
remains dependent on the value of Kp near the shock and the wall. 
This dependence of Kp is surpressed near the shock where h is 
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almost 1 because the values of ttp are nearly the same di spite the 

differences in ft . 

P 

The effect of surface reflectivity, r , on the shock layer enthalpy 
distribution is shown in figures to 3 .4c . If the shock layer gas 

is transparent (i.e., the gas does not absorb) surface reflectivity 
has no effect on the enthalpy distribution because all photons emitted 
by the layer escape. Whether or not a photon is absorbed or reflected 
by the wall is of no consequence. As the optical thickness of the 
layer increases the chance of capture of a photon by absorption in 
the layer is increased. If the surface reflectivity is increased also, 
the probability of capture is increased still further because many 
photons which might otherwise have escaped into the wall are reflected 
back into the layer and are once again subject to capture there. 
Consequently, the enthalpy level is higher near a reflecting wall 
than it would be near a nonreflecting wall . 

It can be concluded from the above, that use of a reflecting 
surface will not reduce the radiant heat transfer rate from the gas 
to the wall by the factor 1 - r w (unless, of cource, the gas is 
transparent) but will reduce it by some smaller fraction. This is 
because the radiant heat flux incident on the wall is larger when the 
wall is reflecting as a result of the higher enthalpy level. In 
addition, the rate of heat transferred to the wall by conduction 
will be greater, also because of the higher enthalpy level. Of 
course, increasing the surface reflectivity always decreases 
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the total heat transfer rate to the wall because the higher enthalpy 
level must lead to an increased loss of energy by radiation through 
the shock in the upstream direction and by convection in a lateral 
direction away from the stagnation point. If the energy balance is 
to be maintained, the rate of heat transferred to the wall must be 
reduced. 

The effects of variations in the parameters €, kp, ftp, and 

on the rate of radiant heat transfer to the wall (normalized by the 

1 3 R 

energy influx to the shock layer, ^ p Vr ) q are shown in figures 
3-5 to 3 * 7 * The rate of radiant heat transfer to the stagnation 
point was calculated with the formula 



el 


- rj J A KpU) B(n) E 2 ^k p T^jdTj 


where the optical thickness 


kp t is given by 


(3-39) 


fcp t(tj) 



(3. to) 


The dashed curves in figure 3*9 indicate the "no decay limits" 
for various values of the Bouguer number. These limiting curves were 
computed by assuming the shock layer to be isenthalpic so that 
Kp(q) = B(q) s 1. Thus, the no decay limit curves are given by the 


formula 
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Figure 3.6.- Effect of the enthalpy variation of the absorption 
coefficient on the rate of radiant heat transfer to the 
stagnation point. 
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where E^kpj is the exponential integral function of third -order. 
This no decay approximation is often used to predict the rate of 
radiant heat transfer when radiation effects are small. Use of this 
approximation always gives an upper bound to the true value of q^. 

A study of figure 3.5 indicates that the no decay limit curve is 
least accurate in predicting the rate of radiant heat transfer in 
the transparent case k p = 0 . This result is expected because the 
enthalpy distribution for the transparent case is the most perturbed 
from an isenthalpic state. Results presented in this figure also 
indicate the importance of absorption (as characterized by the Bouguer 
number kp in reducing the rate of radiant heat transfer from the 
shock layer to the wall. 

The results presented in figure 3*6 indicate the differences in 
Kp, the enthalpy variation of the Planck mean mass absorption coeffi- 
cient, are most important when the optical thickness of the shock 
layer is small. 

Here the radiant heat transfer to the wall is greatest for the 
smallest value of Kp. This, of course, supplements the observation 
(from fig. 3* 3a) that radiation cooling is greatest for gases in 
which Kp is least. The differences in radiant heat transfer to 
the wall brought about by differences in the value of Kp tend to 
vanish as the optical thickness of the layer increases. 

The reduction in radiant heat transfer to the wall due to surface 

reflectivity is shown in figure 3*7* When the shock layer is trans- 

R 

parent, the rate of radiant heat transferred a is in the ratio 
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1 " r w * However, as the optical thickness of the shock layer increases, 
the ratio becomes somewhat greater than 1 - r^ as predicted in an 
earlier discussion of this section. 

The effect of the parameters e, k^, and r^ on the shock standoff 
distance is shown in figures 3.8 and 3.9. The quantity A is the 
ratio of the shock standoff distance in a radiating shock layer to 
that in a nonradiating (or adiabatic) shock layer at the same flight 
conditions, and was computed with the formula 

A= f A h(>])dr] ( 3 * 33 ) 

0 

The results shown in figures 3.8 and 3.9 indicate, as expected, 
that a decrease in enthalpy level (with the consequent increase in 
density level) in a shock layer leads to a reduction in shock stand- 
off distance. 

Nongray results 

It can be seen from figuies 2 .k that the absorption coefficient 
of high temperature air depends strongly on wavelength. This is true 
of all other gases as well. Consequently, the assumption that the 
gas is gray (i.e., that the optical properties of the gas are inde- 
pendent of wavelength) is poor indeed, and has been resorted to so 
frequently in the literature only because of the resulting relative 
simplicity. Fortunately, the small perturbation solution derived 
in this chapter overcomes these difficulties by reducing the 
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absorption integrals in the divergence of the radiant flux to a form 
amenable to direct evaluation. Thus, one need only perform an 
integration over a known, albeit complicated, function of wavelength. 

In view of the current uncertainties, with regard to spectral dis- 
tributions of gaseous absorption coefficients, it was decided to use 
a simplified model for the absorption coefficient of air. Consequently, 
the step function model shown in figure 3. 10 was chosen for use in 
calculations to be reported on herein. The height and width of the 
steps were chosen so that the simple step function model provides an 
adequate representation of the absorption coefficient of air at a 
temperature of about 15,000° K as predicted by Nardone et al. (ref. 1*5) 
and so that the Planck mean absorption coefficient of both distri- 
butions are equal. The relative heights of the nine steps located 
at wavelengths less than 0.113 microns were chosen to be independent 
of enthalpy while the tenth step which covers the wavelength interval 
(0.113, ") was chosen to vary as the 1.28 power of the enthalpy. The 
relative heights shown in figure 3*10 are for h = 1, where h is 
the nondimens ional enthalpy. The enthalpy variation of the step 
heights listed above is consistent with the condition that the Planck 
mean mass absorption coefficient is proportional to the fourth power 
of the enthalpy. 

Shock layer enthalpy distributions were calculated using the 
nongray absorption coefficient model for various values of the 
Bouguer number, kp. A comparison of the results of these calculations 
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Figure 3.10.- Step function model of the mass absorption coefficient of 

high temperature air. 
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with gray calculations using the Planck mean mass absorption 
coefficient is presented in figures 3.11a to 3 -He. 

The maximum monochromatic Bouguer number for the nongray shock 
layers is 186 times the Planck mean Bouguer number. When the Planck 
mean Bouguer number kp is less than about 0.001 (this case is not 
shewn) the shock layer is optically thin at all wavelengths and no 
perceptable difference between the nongray and the gray calculations 
for the enthalpy distribution can be found. When k^ = 0.01 the 
monochromatic Bouguer numbers for several of the steps are order of 
magnitude unity and absorption becomes important in the nongray 
model whereas absorption is still negligible in the gray model. As 
a consequence, of the above the enthalpy distribution for the non- 
gray model lies above that for the gray model. When kp is 
increased to 0.1, the disparity between the two solutions is 
increased still farther. In this case, absorption is very important 
in those regions of the spectrum for the nongray model in which much 
of the energy is emitted. Absorption is still of minor significance 
in the gray model. Wien kp = 1.0 absorption becomes important in 
the gray model but still not to the extent that it is in the nongray 
case. 

Obviously, and not unexpectedly, a gray model which uses the 
Planck mean mass absorption coefficient will not provide an acceptable 
estimate of the shock layer enthalpy distribution for a nongray gas 
unless that gas is optically thin at all wavelengths in which a 
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significant amount of radiation is transported. Nevertheless, it is 
very interesting, and encouraging to note that enthalpy distributions 
computed for the nongray models do not differ significantly in their 
general shape from those that can be computed for gray models. Thus, 
it appears that there is some wavelength averaged absorption coeffi- 
cient (other than the Planck mean when absorption is important but 
tending toward it in the transparent limit) which will provide a good 
approximation to the enthalpy distribution in a nongray gas. 

The rate of radiant heat transfer to the stagnation point has 
been calculated for nongray shock layers. The results are compared 
in figure 3.12 with the results of gray calculations using the Planck 
mean absorption coefficient. The gray approximation provides a 

considerable overestimate of the radiant heating even for values of 

-3 

the Planck mean Bouguer number as small as 10 . It is apparent from 

this result that the tallest steps play a very important role in the 
transfer of energy by radiation. This is not surprising when one 
considers that nearly ho percent of the energy emitted by a particle 
in the shock layer is transmitted in the wavelength intervals occupied 
by the three tallest steps. 

It can be concluded from the foregoing discussion that the 
effective optical thickness (or Bouguer number) of a nongray shock 
layer is greater than that predicted by a gray analysis using the 
Planck mean absorption coefficient. In order to account for this by 
means of average absorption coefficients, it seems proper to follow 
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the advice of Viskanta (ref. 49) and introduce a "mean emission 
coefficient" and a "mean absorption coefficient." As Viskanta 
pointed out, the divergence of the radiant flux is composed of two 
terms, one of which accounts for emission and the second for absorp- 
tion of radiation in an element of volume of radiating media. In 
particular, for this investigation, the divergence of the radiant 
flux may be written (see eq. (2.67)) 

nCC 

€l Ql]“ ‘ 2e J K-tSn) \(r\)XK+ ekp J €K^( Tj) G^(rj)dA 

(3.34) 

where 


<Vi> ' jf A V"'> Vi’> - ^(V)|)av 

(3-35) 

+ 2r w hfe Vi)) J o y Vi’> Vi') E 2 ( k p Vi')) a i’ 

The first term on the right-hand side of equation (3.34) is the local 
emission term. The integration over wavelength can be performed for 
this term using the definition of the Planck mean absorption coeffi- 
cient (eq. (2.68)) so that equation (3*34) becomes 

£l Ci!) = -2€Kp(T)) B(q) 


+ 6k P Jo ^ (n) °A (Tl)dA 


(3.56) 
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Thus, the mean emission coefficient is identical to the Planck mean 
absorption coefficient. The mean absorption coefficient can be 
defined by the formula 

pOO 

J G^(n)dA 

K a ( n ) = -2 (5. 57 ) 

poo 

1 G (T))dA 
A 

Unfortunately, the spectral characteristics of the quantity G,(n) 
which represents the amount of radiant energy incident per unit mass 
on an element of mass located at q, depends on the spectral charac- 
teristics of the radiating media and the boundary surfaces. Therefore, 
the spectral distribution of G-.(ii) will not be the same as that of 
the Planck function B^( tj) which depends only on the temperature at tj 
and in general k ^ k . The primary difficulty involved in the 
determination of the mean absorption coefficient k is that the 
quantity G^( } is not known rigorously until the whole problem has 
been solved. This difficulty does not arise in the use of the small 
perturbation method of this chapter, of course, because the quantities 
in the equations of various order in e corresponding to G, (n) are 
known rigorously from the solution of the lower-order equations. In 
problems where G^(tj) is not known explicitly, it is hoped that it 
will be possible to obtain a reasonable first approximation. 



CHAPTER IV 


OPTICALLY THIN SHOCK LAYERS 

A. The Transparent Approximation 
Under certain conditions , the Bouguer number , which is indica- 
tive of the optical depth of the shock layer, is very small compared 
to unity. 'Then these conditions are met, absorption is unimportant 
and the absorption integrals which are modified by the Bouguer number 
can be dropped from the expression for the divergence of the radiant 
flux vector (see ea. (2.86)). This leads to considerable simplifi- 
cation because only the local emission rate of radiant energy need 
be considered. All of this radiant energy is assumed to escape the 
shock layer and it matters not, insofar as the gas is concerned, 
what path it takes. Consequently, surface reflectivity will have no 
influence on the enthalpy distribution in the shock layer. Since 
only the total rate of radiant energy emitted locally is of interest 
the details of its spectral distribution can be ignored. 

The results of the simplification is the "transparent" form of 
the divergence of the radiant flux vector 

i[i] - -8«p(n) b(h) (it.!) 

Where l[„] is the divergence of the radiant flux vector, k^(tj) 
is the Planck mean mass absorption coefficient, and B(t]) the Planck 
black-body function. The shock layer is termed transparent because 
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the gas is transparent to its own radiation. Use of the transparent 
approximation reduces the governing equations from integrodifferential 
to purely differential form. Several investigators (see, for example, 
refs, li-7) have taken advantage of this simplicity to obtain approxi- 
mate analytic solutions. 

B. The Optically Thin Approximation 

In this paper, a distinction shall be made between the terms 
"transparent" and "optically thin." A layer of gas will be called 
transparent if none of the radiation emitted by the gas in the layer 
is reabsorbed. An optically thin layer is one in which a small 
amount of absorption does occur and the optical depth of the layer 
is small but not zero. In the literature, "optically thin" is often 
used synonymously with "transparent" as defined above. 

P. D. Thomas (ref. 27) expressed concern about the validity of 
the transparent approximation, particularly in the highly cooled 
region adjacent to the cold wall. The transparent approximation is 
based on the assumption that emission is much greater than absorption 
throughout the shock layer. In regions of small enthalpy, emission 
no longer dominates absorption, and when radiation cooling effects 
are large, these regions may extend over a significant portion of 
the shock layer. Even when radiation cooling effects are small, the 
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value of enthalpy adjacent to the wall tends to vanish and absorption 
must become important compared to the local rate of emission. Of 
course, for this case, the region of nonvalidity is very small and has 
no appreciable effect on overall properties such as the radiant energy 
flux to the wall and the shock standoff distance. 

Thomas sought to modify the transparent equations in order to 
take into account this small amount of reabsorption. He did so by 
expanding the Planck function B^(t) which appears in the integrand 
of the divergence of the radiant flux vector in a Taylor series about 
the zero of the argument of the displacement kernal | T ^ " j. 

The expansion is then arbitrarily truncated after the linear term. 
Strictly speaking, this procedure can be used only when the Planck 
function varies slowly within a photon mean free path length. Obviously 
this criterion is not met when the shock layer is optically thin 
(particularly close to the wall, the region of greatest interest, 
where the enthalpy and hence, the Planck function varies rapidly) and 
some doubt must be cast on the -validity of Thomas' analysis. 

It would appear that the effects of small absorption could better 
be discovered through a straightforward expansion of the equations in 
terms of the Bouguer number kp. Such a solution, up to first-order 

* 

An element of gas approaching the wall requires an infinite time 
to reach its destination. Because of this and the fact that the rate 
of energy lost by radiation is proportional to a positive power of the 
enthalpy, the enthalpy of a transparent gas must approach zero as the 
particle approaches the wall. 
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in kp, is presented here. In order to simplify the analysis the 
exponential integral functions E^(x) and E (x) which appear in 
the expression for the radiant flux are replaced by the exponential 
functions e and (1/2 )e , respectively. The particular form of 
the exponential functions was chosen so that the area under the 
curve of E 2 (x) and the approximating exponential are equal for 
the interval (0, x A ), for x A « 1, and so that the value for the 
radiant flux reduces to the proper value in the transparent limit. 

This substitute kernal approximation has been used with considerable 
success in a variety of problems of radiant transfer (see, for 
example, refs. 11, 30, 31, and 50). 

Use of the substitute kernal approximation reduces the expression 
for the radiant flux to the form 


q R (n) 


6 J c 


r A<*>] 


‘ 2k p t A ( t >W ^ 


d| 


- 21 V\(ti) pH. -21: r (|) . . 
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The divergence of the radiant flux vector is 
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Here the monochromatic optical depth is 


T^(r|) = J K A (n)dTi 


(^) 


It was seen in chapter III that an expansion of the governing 

equations in terms of the small parameter e led to a fortuitous 

uncoupling of the energy and momentum equations. Unfortunately, the 

same is not accomplished when the expansion is performed in terms 

of kp. It is frequently pointed out in the literature (for example, 

* 

refs. 4 and 6) that the coupling is quite weak. Advantage can be 
taken of this situation by replacing pv,(t|) = h( rj ) which appears 
in the momentum equation (2.?l) with h, the integrated average of 
h(r]) over the interval ^0. ; that i6, 


h = — [* h( Tj)dT) 

'A J o 


(It. 5) 


(This replacement is tantamount to solving the momentum equation by 
the integral method with f(r]) assumed to be a quadratic function.) 
This approximation has the twin virtues of retaining the coupling, 


The solution to the momentum equation depends on the solution to 

the energy equation through the term ^.(rj) (see eq. (2.71)) where 

2 ^ 
a is order of magnitude X. An analysis of equation (2.7l) indicates 

that the contribution of this term to f(t|) is order h/X. Since X 
is very small (typically 0.06) the effect of the solution of the energy 
equation to f(q) is order 0.25. 



albeit in approximate form, and greatly simplifying the solution to 
the energy equation. 


Ill 



Now the governing system of equations takes the approximate form 


f(n) h’(q) + €l [V] = 0 (h.C) 

2f(q) f"(tj) -[f'(i)f] + a 2 h = 0 (M) 

f(0) = 0 (! .0) 

f (t) ■ 1 < l! - 5 > 

= 1 ( k -H) 


where I h] j is given by equation (k;). When the Bouguer number is 
very small, the absorption integrals in equation (*<o) assume a 
secondary significance throughout the domain of the problem. Neglect- 
ing these absorption integrals reduces the system to purely differ- 
ential form. If, as expected, v/hen the Bouguer number k_p is small, 
the presence of absorption only slightly influences the solution one 
can, to reasonable accuracy, evaluate the absorption integrals using 
the transparent solution for h. The perturbation expansion scheme 
used herein follows the general outline discussed above. Mathematical 


details are presented in appendix C. 
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The zero-order, or transparent, solution is 



2eil A (l - a*)x + a* 

-^ln 

a x 


(»’. 12 ) 


where 


f (x) = (l - a*)x + a x 



- s y h , 


o ''A 




( lj -lfc) 


( J +.15) 



1 + V2X(1 - X) 

1 + fih x(l - X) 


('- 16 ) 
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h 

o 


/ h (x) dx 

Jo ° 


( * 17 ) 


It was shown in chapter II that the Planck mean mass absorption 
coefficient normalized by its value immediately behind the shock can 
be adequately represented by 


h 2 , h > h* 
h < h* 


r 


Kp(h) =/ 


(U.13) 
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'..’here \f~ h* is the value of the enthalpy (depending on the pressure, 
2 oo 

of course) at which the value of the exponent of h changes. The 
constant is obtained by equating the two expressions for 

Kp(h) at h = h* 


vrith the result 


C ± = (h*) 7 2‘ 7 l 


(’- 19 ) 


It was also shown in chapter II that the nondiraensional Planck 
function B(h) is approximately given by the expression 

B(h) = h 2 ' 2 (b.2Q) 

1-Jhen the correlation formulas (U.l8) and (it. 20) are introduced 
into equation (it. 12) the integration on the left-hand side can be 

carried out, and the solution for h (x) given by the explicit 

o 

formula 


h o (x) = ^ 


= + 2 4 2 + Zn h - a*)x + a* V 1 - 


(i: .21a) 


for h (x) > h*, and 
o' - 


( -(7-+1.2) 

7 , + 1.2 \ 

* -(7 +1.2) 

h o (x) = /(h*) V 2 > 

- g-K 

(h*) ^ 1 1 - 1 


2e(7 1 + 1 * 2 ) c 1 n\, ^ (! " a*)x + a*^~ 7 1 +1 * : 

a x j 


for h (x) < h . 
o v 


.21b) 





The first-order solutions which include the effects of 
absorption, surface reflectivity, and nongray radiation are 
presented below 



Jj(b.22) 
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C. The P-L-K Solution 


Careful inspection of the last term on the right-hand side of 

equation (4.22) reveals that the first-order tern (x) displays 

a singular behavior near tj = 0 (where h Q approaches zero). By 

way of illustration consider the case of a gray gas vith *. p (h) = h 7 
2 2 

and B(h) = h' 1 '* . In this case the term in question is proportional 
to the quantity 


7 + 2.2 


'h h 
o 


dh 

7 + 4 . 4 


v4»r 

h' “■ 
o 


7 + 3-4 


-(7+5-4) 


n 


- 1 


(4.26) 


Wear the wall, h Q approaches zero and equation (4.26) approaches 


h- 1 * 2 

o 

7 + 3-4 


which increases without limit. This , seemingly , nnoMali ;iu» i v-.itu v hu; 

can be explained as follows. The first-order solution represents a 

gas which absorbs radiation at a rate determined by the al-sorotion 

coefficient for a transparent gas* while it emits wnerg} at t rate.- 

proportional to the derivative with respect to h of the emission 

o 

rate for a transparent gas. Both the absorption and emission rates 

The magnitude of the incident radiation is independent me 
amount absorbed. 
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*9 

tend to zero as an element of gas approaches the wall. However, the 
emission rate tends to zero much more rapidly than the absorption 
rate. The difference in the limiting behavior of these rates coupled 
with the infinite residence time for an element of gas in the 
stagnation region allows the gas element to absorb an infinite amo unt 
of energy and so the enthalpy of the gas adjacent to the wall becomes 
infinite. 

The difficulty which has arisen as the result of the singularity 
can be avoided through the use of the P-L-K perturbation of coordi- 
nates procedure which transforms the coordinate in such a way that 
the singularity is removed from the boundary (at x = 0) to a point 
outside the domain of the problem (a slightly negative v al ue of x). 
Mathematical details of the application of this method are described 
in appendix C. The P-L-K solutions to first order in kp are 

x = y + kp x* (y) (k.g ,) 

h(x;kp) = h*(y) ( ' . ... ) 

f(x;kp) = f*(y) + kpf*(y) (' .;.y) 

where the starred coefficients in the P-L-K expansions are related 
to the unstarred coefficients in the regular perturbation exp&iiuions 
(see ref. 43) 


f 
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h*(y) » h o (y) 

(Mo) 

<(y) - f Q (y> 

(Mi) 

'^(y) = f-j^y) + x*(y)f^(y) 

(4.32) 

X x (y) = -h^yj/h^y) 

()■ 

\ ‘ v S*/ / 


D. Results and Discussion 


In order to obtain some indication of the accuracy of the opti- 
cally thin shock layer approximation, the results computed for a 
typical case are compared in figure 4.1 with the results computed by- 
means of the small perturbation method of chapter III and the results 
of a numerical calculation performed by Howe and Viegas (ref. 9). The 
agreement among the three solutions is excellent. However, a word of 
caution should be interjected here in order to avoid the implication 
that the numerical results of Howe and Viegas represent the "exact" 
solution to the inviscid, plane-parallel geometry, stagnation f.lui.- 
model. The results of Howe and Viegas include viscosity, hoax con- 
ductivity, and body curvature. The effects of curvature are expected 
to be quite small. The flight conditions (W = 9-75 km/sea. 


ps = 10 atm) were chosen to insure that the boundary layer was very 
thin so that "displacement" effects on the inviscid region were 
minimized. Finally, the thermodynamic and optical properties used 
by Howe and Viegas were obtained from their own correlations while 
the optically thin and small perturbation methods were computed using 
the correlations presented herein. Thus, the comparisons between the 


results from the methods of this paper and the numerical results of 
Howe and Viegas are as much, or more, checks on the validity of using 




z/A 


Figure 4.1.- Comparison of the optically thin and small perturbation 

solutions with numerical results. 
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the inviscid approach and checks on the similarity of two different 
sets of correlations as they are checks on the accuracy of the ana- 
lytical methods of this paper* It is not inconceivable that erors 
due to the various factors mentioned tend to cancel in uhis example. 
Nevertheless, the individual errors due to the omission ox viscosity, 
heat conductivity, and curvature and due to the difference in correla- 
tion functions are expected to be quite small so that the excellent 
agreement can still be interpreted as an indication of the accuracy 
of the methods of this and the preceding chapter. 

The approximate solution derived in the preceding sections of 
this chapter was used to study the effects of the radiation cooling 
parameter, e, the Bouguer number k , the surface reflectivity r , 
and the enthalpy dependence of the absorption coefficient on the shock 
layer enthalpy distribution, the rate of radiant heat transfer to the 
stagnation point, and the shock standoff distance* As in the previous 
chapter, the density ratio X across the near normal portion of the 
shock was fixed at a value of 0.06. In addition, all the results arc 
limited to the case of a gray gas. 

The effect of absorption on the enthalpy distribution is 
indicated by the curves of figure 4.2. The solid curves represent 
the enthalpy distributions in transparent shock layer for e u O.Oi, 
1.0, and 100. The dashed curves represent the enthalpy distributions 
in optically thin shock layers for the same values of the radiation 
cooling parameter e. Values of the optical thicknesses are shown on 
the figure. These results show the expected trend with the enthalpy 
level falling as the radiation cooling parameter e increases . 
Absorption tends to increase the enthalpy particularly in the cooler 
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regions of the flow. Absorption also affects the location of the shock 
reducing the value of t) (the location of the shock in terras of the 
Dorodnitsyn coordinate) because of the decreased density level. 

Although the value of ti decreases, the shock standoff distance A 
increases with increasing optical depth. 

The effect of the enthalpy dependence of the absorption coefficient 
on the enthalpy distribution in transparent shock layers is shown in 
figure In part (a) the absorption coeffici _ent vac* given by the 

relation ftp = h , where 7 takes on the values 5> ^ > and 5* The 
value of 7 determines how the rate of energy emission varies with 
enthalpy across the shock layer. The rate of energy loss by radiation 
will decrease more rapidly as the enthalpy falls if y is large than 
if it is small. Consequently, the enthalpy distribution for a large 
value of 7 lies above that for a smaller value „ This, of course, 
is the same trend exhibited by the small perturbation solutions of 
the previous chapter. In part (b) the absorption coefficient is 
given by the relation Kp = Cl/ where C = (h’ ) l and 
7 = 7p = ^ for h < h , and C » 1 and y = 7,. *. -1 for u h 
This model should be used when the shock layer t rioperattivc-s erv.- in 
excess of about 20,000° K since at moderate altitudes i/ (the 
enthalpy at which the exponent 7 changes value) corresponds to 
temperatures of approximately this value. The effects of /u vying 

These values of 7 are typical for air at temptratu, iu&s 

than about 20,000° K (see chapter II, section £ ). 
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h are shown in figure 4.3b. A decrease in h produces a decrease 
in the enthalpy level because y takes on the smaller value (-1) 
throughout a greater portion of the shock layer. 

The effect of surface reflectivity r^ on the enthalpy dis- 
tribution is shown in figure 4.4. Of course, this effect vanishes 
in a transparent layer. With a small amount of absorption an increase 
in reflectivity brings about an increase in enthalpy level with the 
greatest increases occurring adjacent to the wall. These results 
corroborate the findings of chapter III. 

The variation with the radiation cooling parameter e of the 
rate of radiant heat transfer to the wall for various values of the 
Bouguer number is presented in figure 4.5. Also shown on this figure 
are two limit curves. One of these curves is labeled the "no decay 
limit" and was computed by assuming that the shock layer was 
isenthalpic and transparent. The second limit curve is labeled the 
"available energy limit" because it represents an upperbound to the 
radiant flux on the basis of energy balance. The amount of energy 
entering the shock layer per unit time per unit area of the shock 
surface has been normalized to unity. If all of this energy is 
radiated out of a transparent shock layer only one -half will be 
incident on the wall. 

The curve labeled kp = 0 shows the effect of "decay" in 
reducing the rate of radiant heat transfer to the wall. The remaining 
curves indicate the important effect of absorption (as characterized 





Figure 4.5*- Effect of the Bouguer number on the rate of radiation heat 
transfer to the stagnation point. 7 = 4.0. 
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here by the Bouguer number kp) in reducing the rate of radiant heat 
transfer to the wall. Although values of k_ presented in figure U.5 
are as large as 3) the corresponding shock layers are all optic all y 
thin ^kpT^ « lj. 

The effect of the enthalpy dependence of the absorption coeffi- 
cient on the rate of radiant heat transfer to the wall in a transparent 
shock layer is shown in figure k,6. It is apparent that an increase 
in the exponent 7 (which appears in the correlati on fox inula Up = h^ ) 
magnifies the effect of decay on the rate of radiant heat transfer. 

The effect of the radiation cooling parameter € on the shock 
standoff distance for various values of the Bouguer number kp and 
7 is shown in figures h.j(a) and k. 7 (b). As expected, an increase 
in e reduces the value of A (the ratio of shock standoff distances 
with and without radiation) for given kp and 7 because the 
cooling by radiation tends to increase the density level in the shock 
layer. Increases in kp and 7 for fixed e inhibits the effect 
of decay on A whereas these increases magnified the effects of 
decay on the rate of radiant heat transfer. 

The variation of shock layer optical thickness k^r with the 
radiation cooling pax*ameter e and the Bouguer number kp is 
shown in figure k.8. When the absorption coefficient varies as a 
positive power of the enthalpy, the shock layer optical, thickness 
may be very much less than one even if the Bouguer number is order 
of magnitude one or greater provided that e is sufficiently large. 
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The criterion for the validity of the analysis presented in this 
chapter is that the optical depth of the shock layer be much smaller 
than one for those wavelength regions in which a significant amount of 
radiant energy is transported. It has been shown herein, for the 
case of a gray gas for which the absorption coefficient is proportional 
to a positive power of the enthalpy, that this condition is always 
less restrictive than the condition that the Bouguer number, kp, is 
very much less than one. However, for the more realistic case of a 
nongray gas the criterion stated above is generally more restrictive 
than the condition kp « 1. In mathematical terms the criterion 
implies the inequality 


; ( 1+ r vK f 0 K X ( ’> ) |X V 


B a (!) 


d£ > dX « 1 (k.3k) 


The quantity on the left-hand side of the inequality is the first- 
order term in the expansion of ijjfJ, the divergence of the radiant, 
flux vector, in terms of the Bouguer number, kp. When both k.^ and 
B-^ are proportional to a positive power of the enthalpy an upper- 
bound to the aforementioned quantity can be obtained by replacing 
and B^(£) by their values at tj = r^, immediately 
behind the shock. The result is 


4 * vK \ £ i{\) \{\)^ 


(^• 33 ) 



133 


If the same substitution is used for a gray gas the result is simply 


'( 1 + r v) k p 


\ 


(4.36) 


because 




and 


are k°th id®htically equal to 1. 


When the nonecrav sten function model fv>v +>ip ^}>QnT r n , f i An aa ci-PP*? ^ 5 + 




of air, which was used in chapter III (see fig. 3.10), is used to 
evaluated the quantity (4.35) the result is about 60 times greater 
than the corresponding gray quantity (4.36). Thus, the criterion for 
the validity of the optically thin analysis, in this nongray example, 
is 


60 kp « 1 

for small values of the radiation cooling parameter e . For larger 
values the criterion could probably be relaxed somewhat ^for example 
60 kp « 1 j . As a result of this criterion the px-actical 
applicability of the optically thin analysis ( and consequently of 
all transparent analyses) is seriously restricted. 



CHAPTER V 


THE OPTICALLY THICK SHOCK LAYER 

A. The Optically Thick Approximation 
A qualitative description of the optically thick shock layer has 
been given by Goulard (ref. 5)* He pointed out that this layer is 
characterized by an isothermal region between two thin boundary layers 
adjacent to the shock and the wall. The boundary layer immediately 
behind the shock is a result of the cooling of the hot gas by radia- 
tion through the transparent shock. Because radiation travels only a 
short distance before being absorbed in an optically thick layer, this 
energy loss is restricted to a narrow region which extends approxi- 
mately a photon mean free path. Once this initial adjustment in 
energy has occurred the gas particle is carried into the interior of 
the shock layer by the flow where convection is the dominant mode of 
energy transport. In this region, the enthalpy of the gas is es- 
sentially constant. As the particle nears the cold wall, moving ever 
more slowly as it does so, convection becomes of decreasing importance 
and energy transfer by radiation begins to assume the major role. 
Finally, in the immediate vicinity of the wall all of the energy 
transport proceeds by means of radiation. When the emissive power 
in the interior (or isothermal portion) of the shock layer is larg:? 
the take-over by radiation occurs at greater distances from the wall 
than if the emissive power is small. Thus, the thickness of the wall. 


13k 
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boundary layer depends not only on the optical thickness of the shock 
layer but on the emissive power of the gas as well. 

While the shock layer is optically thick, the boundary layer 
behind the shock is not and so the Rosseland or diffusion approxi- 
mation so commonly used in the study of optically thick gases cannot 
be applied in this region. The Rosseland approximation is valid 
only in regions of an opaque gas which are at great optical distances 
from all radiation boundaries (a perfectly reflecting barrier is not 
a radiation boundary) and in which the thermodynamic and optical 
properties do not vary greatly within a photon mean free path. 

Neither of these conditions are met in the shock boundary layer. 

The conditions of validity for the Rosseland approximation 
might hold throughout much of the wall boundary layer if the emissive 
power of the gas is stiff iciently large. However, the approximation 
must break down optically close to the wall. The use of a tempera- 
ture jump boundary condition as suggested by several investigators 
(refs. 25, 51> and 52) has proven successful in problems of radiant 
and combined radiant and conductive energy transport. Whether or 
not this concept can be applied with equal success to problems of 
combined radiant and convective energy transport has not, as yet, 
been demonstrated. In a region optically close to a radiation 
boundary the temperature predicted through the use of the Rosseland 
approximation and a slip boundary condition represents not the 
temperature of the molecules of the gas, but a sort of average 
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photon temperature. The convective heat flux depends on the 
molecular temperature. Thus, it is not clear that the slip 
boundary condition can be used in a problem of combined radiant 
and convective energy transport. There is a basis for optimism 
when considering the problem of this chapter, however, in that 
the convective flux may be negligible compared to the radiant 
flux optically close to the wall. 

In order to arrive at a solution to the problem of the 
optically thick shock layer, the substitute kernal approximation, 
introduced in the previous chapter, will be used. It will be 
shown that in the interior of the shock layer and close to the 
wall, but not in the shock boundary layer, this method is equiva- 
lent to using the Rosseland approximation with slip boundary 
conditions. The use of this approximation will restrict the 
analysis to gray gases.* 

B. The Substitute Kernal Approximation 
In this and the subsequent chapter, it will be convenient to 
rewrite the energy equation (equation 2.70 of chapter II) with the 
optical path length t as the independent variable, that is 

f( T ) h'(T) + ei[r] = 0 (3.1) 


*This restriction is not a condition for application of these 
approximations, but has been invoked here to avoid the considerable 
additional complication that relaxation of this restriction would 
incur. 


137 


Here, f is the nondimensional stream function, h the nondimensional 
enthalpy, and e the radiation cooling parameter. The divergence of 
the radiant flux vector, I[r] is given by the expression 

t 

I[ T ] = I[n] = - 2B(r) e 2 (0) + ^ s(t) e^t - t |)dt 

+ 2kp r w E 2 (kpT) A B (t) E 2 (k p tj dt (5.2) 

where t is the nondimensional absorption coefficient, B the 

nondimensional Planck black-body function, kp the Bouguer number, 

r the reflectivity of the wall (at r = 0) , r. the value of the 
w A 

optical path length at the shock, and E^ and E 0 the exponential 
integral functions of first- and second-order, respectively. 

In order to simplify the analysis the substitute kernal 
approximation will be used. For the optically thick shock layer, 
the appropriate substitution for E 2 (x) is found to be (3/h) e v ' // . 

Tliis substitution satisfies the conditions that the areas under the 
two functions over the domain 0 £ x S °° are equal and that the 
expression for the radiant flux approach the Rosseland expression 
as x increases without limit. 

If the expression for the radiation flux is differentiated 
twice with respect to T the integral terms can be eliminated 
with the result 
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I"H - lW = “ | B "( T ) (5-3) 


The energy equation (5.1) can then be used to eliminate I[t]. 

(f(T) h' ( t) 3 - | 4 f ( T ) h '( T ) “ f €B ‘'( T ) (5.4) 


This alternate form of the energy equation is a third-order 
nonlinear ordinary differential equation the solution of which must 
satisfy the condition h(l) = 1. Two additional constants of 
integration are introduced by the solution of (5.4). These constants 
are determined by satisfying appropriate physical conditions or by 
satisfying the original integrodifferential equation (5.1). 

An expression for the flux of radiant energy which enters the 
wall can be obtained quite simply. The expression for the flux 
incident on the wall is 




e 



dt 


When the integrodifferential form of the energy equation (5*1) 
is evaluated at t = 0, it becomes (since f(0) =0 is a boundary 
condition) 


B(0) - l kp (1 + r v ) 



1 

2 


V 


e 


dt = 0 
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Thus the flux entering the wall can be written in terms of the value 
of the black-body emissive power of the gas adjacent to the wall, 
that is 

v(rrr)^ B(0) (5.5) 

C. Boundary Layer Analysis 

In teiros of the substitute keraal approximation, the complete 
differential system governing the flow in the stagnation region of 
a radiating shock layer is 


[f(T)h'(T)]" - | € B"(r) - |k| f(r)h'(T) = 0 

(3.6) 

2f(n)f"(T)) - ff'(ii)^ 2 + a 2 h(n) s 0 

(5-7) 

zx 

0 

II 

0 

(5.8) 

f M * 1 

(5.9) 

f'fnA - a 

^ A ' \[ 2X(l - X) 

(5.10) 

h W - 1 

(5.11) 




dt 



(5.12) 

as 0 
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When the optical thickness of the shock layer is such that 
2 2 

kp » 1 and kp » € equation (5*6) becomes asymptotic to the 
equation 

f(T)h*( T ) = 0 (5-13) 

Thus, the enthalpy approaches a constant. It can be shown by 
attempting to satisfy equation (5*12) as a condition, that this 
solution is valid only at large optical distances from both the 
shock and the wall (unless it is a perfectly reflecting wall). It 
also becomes clear that the value of this constant, hereafter 
denoted h^, cannot be determined without knowledge of the shock 
boundary layer. 

If the enthalpy throughout most of the shock layer is constant 
with a value h^, the density will be constant also with a value 

pg. In this case, the momentum equation may be easily solved with 
the result 



A first approximation to the shock standoff distance is 

A K/g g ) x h g X 

R S l + >y2(p 6 /p 2 jX(i . X) 1 + Ji h 2 X(1 . X) 


(3.1b) 


(5. lb) 
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In addition to the region of constant enthalpy in the interior 
of the layer , there are tvo the rm al boundary layers ; one immediately 
behind the shock and the other adjacent to the wall. 

The forms that the energy equation assumes in these boundary 
layers can be determined by means of conventional boundary layer 
techniques. In the vicinity of the shock the "stretched" coordinate 


i = 



( 5 . 16 ) 


is introduced. Close to the shock the quantity f(r) is slowly 
varying and may be adequately represented by the first term in the 
Taylor expansion about | =0, that is 


fU) ~ = 1 (5.17) 

Substitution of equation ( 5 «l 6 ) and equation (5.1?) into the energy 
equation (5 . 6 ) gives the shock boundary layer equation 


h"( 6 ) 


+ | ^ B '(i) - J h(!) = Const 


(5.18) 


Solution of this equation is complicated by the nonlinear term 
(5/2) —• B 1 ( |) . If e is at least an order of magnitude less than 
lip, this term can be neglected and the solution to equation ( 5 . 18 ) 
easily found. This solution is 


h (0 



+ h„ 


£ 


(5.I9) 



The constant can be determined by writing an energy 

balance across the shock boundary layer. This energy balance is 
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1 = ’Wt/> ) 


-h 

e * d| 


(5.20) 


When € « kp condition (5*20) reduces to 


and it is apparent that h^ approaches one and the boundary layer 
ceases to exist. Thus, there cannot be a shock boundary layer with 
a thickness characterized solely by the optical path length in the 
gas. 

An approximate solution to the boundary layer equation (5.18) 
can be obtained if the nonlinear term (3/2) — - B'(f,) is replaced 

by an appropriate linear term, for example 


— — B'(t) « — — — 

2 kp d^ 


B + Bh(0 




where the constant B is arbitrary and represents a mean variation 
of the black body emissive power B with h over the range of 
values of h encountered in the shock boundary layer. The line- 
arized version of equation (5«l3) has the simple solution 


hU) = (l - hje + h. 


/ _ . _ \ 

D-2i; 



1^3 


where 


®i = i ^ B 




i + 


1/ 1 + *) a 


(5.22) 


where e/kp is very much less than unity this solution reduces to 
equation (5*20) . When ej kp is very much larger than unity, the 
solution takes the form 


h(5) = 



(5.23) 


and, in this limit, the thickness of the shock boundary layer is 
determined by the parameter e ^ instead of simply k^". Thus, 
the shock boundary layer can be very much thinner than a photon 
mean free path if the black body radiative power behind the shock 
is large. This effect was shown by Heaslet and Baldwin (ref. 31 ) 
in their study of radiation resisted shock waves. Simply stated it 
means that a particle starting immediately behind the shock loses 
energy at such a rapid rate by means of radiation that it is 
substantially cooled in the time that it takes to travel only a small 
portion of a photon mean free path. 

A value for the constant B can be obtained from the condition 



2 kp 


(i - - f(i 



(5-24) 


This condition was derived by integration of the nonlinear energy 
equation (5.18) between the limits zero and infinity and substitution 
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into the result of the linearized solution (5 .21). In addition, the 

& 

correlation formula B = h was used. It was shown in chapter II 
that 5 2.2. However, the ensuing analysis will be greatly simpli- 

fied, without any significant loss in accuracy, by setting 5=2. 

A second condition is required to uniquely determine the enthalpy 
distribution in the shock layer. The energy balance relation (5. 20) 
evaluated with the aid of the linearized solution provides this 
condition, which is 



(5.25) 


The quantity can be eliminated between the conditions (5.24) 

and (5*25) resulting in an expression for h , the enthalpy level 

in the interior of the shock layer, as a function of e/k p. The 
result of this calculation is presented in figure 5*1* 

The thickness of the shock boundary layer (in terms of optical 
path length) is characterized by the parameter |u^kpj*"\ A plot 
of ^ as a function of e /^-p * s presented in figure 5*2. 

As has been indicated previously, there is also a thermal 
boundary layer due to radiation adjacent to the wall. If this boundary 
layer is thin, which shall be assumed, herein, the dimensionless 
stream function f (t) may be represented by the first few terms of 
its McLaurin expansion 
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f(-r) *» f(o) + rf »(o) 

Employing the kinematic boundary condition ( 5 . 8 ) and the asymptotic 
solution for f(t|) Equation (5.1k) one finds 


f(-r) 




o Joy ( i 




JoJ^hJ k(o) ^ + yfex(i - x) 


= b T 


(5.26) 


Of the several approximations introduced in the analysis of this 
chapter this is perhaps the poorest because the requirement that the 
wall boundary layer be thin with respect to the optical path length t 
does not necessarily imply that it is thin with respect to either the 
Dorodnitsyn coordinate q or the geometric coordinate z. 

Substituting this expression into the energy equation (5.6), 
introducing the "stretched” coordinate 


, ■« 1/2 

1/2 


fu \ 

Uft I I 


and neglecting terms of order kp yields the boundary layer equation 


b "(0 + h«(o =0 (5.28) 

In general, equation (5.28) is nonlinear and does not possess an 

analytic solution. A simple approximate analytic solution can be 

• # 

obtained by replacing the quantity h'U) with hB ' ( C ) , where h 
is an as yet undetermined constant. This substitution reduces 
equation (5.28) to the linearized form 



where 


0) 



(5.50) 


xne solution to equation (p.29) is easily found with the result 


B(S) - » w + (b 2 - B u )erf Lfe 8 (5.31) 

The quantities B^, the nondimensional black body emissive power 
of the gas adjacent to the wall, and ^ (because it contains the 
arbitrary constant h) are still unknown. One condition for evaluating 
these quantities can be obtained by integrating the nonlinear boundary 
layer equation (5*28) with respect to £ between the limits zero and 
infinity. In performing this integration, it is convenient to 
eliminate the term B"(£) in equation (5.28) with equation (5.29). 

Then it is found that 



h_ 


B_ 



1 

+ h 

w 


(5-32) 


Here B is the nondimensional black-body emissive power of the gas 

in the interior of the shock layer and h is the nondimensional 

w 

enthalpy of the gas adjacent to the wall. The last equality in 
expression (5-32) holds because it has been assumed that B = h . 
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The second condition is obtained by evaluating the integral 
condition (5*12) at t = 0. Written in terms of the boundary layer 
coordinate £ this condition becomes 



(5-53) 


Substitution of the linearized solution for B(£) into equation (5-33) 
and integration yields 



(5. 3*0 


Equations (5*50) and (5*32) can be used to eliminate from 

equation (5*3*0 yielding a transcendental equation for the value of 

the black body emissive power at the wall 

The variation of B as a function of the radiation cooling 

w 

parameter to Bouguer number ratio, G /^p **or various values of kp 

v 

and the exponent 7 (from the correlation formula k •= h 7 ) is 
shown in figure 5*3* This curve has particular significance because 
the ratio of radiant heat transfer to the cold wall is directly 
related to B through expression (5*5)* The variation of the 
quantity ^ ( e< l* (5 *30) with these same parameters is presented 
in figure 5.4. 
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D. The Rosseland Approximation 

The Rosseland or diffusion approximation is frequently used in 
the study of problems in which the medium is optically thick. As 
was pointed out earlier this approximation is not valid in regions 
optically close to a radiation boundary or in regions in which the 
optical and thermodynamic properties vary significantly within an 
optical path length. Some investigators have attempted to correct 
the former deficiency through the use of temperature jump boundary 
conditions and have achieved considerable success in problems of 
pure radiant or combined radiant and conductive energy transport. 

In this section, an attempt will be made to use the Rossel a nd 
approximation and temperature jump boundary conditions to anal yze 
the optically thick shock layer. It is hoped that this exercise 
will provide some insight into the attributes and limitations of 
this approximation in problems of combined radiant and convective 
energy transport. 

With the Rosseland approximation for the divergence of the 
radiant flux, the energy equation becomes 

3 *P 

B"(t) + | “ f(T)h'(r) = 0 (5.36) 

This equation is the same as equation (5 .6) except for the omission 
of the third-order differential term (f(T)h'(T)] . 

In the interior of the optically thick shock layer, equation 
(5*3 6) reduces to 
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f(T)h'(r) a 0 

provided kp » — . This result is identical to the result obtained 
by means of the substitute kernal approximation. This agreement is 
not surprising because the diffusion approximation is known to be 
valid in this region. Of course, the value of the constant enthalpy 
in the interior of the shock layer cannot be determined until some- 
thing is known about the shock boundary layer. 

If the usual type of boundary layer analysis is applied to the 
energy equation in the Rosseland approximation (5 .36) for the 
neighborhood of the wall the result is 

B"U) + | (x ) b 5 h '^) - 0 

This equation is identical to the wall boundary layer equation in 
the substitute kernal analysis. Two boundary conditions are 
required to completely determine the solution to this equation. 

One of these conditions is 

lim B(C) = Bp 
C-»" 2 

Here B^ is the black-body emissive power of the gas in the interior 
of the shock layer. The second is the jump boundary condition, 
written in terms of the black-body emissive power B rather than 
the temperature. 



l6o 


B w = CB'(O) 




(5.37) 


The second equality follows from the expression for the radiant fl ux 
in the Rosseland approximation. The constant C is usually 
evaluated by requiring the flux to be correct in the black-body limit 
(see, for example, ref. 29 ). However, it is noted that condition 
(5-37) is identical to the condition used in the substitute kernal 
approximation (i.e. (5*33)) if C is chosen to be 


1 + r 


V 


M 1 - r v 


(5.38) 


Thus, the results obtained in the wall boundary layer by the two 
methods are identical if C is chosen to satisfy ( 5 . 38 ). It has 
been shown that the two methods also give identical results when 
applied to the problem of combined radiative and conductive energy 
transport between two infinite parallel plates separated by a radi- 
ating and conducting gas (ref. 29 ). 

If the usual boundary layer analysis is used to obtain the 
boundary layer form of the energy equation in the Rosseland 
approximation for the neighborhood of the shock the result is 


B '(S) - §(^T^ h (S) =• Const ( 5 . 39 ) 

This equation is not identical with the shock boundary layer equation 
in the substitute kernal approximation because of the omission of the 
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third-order differential term. Inspection of equation (5-39) 
indicates that any solution other than the trivial solution 

~ 1*2 will not tend to a constant as 5 becomes very 

large. Thus equation (5-39) cannot be forced to simultaneously 

satisfy the conditions h(o) = 1 and lim h(|) = h 

£->«> 2 

Apparently then the jump boundary condition at the shock must be 
h(o) = hg, but this result leads nowhere as there is insufficient 
information to accurately determine hg. 

It must be concluded then that the Rosseland approximation 
with slip boundary conditions is not sufficient by itself to be 
used in the analysis of the complete shock layer. It can be used 
in the combined radiation and conduction problem because the two 
separate energy fluxes are represented by similar mathematical 
models and may be treated as an equivalent radiation alone or 
conduction alone problem. Even in this case, it is not possible 
to obtain a temperature distribution nor to separately determine 
the radiant and conductive fluxes optically close to a boundary. 
In the combined radiation and convection problem, this inability 
to determine a temperature distribution or to separate the radiant 
and convective fluxes optically close to a boundary (such as a 
transparent shock) where convection is important presents a 
serious obstacle to solution because the convective flux depends 
on the unknown temperature (or enthalpy) distribution. 
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E. Results and Discussion 

The analysis presented in this chapter applies only when the two 
enthalpy boundary layers are thin compared to the shock standoff 
distance with these distances expressed in terms of the Dorodnitsyn 
variable r\. When ej kp, the ratio of the radiation cooling 
parameter and the Bouguer number, is much less than one, the 
thickness of the shock boundary layer is characterized by the inverse 
of the Bouguer number, k^ 1 , while when is large the shock 

boundary layer thickness is characterized by the inverse of the 
radiation cooling parameter, e’ 1 . The thickness of the wall boundary 
layer is characterized by the parameter . Thus, the most 

restrictive conditions on the applicability of the optically thick 
analysis are that kp »1 for e small and kp » e 1 / 2 for e 
large . 

Several enthalpy distributions were calculated with the formulas 
developed in the preceding section. The results are presented in 
figures 5.5a and 5 -5b. The previous discussion of the effects of the 
parameters on the thicknesses of the boundary layers is substantiated 
by these results. The effect of the Bouguer number, kp, and the 
radiation cooling parameter to Bouguer number ratio, € on the 
shock layer optical thickness kpT^ is also shown. The effect of 
€ / k P de P ends on the enthalpy variation of the absorption coefficient. 
In the cases shown the absorption coefficient is proportional to the 
fourth power of the enthalpy and an increase in € /kp brings about a 
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| Figure 5.5.- Effect of the parameters e/kp and kp on the enthalpy 

L distribution in an optically thick shock layer. 
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reduction in the shock layer optical depth. The value of the 

enthalpy of the gas adjacent to the wall (which is related to the 

ratio of radiant heat transfer to the wall through equation (5.5) 

2 

and the correlation formula B = h ) decreases with increasing 
e^kp and/or kp. 

The effect of 7 (where 7 is the exponent in the correlation 

y 

formula k * h ) and the surface reflectivity r on the character 

w 

of the wall boundary layer has not been shown but can be readily 

deduced from the curves of figures 5.5 and 5.k. Increasing 7 tends 

to reduce the optical thickness of the wall boundary layer and 

increase the value of the enthalpy of the gas adjacent to the wall. 

It can be shown that the wall boundary layer thickness expressed in 

terms of the Dorodnitsyn coordinate tj is only slightly effected 

by a change in 7. Increasing the surface reflectivity r tends 

w 

to increase the optical thickness of the wall boundary layer and 
increase the value of the enthalpy of the gas adjacent to the wall. 
When expressed in terms of the Dorodnitsyn coordinate t} the 
boundary layer thickness decreases with increasing r . These 
results are consistent with the results of Chapter 5. 

The manner in which the rate of radiant heat transfer to the 
wall, q^, depends on the radiation cooling parameter to Bouguer 
number ratio, e/kp, the Bouguer number, kp, the variation with 
enthalpy of the absorption coefficient ( through the exponent 7 of 
the correlation formula k = h?) , and the surface reflectivity, r^, 
is indicated in figures 5»3& to 5 -3b. For fixed values of kp, 7, 
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and i* w the rate of radiant heat transfer to the wall, q^ , 

t> 

would 

become asymptotic to the available energy limit of l/2 as 

ej kp increased without limit. As the Bouguer number, kp, 

increases (hence increasing the shock layer optical thickness), 

all other parameters remaining fixed, the value of decreases 

and becomes asymptotic to zero. This is the same trend exhibited 

in the problem of radiant energy transfer between two plane 

parallel walls separated by an absorbing and emitting, but 

motionless and nonheat conducting gas (see, for example, ref. ll) . 

Increasing 7 while holding the other parameters fixed results 

R 

in an increase in q^. This trend is the reverse of that for a 

transparent shock layer (see fig. 4.6). The results of Chapter 3 

(see fig. 3*6) show that this reversal occurs at intermediate 

values of the Bouguer number kp. Finally, it is apparent from 

figure 5.3 that an increase in surface reflectivity r, 

w 

values of the other parameters results in a decrease in 

p 

change in with r^ satisfies the inequality 



which agrees with the physical argument presented in section 4 of 
chapter III. 


for fixed 
R 


increases with increasing e/k p. It appears that 



CHAPTER VI 


THE RADIATION DEPLETED SHOCK LAYER 


A. The Strong Radiation Approximation 
When the radiation cooling parameter e is very much greater 


+>. 


i on both one snd y the Bougizsr number s ed ^ 9, p8.2rti.d0 


leaving the shock with an initial enthalpy of i will emit 
radiation at such a rapid rate that it will lose a significant 
amount of its energy before traveling the distance of a photon mean 
free path. Because this energy is emitted in a region of snail 
optical thickness adjacent to the transparent shock most of it will 
escape from the shock layer, and the enthalpy level within the shock 
layer will be quite small in comparison to the value at the shock. 

In fact, as will be shown subsequently, the zero -order solution for 
the enthalpy in the interior of the shock layer is identically zero. 

It is for this reason that the term "radiation depleted shock layer" 
has been coined. Of course, the narrow region adjacent to the shock 
in which the large change in enthalpy occurs can be described as a 
boundary layer and boundary layer techniques can be applied to obtain 
solutions in it. 

The conditions which must hold in order that there be a radiation 

2 

depleted shock layer, that is € very large and kp not too large, 
occur at high altitudes for rather large objects (shock radius greater 
than 1 meter) entering at extremely high speeds (entry speeds close 
to 70 Ion sec ) . It is not at all clear that the requirement for 
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chemical equilibrium can be satisfied under these condition, particu- 
larly in view of the existence of a shock boundary layer in which 
large changes occur over a short distance, and hence, a short time. 
Dispite this objection, the solutions for the radiation depleted shock 
layer represent an interesting limiting case which should lead to an 
increased understanding of the radiating shock layer and provide a 
firm base for extension into areas of more practical concern. 

B. Analysis 

Once again, as was the case for the optically thin and optically 
thick shock layers, analysis can be facilitated through the use of 
the substitute kernal approximation. In this case, the energy equation, 
written in terms of the optical path length, is 

[f(T)h«(r3 " - | € B"(t) - | k^ f ( t ) h * ( t) = 0 ( 6 . 1 ) 

Here f is the nondimens ional stream function, h the nondimensional 
enthalpy, B the nondimensional black -body emissive power, e the 
radiation cooling parameter, kp the Bouguer number, and t the 
optical path length. It should be remembered that use of this equation 
restricts the analysis to gray gases only. The boundary conditions on 
the enthalpy are, as before 

h^r A j = 1 (6.2) 


and the integral condition 
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The particular form of the Substitute prrrnl 

E.(t) - 


V, £5k>*^i 4’V\<^4- -? ^ 

'-'J v ^'-u- av>a vaicx ^ 


4 


- — k_T 

2 h 3 

e , vas chosen for simplicity. Somewhat greater 


accuracy night be achieved by letting the constants depend on the 
optical depth kpT^. However, it vas not felt that this procedure 
would lead to a better understanding of the radiation depleted shock 
layer. 

The momentum equation, in terms of the Dorodnitsyn coordinate, 
and the boundary conditions on the nondimens ional stream function 

are 

2f(q)f”(q) - [f'(q)j + a 2 h = 0 (6.4) 


f(0) = 0 

(6.5) 

f M - 1 

(6.6) 


f 




a 

y2x(i - x) 


(6.7) 
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When the radiation cooling parameter e is very much greater 
2 

than one and k the energy equation (6.1 ) admits the asymptotic 
solution 


E(t) = C x + C 2 t (6.6) 

Substitution into the asymptotic form of the integral condition 
(6o) reveals that each of the two constants must be identically zero. 
Thus , in the interior of the shock layer B(t) and h(r) are aero. 

In this case, the density is infinite and the momentum equation can 
ce readily solved for f(q) with the result 


f( "> ’ 

and, of course , the shock standoff distance tends to zero. 

In order to investigate the shock boundary layer, it is 
convenient to introduce the "stretched" coordinates 


( 6 . 9 ) 


and 



( 6 . 10 ) 


( 6 . 11 ) 


into the energy and momentum equations, respectively. Performing 

the usual manipulations (details are presented in appendix D) shows 

that the boundary layer is characterized by the parameter and 

it would seem proper to expand both the boundary layer and asymptotic 

solutions in power series of this sm all parameter. However, f (q) 

a 
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(where the subscript a indicates the asymptotic solution valid far 

-1 


from the shock) is not analytic in e “ near 

:-l/2 


e "*" = 0, but is 


Consequently, the solutions must be expanded 

-V 2 . 


analytic in e 
as power series in e 

The lowest -order form of the energy equation in the boundary 
layer is 

( dh 


V 


! ’b + 2 B b <*> = C 1 + C 2 S 


(6.12) 


The subscript b Q has been used to denote the zero-order boundary 

layer solutions. This equation must satisfy the boundary conditions 

that both ( £) and B’ ( £ ) vanishes as £ increases without 
o o 

limit. Thus, the constants and C ? are both identically zero. 

The third condition to be satisfied is 

B b (0) = 1 (6.13) 

o 

The solution to equation (6.12) subject to the boundary conditions is 


£ 



r 1 fm 

\ |sB 


dB 

B 


( 6 . 1 * 1 ) 


Solution of the momentum equation in the boundary layer gives 
the zero -order form of the nondimens ional stream function 


f b (?) = i 


(6.13) 
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These zero-order solutions can be used to generate solutions 
of higher order. Mathematical details are presented in appendix D. 

In general, the analysis follows the procedure outlined by Van Dyke 

(ref. ;>6). 

A complete listing of these solutions up to second-order in the 
_ 1/2 

parameter e ' is presented below. 

Zero-order solutions. - 

B (t) = 0 (6.16) 

a 

o 

v n) ■ (tf (6 - rn 



(6.16) 


f b (C) = 1 (6.19) 



2 V2X(l - X) 

a 


1 + y 2 x(i - x) 


(6.20) 



( 6 . 21 ) 


First-order solutions . - 


B (r) = 0 

a l 


( 6 . 22 ) 
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(i) = a \ JwQt) 



(6.25) 


Bh (I) = 0 (6. 2k) 

f b = 0 ( 6 . 25 ) 

1 




17 ^ 


♦ 


*- 


'/here 




k P 





/n. V 




>■ r 





(6.30) 



+ 


£T. +€^ ' j +T. 

' A o A 1 ^2 



( 6 . 31 ) 


Radiant heat -flux and standoff distance . - The total radiant heat 
R 

flux to the wall and the ratio of the shock standoff distance to 

the shock standoff distance for radiationless flow A are given by 
the following expressions. 


R 


1 - r 


w 


1 + 


i 1 - r v) 


k. 


P A 


( 6 . 52 ) 
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S' 

( K p i(o) ) 

h(0) + h(0) i-S, 


7 

' 1 % 

v (o) 

1 + J "p T i 

1° 

K P ^ - 

2 

L o_ 



( 6 . 33 ) 


C. Results and Discussion 

In the analysis of the preceeding section, it was convenient to 
use the black-body emissive power B rather than the enthalpy h as 
the dependent variable. This choice necessitated the assumption that 
the thermodynamic and optical properties (in particular h and k” 1 ) 
be analytic functions of B in the interval (0,1). Unfortunately, 
this condition does not hold for the correlations of chapter II 
written in terms of B in the limit as B approaches zero. This 
difficulty can be circumvented through the use of analytic approxi- 
mations to the correlating functions. For example, the enthalpy 
might be approximated by the function 


h = (B + B*) 1 / 2 


(6.3*0 


where B* i S a very small positive number. Use of formula (6. 3^-) in 

expression (6.13) results in the following solution for B^ (l) 

o 


B b U) = B* 


1 + Ce 


1 - Ce 




1 


( 6 . 35 ) 
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where 

C = l/l + B* - l/~B* 
l/l + B* + V? 

POI* voliiac r\-P t -4-Vv o 4- ■? «-. 

* v * t>'*' * •»/ ^ y uuu u» J.O 

i » 1/3 yi* 

the value of B fc ( £) is directly proportional to B? It is 
* 0 

then that B should be chosen sufficiently small to insure 

B b (|) is nearly independent of B* for values of B, as 
o o 

(0). 

S 2 

Because of the unlikelihood of establishing local thermodynamic 
and chemical equilibrium in a physical shock layer under those 
conditions for which this model analysis is supposed to apply, it 
would be somewhat superfluous to present the results of detailed 
calculations for the enthalpy distribution and shock standoff 
distance. Suffice it to say that the shock layer is characterized 
by an enthalpy boundary layer immediately behind the shock the 
thickness of which (in terms of the Dorodnitsyn variable q) is 
characterized by the inverse of the radiation cooling parameter e”\ 
It should also be pointed out that the shock boundary layer is always 
very much less than a photon mean free path and hence is always 
optically thin. The enthalpy level in the interior of the shock 
layer is order of magnitude k_,/e. The ratio of the shock standoff 


( 6 . 36 ) 


clear 
that 
small as 



distance to the shock standoff distance for radiationless flow is 


order of magnitude e 


Curves representing the magnitude of the radiant heat flux which 
is absorbed by the wall are presented in figure 6.1. In the 

optically thin limit fkp « 1J the radiant heat flux approaches 

the "available energy limit" of ^1 - r vjj 2 ' As the thickness 

of the wall increases and absorption becomes more important, less of 

the energy emitted from the gas in the shock boundary layer in the 

direction of the wall is able to penetrate the shock layer and reach 

the wall before being absorbed. Part of what is absorbed is then 

reradiated in the forward direction and escapes from the shock layer 

through the transparent shock. Finally, as k^ t. tends to infinity 

o 

none of the energy emitted in the shock boundary layer reaches the 
wall and the radiant flux incident on the wall vanishes. 




CHAPTER VII 


radiating shock layers 


A. Discussion of the Approximate Solutions 
Four different approximate stagnation point solutions for an 
inviscid, radiating shock laver VPrp nhf*.pn r> p>rl i n r%rv 

^ V w ~ ~ — v^vuj.llg, X-AitAjL/ V-rV^J. U 1 

bach one represents a limiting case for some combination of the 
radiation cooling parameter e and the Bouguer number kp. The 
regions of validity of the approximate solutions are depicted in 


figure 7*1. The boundaries as drawn pertain only to a gray gas 
with constant absorption coefficient. It would be necessary to redraw 


the boundaries for each nongray gas and for every change in the 


enthalpy dependence of the absorption coefficient. As was pointed 
out in chapter IV, the criterion for validity of the optically thin 
solution is that the gas layer be optically thin in all wavelength 
intervals in which a significant amount of energy is transported by 
radiation. For a gray gas, this means kpT « 1. Thus, the 
boundary is not specified completely by kp but varies with e as 
well (since kpT^, the shock layer optical thickness depends on e 
as well as kp when the absorption coefficient is a function of the 
enthalpy) . When applied to a nongray gas, the criterion for 
validity of the optically thin solution is always more restrictive 
than the condition that the Planck mean optical depth be small. Thus, 
the boundaries for all nongray gases will be displaced to the left 
of the boundary for the "Planck-equivalent" gray gas. A "Planck- 
equivalent" gray gas is one in which the wavelength independent 
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absorption coefficient is equal to the Planck mean absorption 
coefficient in the nongray gas. 

The location of the boundary for the small perturbation 
approximation depends on the radiation cooling parameter € and 
the enthalpy variation of the absorption coefficient. The value 
of £ for which the solution will yield results of a given accu- 
racy is reduced with an increase in 7 (where 7 is the exponent 

y 

in the correlation formula K p = h ) , because of the reduced 
accuracy of the truncated expansion for tc p (equation (B-43) of 
appendix B) . Since the small perturbation solution was shown to 
be correct to second-order throughout most of the domain of the 
problem the condition for establishing the boundary is e 2 « 1.* 
The location of the boundary does not depend on the wavelength 
dependence of the absorption coefficient. 

The most restrictive condition limiting the validity of the 

optically thick analysis for moderate values of the radiation 

cooling parameter e is the thickness of the enthalpy layer 

adjacent to the shock. This thickness is characterized by the 

inverse of the Bouguer number l^kp. Thus, the criterion for 

validity is kp » 1. For larger values of € the condition 
( I 2V-/2 

[e/kpl « 1 becomes more restrictive and must be used to 


establish the boundary. This latter condition insures that the 
enthalpy boundary layer adjacent to the wall is thin compared to 


This condition holds when 7 = 0. When 7 = 4 the proper 
condition becomes 10 « 1. 
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the shock standoff distance. The analysis presented in chapter V 
is restricted to the case of a gray gas but could be extended 
rather simply to the case of an absorption coefficient with a step 
function dependence on wavelength for which the step heights are 
either a(h) Kp(h) or zero. There is no restriction to the 
number of width of the steps. The only changes that would appear 
in the formulas would be the substitution of a(h) Kp(h) for Kp(h) 
and B(h)/a(h) for B(h) . The boundary to the region of validity 
of the optically thick shock layer analysis would be displaced to 
the left for this particular class of nongray gases. 

The region of validity of the radiation depleted shock layer 
analysis is restricted by the conditions e » 1 and e » k^. 

The first condition insures that the thickness of the enthalpy 
boundary layer adjacent to the shock is small compared to the shock 
standoff distance, while the second condition insures that radiation 
is the preponderant mode of energy transport within the shock layer. 
Like the analysis of chapter V, the radiation depleted shock layer 
analysis is restricted to gray gases but can be extended to the 
nongray model absorption coefficient with multiple steps of uniform 
height. Use of such a nongray model would cause a leftward shift 
in the boundaries to the region of validity. Of course, the regions 
of validity of both the optically thick and radiation depleted shock 
layer analyses must vanish for all other classes of nongray gases. 

In order to relate the radiation shock-layer regimes to the 
problem of entry into the atmosphere of the earth, several 
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trajectories are indicated on the e - kp map presented in 
figure 7.2. The arrows indicate the direction of increasing time. 
Trajectories 1 and 2 represent iron spheres of radius 1 meter and 
1 centimeter, respectively, entering vertically with an initial 
velocity of 70 km/sec.* Trajectories for all other objects of 
the same size and lesser or equal initial velocities must lie below 
them in the e - kp space. The third trajectory corresponds to 
the entry of a round-trip Martian probe which would encounter some 
of the more severe heating conditions of the currently envisioned 
class of m a nma de objects. It is apparent that the small pertur- 
bation approximation has considerable utility for the analysis of 
radiation effects on the entry of manmade objects. It also appears 
that the optically thin shock layer analysis might enjoy wide 
applicability. However, in the more realistic case of a nongray 
gas the boundary would be shifted to the left one or two orders of 
magnitude in the Bouguer number kp, considerably reducing the 
practical usefulness of the optically thin approximation. The 
optically thick and radiation depleted shock layer analysis would 
seem to be nearly devoid of direct practical usefulness, both 
because of the inacessibility of the proper magnitudes of the 
parameters e and kp to objects of interest and because of the 
restriction of these analyses to the gray case (and the simple 

■ft *** “ ** — - — — — ■ ■ ■ ■ 

No account has been taken of mass loss of these spheres due 
to ablation. 
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nongray model absorption coefficient with multiple steps of 
uniform height) . 


B. A Model Earth Entry Environment 
The four approximate solutions can be used to compute the 
radiant heat transfer to the stagnation point over a wide range 
of the radiation cooling parameter e and the Bouguer number kp. 
The results depend on the particular gas, the surface reflectivity, 
and the size of the object and must be recomputed for every change 
in these variables. Actually, the size of the object is important 
only if the exponent 7 (which appears in the correlation formula 
Kp = h ) varies throughout the e - kp space. In this event, the 
value of h at which a change in 7 occurs depends on the para- 
meter kp =■ P g Kp A^ which is influenced by the body size through 
s 

the radiationless shock standoff distance A, . 

A 

Contours of constant values of T (the temperature immedi- 

s 

ately behind the shock), p k (The Planck mean volume absorption 

s s 

coefficient immediately behind the shock), e /kp, and X (the ratio 
of free- stream density to the density immediately behind the shock) 
on plots of ambient density ratio versus free- stream 

velocity W ro up to 70 km/sec are presented in figures 7* 5a through 
7.5d> for a model earth entry environment. This entry environment 
was obtained by combining the thermodynamic and optical property 
correlations presented in chapter II, section E with the strong 
normal shock relations. The resulting formulas are: 
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T = 1.058 x 10 5 (— ) W 1 * 28 , °K 

S V°SL ) " 


( 7 . 1 ) 


X = ( — ] « 6.95 x 10 

\ o / \ » 
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08 


( 7 . 2 ) 


7 - = 1.852 X 10' 16 T 4 
*P s 



( 7 . 5 ) 


P Kp = 7 . 9 ^ X 10 
s 


-26 f w . 1 ^V-25 m 6.0 - 0.5 Log 10 fx" 1 — 1 __-l 


St, 


p siJ' cm 


for the lover temperatures (less than about 20,000° K) and 


\7 • 4a) 


o s K p _ = 9.35 x 10 2 /V 1 y - 507 T s - °* 39 + °- 21 Log 10 ( X_1 o SL ), cm - 1 


( 7 . 4 b) 


for the higher temperatures. 

The values of X and 7 (the exponent in the correlation 
y 

formula K p = h ) do not vary greatly over a rather extensive range 
of ambient densities and velocities. Consequently, it was decided 
to fix these quantities at the constant values, X = 0.06 and 
7 =4.0, for the discussions which follow. 
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C. Radiant Heat Transfer 

The rate of radiant heat transfer to the stagnation point of a 
R 

blunt object, q, was calculated by the four approximate methods 

for a wide range of the radiation cooling parameter e and the 

Bouguer number lip. The results are presented in figure 7.4 as a 
R 

plot of against kp for various values of the ratio e/k p. 

This ratio, sometimes known as the inverse of the Boltzmann number, 
was used because it is what might be termed an "environmental 
parameter," that is a parameter dependent only on free- stream 
conditions (ambient density and velocity) and not on body geometry. 
The Bouguer number kp, on the other hand, is directly proportional 
to the body nose radius for a given set of free-stream conditions. 
Thus, each curve in figure 7*4 can be thought of as representing the 
effect of body nose radius on radiant heat transfer at a given 
trajectory point. 

For the purpose of calculating the results presented in figure 
7.4, the shock layer gas was assumed to have a gray mass absorption 
coefficient which varies as the fourth power of the enthalpy. The 
surface of the object was considered to be nonreflecting. The dashed 
portions of the curves do not represent computed data, but rather 
represent arbitrary connections across regions in which the various 
approximate solutions are invalid. 

The radiation cooling parameter e is equal to the radiant flux 
leaving each side of a transparent, isenthalpic gas slab in which the 
nondimensional enthalpy tak.es the value one. Hence, this product 
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represents an upperbound to the rate of radiant heat transfer to the 
stagnation point (or wall), q^. When this product is small, the rate 
of energy loss through radiation is small, and the average intensity 
is only slightly perturbed from the isenthalpic value. However, as 
e increases (the Bouguer number kp remaining very much less than 
one) the increased energy tost by radiation is reflected in 
decreased levels of enthalpy and average intensity. Hence, 
becomes a decreasing fraction of €. Finally, as e becomes very 
large, <*p still small) nearly sill of the energy is removed from 

p 

the shock layer by radiation and q^ which represents the rate at 
which radiant energy leaves one side of the transparent layer, 
approaches the physical maximum of l/2. 

As kp increases toward and beyond unity, absorption becomes 
important and this mechanism, which tends to inhibit radiant energy 


R 


transfer, halts the increasing trend of q^ with kp. As kp 


R 


and 


continues to increase, the trend is reversed and q^ decreases 

becomes asymptotic to zero. Consequently, the curves of rate of 

, R 

radiant heat transfer to the stagnation point against Bouguer 

number kp for constant values of the ratio e jk^ have maximums 

the locations and heights of which depend on €^kp. It can be 

inferred from this that for every altitude and velocity in this 

simple model atmosphere, there is a finite value of nose radius for 

which the rate of radiant heat transfer to the stagnation point will 


be a maximum. 
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In order to obtain some understanding oi the effects of 
radiation cooling, gray absorption, and spectral absorption of the 
rate of radiant heat transfer to the stagnation point a series of 
calculations utilizing various approximations were performed. The 
results of these calculations corresponding to a free-streara velocity 
O-f 14 .2 km /sec and an altitude of 32 . 4 km are plotted against body 
nose radius in meters in figure 7. 5- The curve labeled 1 was 

computed by assuming that the shock layer was both isenthalpic and 
nonabsorbing. In this case the rate of radiant heat transfer to 
^he stagnation point is given by the simple expression, = e. 

This approximation was used in the early estimates of radiant heeling 
(refs. 1 and 2). Curve number 2 was computed by assuming that the 
shock layer was isenthalpic and contained a gray, absorbing gas. The 
effect of gray absorption is seen to be small (under the conditions 
of this example) for a nose radius as large as 0.1 m. The third 
curve was obtained using the transparent approximation discussed in 
chapter IV. This assumption of a nonabsorbing but radiation cooled 
shock layer is frequently employed in the literature (see, for 
example, refs. 3-7). For this example, at least, the effect of 
radiation cooling is more important then the effect of gray absorption 
for nose radii of 0.1 m or less. Curve number 4 contains the effects 
of both radiation cooling and gray absorption. These combined effects 
are included in the numerical solutions of Howe and Veigas (ref. 9). 

It can be seen that for small nose radii (less than about 0.1 m) gray 
absorption has little effect. However, gray absorption plays an 
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increasingly important role as the radius increases. The final curve, 
mannered 5, includes the combined effects of radiation cooling and 
nongray absorption. The absorption coefficient used in these 
calculations was the step function model introduced in chapter XIX 
(see fig. 3.10). The curve is limited to small values of nose ]ro.dius 
because of the restricted region of validity of the small perturbation 
method with which this curve was computed. “It is very apparent from 
these results that nongray effects cannot be ignored if one wishes to 
obtain a realistic evaluation of the radiant heating of objects 
during entry at hyperbolic velocities. 

The analysis of this paper has been restricted to a shock layer 
with plane- parallel geometry. The largest effect of this assumption 
is felt in the calculation of the rate of radiant heat transfer. 

Koh (ref. 19) has shown that the plane-parallel geometry assumption 
can lead to an overestimation of q_. by no more than 15-percent 
when the gas is transparent to its own radiation and when the shock 
standoff distance to shock radius ratio is no greater than 0.05 . As 
the Bouguer number kp increases, the size of the error decreases 
and vanishes when the shock layer becomes optically thick. Because 
the effective optical thickness of a nongray shock layer is greater 
than that for a Planck equivalent gray gas, the error due to 
geometry will be smaller for a given Bouguer number in the more real- 
istic nongray case. 
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D. Convective Heat Transfer 

Even though the analysis of this investigation is based on the 
assumption that the gas in the shock layer is inviscid and nonheat 
conducting, it is possible to draw some conclusions regarding the 
coupling between radiant heat transfer and convection heating. The 
convective heating rate (sometimes referred bo as the aerodynamic 
heating rate) is, in the case of a laminar boundary layer, the rate 
at which heat energy is transferred to the body surface by means of 
conduction. 

To first order in the boundary layer naram.-'ter Pe -1 ^* (see 
section D of chapter II) the convective heating rate is proportional 
to the enthalpy difference across the conduction boundary layer. 

If the wall is cold (as has been assumed throughout this investi- 
gation the enthalpy of the wall can be neglected and the convective 
boating rate becomes proportional to the enthalpy at the outer 
edge of the boundary layer. The location of the outer edge depends 
"pon the Peclet number. Since it has been assumed throughout this 
investigation that the viscous boundary layer is thin (in terms of 
both the Dorodnitsyn coordinate and the optical path length) the 
location of the edge of the viscous boundary layer will be arbi- 
tratily specified as Tyry ^ 0.05 for both the small perturbation 
and the transparent solutions. The rapid change ir. enthalpy near 
the wall, particularly for the transparent approximation which 

The inverse square root of the Peclet number. 
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gives a value of zero for the enthalpy of the gas adjacent to the 
w all , necessitates choosing an edge location other than zero. For 
the optically thick and radiation depleted shock layers, it is 
more convenient to specify the edge of the viscous boundary layer 
in terns of the normalized optical path length t. The variation 
of enthalpy near the vail is quite small in the case of the 
radiation depleted shock layer. Consequently, the edge of the 
viscous boundary layer can be considered to be located at r = 0 
for this case. A vail boundary layer due to radiation has been 
shown to exist in the optically thick shock layer. This wall 
boundary layer is always thicker than a photo luGS n free path, and, 
of course, is very much thicker than the optically thin viscous 
boundary layer. Therefore, t = 0 can be considered as the edge of 
the viscous boundary layer for this case also. Values of the 
enthalpy h^ at tht edge of the viscous boundary layer have been 
determined from the four approximate solutions for a wide range of 
the ratio of the radiation cooling parameter to the Bouguer number 
e / k p and the Bouguer number kp. The results are shown in figure 
7 . 6 . The dashed portions of the curves represent arbitrary 
connections across regions of nonvalidity. 

The quantity h g is a rough approximation to the ratio of the 
convective heating rate for a radiating shock layer to that for a 
nonradiating shock layer. When radiant energy transport is important, 
the convective heating is reduced from the radiationless value 
(h g = l). The effect becomes larger as both e /k> increase. It is 
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Interesting to note than the convective heating continues to decrease 
for increasing kp even vhen the shock layer is optically thick and 
the rate of radiant heat transfer is decreasing as a result of ab- 
sorption. Even though the total heating rate (radiant plus convec- 
tive) cannot, be deduced from an inviscid analysis, it is apparent 
that the total heating rate decreases vith increasing shock layer 
optical thickness for all values of kp at least as large as the 

value for maximum rate of radiant heat transfer to the stagnation 

. . R 
noint a . 

~T.? 

Of course, the results of figure 7.6 only give an order-of- 
ragnitv.de estimate of the radiation- convection heating coupling, 
riot included are the effects of variable transport properties, 
enthalpy grndient at the edge of the boundary layer, and differences 
in the character! Stic Reynolds and Prandtl numbers between the 
radiating and nonradiating cases. Also no account has been taken 
of the effect of radiation in the boundary layer. In the cooled 
region of the boundary layer adjacent to the vail the gas will 
absorb more radiant energy than it will emit. This will tend to 
increase the slope of the enthalpy distribution adjacent to the wall 
thereby increasing the convective heat transfer somewhat. 

The effects of radiation cooling, gray absorption, and spectral 
absorption on the ratio of convective heating rate for a radiating 
shock layer to that for a nonradiating shock layer, h , is shown in 
figure 7 . 7 • It is apparent that radiation cooling plays the major 
role while absorption (both gray and nongray) tends to reduce the 
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effectiveness of radiation cooling. The calculations for curves 
1 and 2 ignored cooling. Consequently, no reduction in the cal- 
culated value of the convective heating rate was obtained. Curve 
3 includes radiation cooling and ignores absorption. Thus the 
reduction in the calculated value of the convective heating rate 
is maximized in this approximation. Finally curves 4 and 5 
indicate that absorption inhibits the effectiveness of radiation 
cooling, and since absorption is more important in a nongray gas 
than it is in a Planck equivalent gray gas the rate of convective 
heating will be greater in the nongray case. 


E. The Role of the Radiation Cooling 
Parameter and the Bouguer Number 
The radiation cooling parameter e admits of several physical 
interpretations which are useful in the understanding of the radi- 
ating shock layer. Of these, one of the most useful is the 
following: 


€ = 


I rate of emission from 
l element of volume of gas 
Vemerging from shock 


/time required by element of 
/ volume to traverse distance 
l A. at rate of emergence from, 
Vshock 


2 (energy of element of volume upon emergence from shock) 


Here A^ is the shock standoff distance in a nonradiating (or 
adi abati c ) flow . 

It can be seen from this interpretation that the radiation 
cooling parameter is indicative of the slope of the enthalpy 
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distribution immediately behind the shock. In fact, in the trans- 
parent limit there is a direct relation between e and the initial 
slope. That is, 


ah/dfl/Ta) = 2e 

(see, for example, chapter IV, figure 4.2) 

In the case of an optically thick shock layer, the initial 

enthalpy gradient is reduced by absorption. However, a lower bound 

to the gradient is the value € (one-half the transparent value) 

because the emergent elemetary volume will emit at least twice as 

much energy as it absorbsj it emits energy at a rate proportional 

to the Planck function at the equilibrium shock temperature, T , 

s 

in both the upstream and downstream directions while it absorbs 

energy at a rate at most proportional (by the same factor; the 

monochromatic volume absorption coefficient) to the Planck function 

at temperature T from only the downstream side, 
s 

A physical interpretation of the Bouguer number is given below: 



(radiationless shock standoff distance, 

''Planck average photon mean free path in gas emerging 
. from shock 


Only when conditions do not vary greatly across the shock layer will 
the Bouguer number be indicative of the Planck mean optical thickness 
and only when the gas is nearly gray will the Planck mean optical 
thickness be indicative of the various important monochromatic 
optical thicknesses. Consequently, critical values of the Bouguer 
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number are subject to a number of influences; among them, the 
enthalpy and spectral variation of the absorption coefficient and 
the value of the radiation cooling parameter. For example, the 
value of kp for which absorption first becomes important is 
about 0.1, when the radiation cooling parameter is very much less 
than one and the absorption coefficient is independent of wave- 
length. When the absorption coefficient varies spectrally as 
shown in figure 3.10, chapter III, and when e is very small, 
absorption begins to become important for Bouguer numbers as small 
as 0.001. With e about 10 for a gray gas absorption is important 
for values of the Bouguer number greater than about three. Dispite 
these drawbacks, the Bouguer number as defined in this investigation 
is about the best a priori indicator of the importance of absorption 
that can be obtained. 

When the radiation cooling parameter e is very much less than 
one, an elementary volume of gas vill lose very little of its energy 
by radiant emission in the time required to traverse most of the 
shock layer (of course, it takes an elementary volume of gas 
travelling along the stagnation streamline an infinite time to 
reach the wall) . Hence, radiation cooling of the shock layer will be 
slight. When the radiation cooling parameter is very much greater 
than one, an elementary volume of gas will emit energy at such 
a rapid rate that the energy of the volume will be reduced a 
significant amount before it leaves the vicinity of the shock. 
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This is true whether the shock layer is optically thick or optically 
thin (that is, regardless of the size of the Bouguer number) . This 
physical argument is used to establish the exi stance of the thermal 
boundary layer behind the shock in the radiation depleted shock 
layer (chapter VI) . If the shock layer is optically thick, the 
reduction in enthalpy will continue only so long as the elementary 
volume is within about a photon mean free path of the shock. Beyond 
this point, the elementary gas volume receives radiation from all 
sides and begins to establish a condition of radiative equilibrium 
with its surroundings . The energy lost during the time required by 
the elementary volume to travel a single photon mean free path is 
characterized by the ratio of the radiation cooling parameter to 
the Bouguer number, € (it was shown in chapter V that the 
enthalpy level in the interior of an optically thick shock layer 
was characterized solely by the parameter, e J kp) . 

Within the interior of an optically thick shock layer, radiation 
heat transfer can be treated in a manner analogous to conductive heat 

transfer. Thus, one would expect that a parameter analogous to the 

* 

Peclet number could be constructed which would suggest the nature of 
the enthalpy boundary layer adjacent to the wall. Such a parameter, 
which is a ratio of the importance of convective to radiative heat 
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transfer is given by the grouping k?/e.* Since the thickness of 

-l/2 

the enthalpy boundary layer is characterized by Pe ' in the 
conduction problem, one expects, by analogy, the thickness of the 
enthalpy boundary layer adjacent to the wall in an optically thick 

. The 

importance of this parameter (in a somewhat different form) and its 
analogy with the Peclet number was pointed out previously by Goulard 
(ref. 21). 

The importance of the surface reflectivity, r , depends on the 
importance of absorption in the shock layer. When absorption is 
negligible, the effects of surface reflectivity are negligible 
because all photons originating within the shock layer will escape 
the layer and it matters not whether some of these photons are 
absorbed by the cold wall or reflected by the wall into the free 
stream. However, when absorption is important, the reflected photons 
have a large probability of being recaptured in the shock layer. Thus, 
an increase in surface reflectivity tends to raise the enthalpy level 
of an absorbing gas in the vicinity of the wall . 

In this section it was shown that both the radiation cooling 
parameter e and the Bouguer number kp play prominent and inter- 
related roles in determining the character of the radiating shock 

In the optically thick shock layer analysis of chapter V the 
gas was assumed to be gray. Thus the fact that the Bouguer number 
was based on a Planck mean absorption ceofficient was of no con- 
sequence. However, in the case of a nongray gas, it would probably 
be more correct to use a Bouguer number based on a Rosseland mean 
absorption coefficient. 


radiating shock layer to be characterized by 
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layer. Further, it was shown that the spectral variation of the 
absorption coefficient greatly influences the role of the Bouguer 
number. In general, then, one cannot ignore either of the processes 
of radiation cooling and nongray absorption. 



CHAPTER VIII 


SUMMARY AMD CONCLUSIONS 

A mathematical model for the stagnation region of a radiating 
shock layer was derived in this investigation subject to the following 
conditions: (l) the gas in the shock layer is in local thermo- 

dynamic and chemical equilibrium, (2) the body geometry is axisym- 
metric , ( 5 ) there is no mass addition to the flow from the body 
surface, (4) the thicknesses of the shock and the viscous boundary 
layer are small in comparison to the shock standoff distance, and 
( 5 ) absorption in the free stream ahead of the body is negligible. 

The divergence of the radiant flux vector, which appears in the energy 
equation, was formulated to include a wavelength varying absorption 
coefficient. The body surface was considered to be cold and to reflect 
diffusely and independently of wavelength a fraction r of the 
incident radiation. The results of a boundary layer analysis indicate 
that the equations for the flow in the inviscid region are independent 
of the boundary layer equations only when the boundary layer is 
optically thin or optically thick. It has been assumed throughout 
this study that the boundary layer is optically thin. Simple corre- 
lation formulas for the thermodynamic and optical properties of high 
temperature equilibrium air were developed and used herein. 

The general form of the governing system of equations was found 
to be integrodifferential in character. The solution of this system 
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is extremely difficult to find even with numerical techniques and high 
speed electronic computing machines. The approach of this investigation 
was to take advantage of the simplified form to which the governing 
equations were reduced when the radiation cooling parameter e and 
the Bouguer number kp took on limiting values and obtain approximate 
analytic solutions if available. It was found that the general problem 
reduced to a singular perturbation problem in each of the four cases 
studied. A small perturbation solution valid when the energy lost to 
the shock layer by radiation is small (i.e., when the radiation cooling 
parameter is small) is described in chapter III. The Poincare- 
Lighthill-Kuo perturbation of coordinate method was used to obtain a 
uniformly valid solution. This solution was used to study radiation 
cooling, absorption, effects of surface reflectivity, and effects of 
nongray optical properties. 

An optically thin shock layer method of solution, discussed in 
chapter IV, utilizes an expansion in terms of the Bouguer number kp 
to reduce the governing system to purely differential form. Again 
it was necessary to resort to the P-L-K method to obtain a unifo rml y 
valid solution. This solution was used to study radiation cooling, 
absorption, and the effects of surface reflectivity. 

The optically thick approximation, valid when the optical 
thickness of the shock layer is very large (i.e., the Bouguer number 
very much greater than l) was used to obtain the solutions of 
chapter V. The governing equations were reduced to differential form 
through the use of a substitute kernel approximation. Two thermal 
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boundary layers were seen to exist; one adjacent to the shock and 
the other adjacent to the wall. It was noted that the Ross eland 
approximation together with a properly specified temperature jump 
or slip condition at the wall reduces the governing equations to the 
same form as the substitute kernal approximation in the interior or 
isenthalpic portion of the shock layer and in the wall boundary 
layer. However, the Ross eland approximation with slip conditions 
was found to be inadequate for analyzing the shock boundary layer. 

The optically thick solutions were restricted to gray gases but were 
used to study radiation cooling, absorption, and the effects of 
surface reflectivity. 

The radiation depleted shock layer was analyzed in chapter VI. 
This approximation is valid when the rate at which energy is radiated 
away from the shock layer is nearly equal to the rate at which energy 
enters the shock layer ( i . e . , the radiation cooling parameter is very 
large) so that the enthalpy level is very much less than the radia- 
tionless value. The substitute kernal approximation was used to 
reduce the governing system of equations to differential form. The 
method of matching of inner and outer expansions was used to obtain 
solutions valid in the thermal boundary layer adjacent to the shock 
and in the interior of the shock layer. These solutions were 
restricted to gray gases but were used to study radiation cooling, 
absorption and the effects of surface reflectivity. 
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It is apparent from the results presented in chapters III 

through VI, that radiation cooling first becomes important when the 

rate of energy lost by radiation from the shock layer is only about 

1 percent of the rate with which energy enters the shock layer. 

Absorption in a gray gas begins to become important for shock layer 

optical thicknesses greater than about one-tenth. An increase in 

the surface reflectivity r from zero reduces the radiant heat 

w 

transfer by a factor of roughly 1 - r , and increases the heat 

w 

transfer rate to the wall by conduction because of an increase in 
enthalpy level near the wall. 

The results of some nongray calculations are presented in 
chapter III. The Planck mean absorption coefficient can be used to 
compute the enthalpy distribution and the radiation heat transfer 
rate to the wall as long as the optical depth of the shock layer is 
very much less than 1 in all wavelength regions in which a significant 
amount of radiant energy is emitted. For larger optical thicknesses 
nongray effects are very important. 

The various approximate solutions were used to compute the rate 
of radiant heat transfer to the stagnation point of blunt objects 
traversing an optically gray model earth atmosphere. The results of 
this computation indicate that at every altitude and velocity there 
is a finite value of body nose radius for which the rate of radiant 
heat transfer to the stagnation point is a maximum (this result is 
contrary to the earlier results, based on the assumptions of an 
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isenthalpic and transparent shock layer, which indicated that the 
heating rate was directly proportional tc nose radius). A signifi- 
cant reduction in the computed value of the radiant heating resulted 
upon taking the nongray character of air into account. This served 
to emphasize that the nongray character of gases plays a very real 
and important part in problems of radiation gas dynamics. 

In general, the coupling between radiant and convective heat 
transfer is such that increases in the rate of radiant heat transfer 
result in decreases in the rate of convective heat transfer to the 
body surface. Of course, the amount by which the total heating rate 
is affected cannot be determined from this inviscid analysis. 

It is hoped that the methods used in this investigation will 
point the way to simplified methods for treating the general problem. 
For example, the study of nongray absorption coefficients by means 
of the small perturbation method may lead to the difinition of an 
approximate mean absorption coefficient through which the general 
nongray problem can be reduced to an equivalent gray problem. As 
was pointed out previously (chapters V and VI) the integrodifferential 
system of governing equations for gray gases can then be reduced to 
purely differential form through the use of the substitute kernal 
method or other available methods (see for example, refs. 3k and 53). 

Obviously such simplifications are urgently needed if current 
analyses are to be extended to include the important effects of 
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chemical nonequilibrium, absorption in the free stream ahead of the 
shock, and the injection of foreign species into the shock layer 
due to ablation of the body surface. 
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APPENDIX A 


THE VISCOUS BOUNDARY LAYER 


In this appendix, a boundary layer analysis will be performed 
on the integro-differential system (2.51) to (2.58) to determine the 


x ux iu wx uaic c^uauxuaS mix mic xuviduiu a crgxOIi emu. uxjlc viSCuUS 


boundary layer and to determine under what conditions such a sepa- 
ration can be achieved. For convenience, the system will be 
rewritten here 


f(n) h'h) + * 2 p^h) h'(il)] + €I[ti] = 0 (Ai-l) 

27? [?*(h) f'*(ti)] + 2f(T]) r(tl) - ff(Ti)] 2 + a 2 h(Ti) = 0 (A*2) 


f(0) = 0 (A- 3 ) 

f'(0) = 0 (A-4) 

f(n A ) = 1 (A-5) 

J = - y~= = = (A-6) 

h(0) - h w (A-7) 

h(n A ) = 1 (A-8) 
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where 

€ = rpg *p a a (a-9) 

s 

3**(h) = Pr g 3> 2 (h) (A-10) 

and A, introduced here for convenience of notation, is the inverse 
square root of the Peclet number. 

When the parameter A is very much less than one, a per- 
turbation type solution can be attempted. However, the energy 
and momentum equations each lose the most highly differentiated 
term as A vanishes. As a result, neither the zero-order (in the 
small parameter A) solution for f(q) nor that for h(rj) can 
satisfy all the boundary conditions. In particular, the conditions 
f*(0) * 0 and h(0) = 0 must be relaxed, and the perturbation 
solution will not be valid as tj approaches zero. Thus, this 
pi~obl6in is a singular perturbation problem of the "boundary layer" 
type (refs. 36 and 5^) • 

In order to obtain the boundary layer form of the equations, 
the "stretched" coordinate | = A ^ is introduced where a is an 
as yet undetermined constant. It is also convenient, to avoid 
confusion, to introduce the change in notation 



♦ft* 

II 

^5 

(A- 11 ) 

0 

g'U) * f*(n) 

(A- 12 ) 


j[i] . i[i] 


(A-13) 
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(A- 12) is v/ritten in this particular form because it is f ' (tj) and 
not f(n) which fails to satisfy the boundary condition at t] = 0. 

When the stretched coordinate t and the definitions (A-ll) 
to (A-13) are introduced into system (A-l) through (A-8), the 
only choice for a which will retain the most highly differentiated 


terms without loss of the most significant terms in the "unstretched" 

p 

problem is a = 1. Thus, A and not A is the significant small 
parameter and the stretched coordinate is 


5 = A _1 r) 


Perturbation solutions are now sought in the forms 

00 

= V A n z n (fi) 


n=o 


(A-lk) 


(A- 15) 




\ Tv" K (I) 

L . i 

n-o 


_ n ■ 


in the boundary layer, and 

00 

h(n,A) * £ A n h n (n) 


n=o 


f(n,A) A n f n (n) 
n=o 


(A-l6) 


(A-17) 


(A-18) 


in the inviscid region. 
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It shall be assumed that all functions of h (and i) are 

analytic about the value h Q (and i 0 ) so that they may be 

expanded in Taylor series about h =. h and i = i in the 

o o 

following manner: 


F(h) = F Jh Q + + A 2 h 2 + ••• j = F(hj + AF^h^J^ + A" 


FA 


+ Ito + ••• 


The exi stance of the expansions 


00 


M - 

L X ” T nM 


n=o 


rv* 

I \ J „H 


n=o 


00 



n^o 


o) h 2 

(A-19) 


(A- 20) 


(A-21) 


a 


(A-22) 



225 


v - y v (ju23) 

&> 

is also assumed without, for the present, specifying details of the 
terms I n [.iQ and J [s]. 

Furthermore, to insure compatibility of the boundary layer and 
inviscid solutions, it is necessary that the inner boundary con- 
dition on the outer solution be written in the form 

f(5) * 0 (A-24) 

where 6 (the displacement distance) is specified by the matching 
condition 

lim g(|) * A" 1 f(rj) (A-25) 

00 

The quantity 6 depends on A and must be written in expanded 
form 

A” 5 n (A-26) 

The term 5 was chosen to be zero because 5 is order A. 
o 

The system which describes the solutions valid in the inviscid 
region can be obtained by substituting expansions (A- 17) - (A-23) , and 
(A-26) into system (A-l) - (A-8) . The result is an infinite power 

series in A the sum of which is zero for all values of A. The 
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only such series is one for which the coefficient of each of the A n 
terms is identically zero. These coefficients yield a set of recur- 
sive integro- differential systems. The system of zero order is 


f 0 (n) h£(i) + eijrj] = 0 

(A- 27 ) 

2f 0 <i> «J(i) -[f;<n)] 2 ♦ 4 h 0 (n) " 0 

(A-28) 

f o (0) . 0 

(A -29) 

1 

(A -30) 


(A-31) 

h o(\) * 1 

(A-32) 


The system which describes the solutions valid in the boundary 
layer can be obtained by substituting expansions (A -15), (A-l6) , and 
(A-19) - (A-23) into system (A-l) - (A-8) . As for the inviscid 

case, this procedure results in a set of differential systems. 

The zero- order system is 

O 5 '] + «o< 6 > J o (5) + M 0 W *° 


(A-33) 
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2 [^( l o) <<»]' + - [s;(t>] 2 * 4 

i 0 (i) - o 


(A-54) 

g o (0) * 0 

(A-55) 

6q(0) = 0 

(A-3 6) 

lim g'(g) = f'(0) 
!-»» 0 ° 

(A-37) 

II 

o 

(A-38) 

lim l (|) = h (0) 
| -»00 ° 

(A-39) 

iii ww~yiuv.i' i/ci'iu U i wHg uxo^iav-ciucii u uxSb&UlUC 0 ^ y IS 

found from the matching condition 


g o< s > - ( 6 - 8 i) 

(A-40) 


It is apparent that the zero-order solution for the boundary- 
layer equations depends only on the inviscid enthalpy level in the 
vicinity of the wall and not the enthalpy gradient. The enthalpy 
gradient will, of course, have an effect on the first order boundary 
layer solution. Thus, if the enthalpy gradient is very large, as it 

can be for a radiating shock layer, the boundary layer solutions must 
be carried out to first order in A. 
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The divergence of the radiation flux (see eq. (2.63)) includes 
integrals which extend over the whole domain of the problem. It is 
convenient to separate each of these integrals into two integrals 
as follows: 


\ A 

/ 0 * ■ M) 4t * ■ 1 / ' **0“ 


\\r * J * ">n u ?v]rn ti’x “Air a 

t a 

(A-Ul) 

# 

0, 


f o ■ Nr 


and 


T ?v T a 

f o A b 4 t ;\) dt A 0 f r * A b a h ( t ' A )] E2 ( kp **) 


dt., 

/\ 


( a«U <? \ 


+ A 


u 

f B >T i ( , *xTK( k p N ss >. 

' O I v y | / 


where 


< ■ *>* - * / V 1)d6 


(A-i» 3 ) 


* 

is the monochromatic optical thickness of the boundary layer, i is 
the thickness of the boundary layer in terms of the stretched 
Dorodnitsyn variable |. 
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It is convenient to redefine the monochromatic optical path 
length as 



In order to expand equation (2.65) as a power series in A, it 

is necessary to expand the exponential integral functions and all 

functions of h (and i) as well. Expanding the optical thickness 
. * 

yields, for t. 1 t 



(A-!)6) 


(A-47^ 
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The exponential integral functions can be expanded in the Taylor 
series 

E n (x - y) ■ eT(* 0 - y 0 ) + ^ - y,) + • • .1 

(A-U8) 

a E n ( x o ' y ») ' ^ ' y l) E n-l( X o * y o) * ••• 

If the argument is order A 


EJM = E n (0) - ^ xE^tO) + ... 


n-2 


, t \n-2 ^n-2 x' -n t*\ 
* ( - l) X (n - g JT E S (0) 


- (- 1 )”' 1 t*”' 1 !nX ) Tir^TTT 


(A-49) 




n-1 


where y is Euler's constant (y = 0.57721b). Use of this expansion, 
while it avoids any dependence of the terns I n £tJ on A, introduces 
terms of order A in A into the boundary layer solutions. 

Incorporating the various expansions into equation (2.63) and 
separating the result into powers of A and A in A yields the 


zero-order expressions 



- 2 “pC h o<’'0 B [ h o('> ) ] 





231 



t 

cd 


i 



The second of these expressions, which contains only definite 
integrals, is valid only when expansion (A-ho) holds. But 
equation (A-43) converges in the first few terms only if the 
arguement Toe (or in the terns of this problem A kp a) is small 
compared to 1. Thus, expressions (A-U 9 ) and (A-50) can be used 
only when the boundary layer is optically thin, that is, 



APPENDIX B 


THE METHOD OF SMALL PERTURBATIONS - 
MATHEMATICAL DEVELOPMENT 

In this appendix, the method of small perturbations is used to 
obtain a solution to the integro-aifferentiai system of equations 
governing the flow in the inviscid region of a radiating shock layer. 
Mathematical details which are not considered to be appropriate to 
the main test (chapter III) are included herein. 

The system of equations to be treated are presented below. 


f(t)) h'(n) + ei[n] = o 

(B-l) 

2 f(n) f”(r)) . [r(r,)] 2 + a 2 h(n) * o 

(B-2) 

H\) “ 1 

(B-5) 

O 

It 

o 

(B-4) 

* 1 

(B-5) 

f ' w * ^<rn o 

(B-6) 
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The conventional perturbation procedure .- If it is assumed that 
the functions h( e) and f(n;e) are analytic in the vicinity of 
6 * 0 they may be written in the expanded form 

00 

^ e n h n (rj) (B-7) 

n=o 


f(nje) = £ e 11 f n (n) (B-8) 

n=o 


It is anticipated that the first few terms of these expansions will 
provide an accurate estimate to the solution of system (B-l) - (B-6) 
when the parameter e is small compared to unity. 

The integral terra I [t|] and the constant also depend on 
the parameter e through their dependence on the functions 
h(T;€) and f(T);e). These quantities will also be assumed to be 
analytic functions of e near e = 0, so that 

00 

I[nje] = £ e n I n [tj] (B-9) 

n=o 


00 



(B-10) 
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Substituting the expansions (B-7) - (B-JO) into system (B-l) - 
(B-6) yields 


( f o h ;> + € {w ♦ vi + yi} 


-2 L, ... 


+ fc r 2 n o + r i°i + x o n 2 * J ’iL T1 -l + •*' * 0 


k f o - K] 2 + a ' h 


+ a 


c] + e ( 2f o f i 

h l} + ^o r 2 


2f’« + 2f"f, 
o 1 o 1 


- 2f'f" + 2f"f 
2 o 2 o' 2 


+ 2f f 


"x - Pi] 


+ \ h °(\)} 


f>(\) • x } + 4 w 
+ 4 2 w + \ %\) 


(B-ll) 


(B-12) 


+ t) h' n A i 

*2 °OoJ 


(B-13) 


{ f 0 (°)} + K f !(0)} + e2 { f 2 (0) ) + “* 


0 


(B-14) 
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(MX) - * VK\)} 

+ £ t f <\) * \ fi (\) + \ f i(\) 

+ 2 \ f 'o{\)J+ ••• - 0 (B-15) 

k \) ■ e fe\) + \ %'4 
+ £2 k\) * \ %\) ♦ \ *;(\) 

<E-l6) 

Since the small parameter e is arbitrary system (B-ll) - 
(B-l6) can be satisfied only if each coefficient of each expansion 
in e is identically zero. This leads to a recursive set of 
purely differential systems. 

The zero- order system is 


h* =s 0 
o 


2f f" - rff + a 2 h 
o o L oj c 


= 0 


(B-17) 


(B-0l8) 
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M\) • 1 

f o (0) = 0 

f A ) = 1 

V o / 



a 

V2X(l - X) 


(B-19) 

(B-20) 

(B-21) 

(B-22) 


The solutions to this system are easily found with the result 


h = 1 
o 


(B-23) 


f * (l - a) Tj + a rj 


(B-24) 


TL = 1 

% 


fB-25) 


The systems of first and second-order may be written in the 
general form 


f h* + f.h* + I = 0 
on 1 n-1 n-1 


(B^26) 


f f" 
o n 


f'f + f ,: f 
on on 


<b ( tj ) 

n ' ' 


(B-27) 


h (1) . 0 
nV ' r n 


(B-26) 
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f n (0) * 0 

(B-29) 

”(\) ‘ 0,1 

(B-30) 

»(\) • *» 

(3-51) 


Equation (B-26) can be integrated directly to obtain 


h (n) - i ♦ f 1 ^ VA (!t) 

n ^n J n fJZ] 


dx 


(B-32) 


The abriged version of equation (B-7b) admits the pair of linearly 
independent solutions tj + a/2(l - a) and tj 2 . Following Ince 
(ref. 55) the complete solution is found to be 


o i _ 1 & (x) dx 

- 1 [a(i - a) n ♦ a]/ — s- -5 

- O (_v ju - aj X T aj 


1 „2 f' 1 2(1 - a) X - a * , x 3 

- o ’l / -37=* L o Mx) *c 

2 '-'r, x 2 [(l - a) x + a] 2 n 


I < 2 * “> fj 


2 Z' 1 4,6 

o Ql - a) x + aj 2 


(B-35) 


Substituting this expression into condition (B-3l) provides a 

relation for the determination of tj. , that is 

A n 


(B-J4) 


r l 4 (x) dx 

29 + / 2 • « t 

n [(1 - a) x + a] 2 n 

The quantities $ n (T}), 9 n > and i|r are 


V 1 *) * ■ 1 

(B-35) 

® 2 (n) • - 2 + [q] 2 - 1 »%(!!) 

(B-36) 

o 

II 

(B-37) 

rv^ 

II 

Jf 

H 

O 

I 1 

(B-38) 

\ ■ - (2 - a)^ 

(B-39) 



* 1 (x) 1 

[(1 - a) x + a] 2 J\ 


(2 - a)^ - (1 - a) ^ 


(B-40) 


« - 2 (1 - a)tj 


(B-4l) 



2 k 0 


\|f 2 = -|2(2 - a)^ + (2 - aj 


pi Mx) dx *\ 

’j 0 f(l - a) X 7 7pj\ ' 2(1 ' a ’\ 


[(l - a) x + a] J ^1 


(B-42) 


The divergence of the radiant flux l[q] may be written in 
expanded form by substitution of (B-7) and the expansions of the 
quantities K^(h), B^h), t^(tj$€), ana E^jj^TjjejT) into expression 
(2.86) of the text. For completeness the expanded forms of k ^, B^, 
and E are written down here. 


K^(h) s k^(i) + e k^(i) h^n) + e 2 ^(l) h 2 (l) 

+ ^ *^(1) h^(q)J + (B-4-3) 

B x (h) = B^(l) + e B^l) lu(n) + e 2 ^!) h 2 (Ti) 


+ | B^d) hj(n) 


(B-44) 


r n 

= hp J K A (n) dq 


ak F 


^(i)n 


+ e 



(B-45) 
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E n0\^ " T A^3 ~ E n( lc P [ K ^ l) ^ ‘ + € ^( 1 ) J V x)dx + “3) 

- E n|> V 1 ’ <" - 5 >] 


eV * ( 1 ) 

**P ”A'“' 


\J | 


Vi fv^Hv 
1 '“ - * 


V. Hr ^ (l \ { n 

~n-lLTF 


ij] + •• 

(B-46) 


The following property of the exponential integral functions was 
used to obtain (B-46) 


E „-1 W - • s V*> 


With these expansions in hand, expressions for the terms 
and can be obtained. The results are 








(B-47) 
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(B-:48) 


In these expressions, the notation has been simplified somewhat 
by omitting the argument 1 in the terns B^, and and 

by introducing the quantities 


lt 7v = l! P ,t X 


(B-l *9) 


r o x * r v C 1 - 2 (B ’ 50) 


/ c M 5)e 2 M 35 


+ b> ijka r h, ( 6 ' d I + n A 

-L /\ /\ c.\ /\yj 


* \ /' 'V yj 


(B-51) 


The P - L - K solution .- It has been pointed out in the text 
that the first order solution for the enthalpy distribution has a 
logarithmic singularity at the point tj = o and the second order 
solution behaves like the logarithm squared. As a consequence the 
assumed expansion diverges as the origin is approached and the 
small perturbation solution is not uniformly valid. In order to 
obtain a solution which is uniformly valid throughout the domain of 
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the problem, the Poincare - Lighthill - Kuo method (see ref. 47) 
will be used. In this method the independent variable as well as 
the dependent variables is expanded in a McLauren series of £ . 

For this problem 

T = x + e tj*(x) + e 2 tj*(x) + ••• (B-52) 

h(nje) = hj(x) + e hj(x) + e 2 h*(x) + ••• (B-55) 

The superscript * has been used here to distinguish between the 
coefficients in the P-L-K expansion and the coefficients in the 
conventional expansion (equation (B-7)). The quantities f(Tj) 
and l[n] may also be expanded in terms of x as follows: 

00 00 

f(n) = ^ € n f n (n) ■ ^ 6 n + eq*(x) + •••] 
n=o n=o 

- f Q (x) + e[f 1 (x) + n*(x) r(x)] + e 2 [f 2 (x) + t)*(x) fj(n) 

+ t' Q (x) + | q* 2 (x) r(x)] + ••• (b-54) 

Similarly, 

I [ T l] - y*] + € { I lL x l + V x > y*]}* (B-55) 

When expansions (B-52) - (B-55) are substituted into equation (B-l), 
a set of equations for h*(x), h£(x), h*(x), and so forth result. 



The quantities n*(x), ^(x), and so forth and their first derivatives 
also appear. These quantities are arbitrary and should be chosen in 
such a manner as to reduce the strength of the singularities in the 
higher order terms, h*(x) , so that these singularities are never 
stronger than that of the lowest order singular term (in our case 
the first order term) . Pritulo (ref. 48) has shown that the 
coefficients in the expansion of h(Tjje) in the P-L-K method are 
related to the coefficients of the conventional expansion in the 
following manner: 


h *( x ) a h o (x) (B-56) 

h*(x) = h^x) + tj*(x) h^(x) = h^x) (B-57) 

h|(x) = h 2 (x) + TJ*(x) h£(x) + TJg(x) h^(x) + | TJ*(x) h”(x) 

= hg(x) + TJ* h^(x) (B-58) 

Now, instead of choosing differential equations for the T)*(x) in 
order to satisfy the criterion previously states, one can choose the 
values of the tj*(x) directly. In this case, an obvious choice is 
simply* 

* * 

This choice satisfies the condition tj^(1) = t] . 


* / X 

which gives h g (x) 5 0. 

The result of transforming the independent variable by means of 
formula (B-52) is to remove the singularity from the domain of the 
problem. That this is true can be seen by noting that the condition 
tj = 0 does not imply x = 0 but rather (for this problem) implies 
that x is some small positive number x q . Hence, to first order 

'ft 

h(rj; e) = 1 + €h^(i}) is nonsingular in the domain 0 $ t| 5 !• 



APPENDIX C 


OPTICALLY THIN SHOCK LAYERS - 
MATHEMATICAL DEVELOPMENT 

With the substitute kernal approximation the divergence of the 
radiant flux vector can be written 


**,>*,> - 


+ r e 
w 




/ ^ (I) V 


l)e 


-2kpT A ( |) 


(C-l) 


d|> dA 


The monochromatic optical path length kpT^(ii) is given by the 
expression 


k p T ^( T l) - kp / 


(C-2) 


where kp is the Bouguer number 


kp =P s Kp a a (C-3) 

s 

The approximate governing system presented in chapter IV is 


f(Tj)h'(n) + «i[Y]« o 


2f(ri)f"(n) - 


f(ll) 


+ a 2 h = 0 


(C-U) 

(C-5) 
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f(0) * 0 

f (^ ■ 1 


V A J 1 + j2X(l - X) 

f 

h ( n A) ■ 1 

The quantity h is defined "by the expression 



h(Odft 


2k7 

(C-6) 

(C-7) 
(C-3 ) 

(C-9) 

(c-io) 


The conventional perturbation procedure .- If the functions 
h( tj j kp) and f(rj; kp) are analytic in the vicinity of kp = 0*, 
they may be written in the expanded form 

00 

h('ii; kp) = y kj h n (t]) (C-ll) 

n^O 

00 

f(n; kp) = £ h n (rj) (C-12) 

n=0 

It is anticipated that the first few terms of these expansions will 
* 

It is assumed herein that such is the case. 
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provide an accurate estimate of the solution to the system (C-4) to 
(C-10) when kp is small compared to unity. 

In addition, all quantities which depend on the parameter kp 
either directly or indirectly must be expanded in terms of kp. For 
example, a function 3Lh(Ti)] 

becomes 


S[t>(n>] * 3[\(n) + kj h^) . ...j 

= + k p ^k(n)]M^) + ••• (c-u) 

= :* o 0l) + kp ^(qjh^q) + ... (C-l4) 

The quantities £ o (t)) and ^ 0 (t) have been introduced to simplify 
the notation. The constant t) is given by the expansion 


%= \ q + \ + **• 

Substituting the expansions (C-ll) to (C-15) into system 
(C-4) to (C-10) yields 


(C-15) 


{< 


f 0 (n)h^(q) - 2c Kp (t))B o (t)W 


+ + f ]> )h o (ri) * 2€ |\. (n) B 0 (n) 


+ ^ (ti)b o (ti r ^(n) 
° J ^ 


+ 26 f 1 + r J f K (n) f ° \ (6)B> U)di]d7^+ ... - 0 
' ' ^ O o “ O o O J 

j -» £\ 


(C-16) 


2k9 


,2f 0 (i)f"(n) 




+ a hA + fcj > <2f 0 (ri)I-"(n) 


- + sr^CnJf^tn) + a 2 hA + 


... = 0 


f . (0) + k^f,(0) + ... = 0 
o r x 



(C-17) 

(C-13) 

(C-X9) 

(C-20) 

(C-21) 

(C-22) 


Since the small parameter is arbitrary system (C-l6) to (C-2l) 
can be satisfied only if each coefficient of each tern is identically 
zero. This leads to a recursive set of purely differential systems. 
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The zero -order system is 


- 2€Kp (ri)B o ( q) = 0 


2f 0 (r 1 )f; , (a) - fjOl) ‘ 


+ a 2 h =0 

O 


f o (0) = 0 


f o (\> =1 

o> 


f ' tj a 
o 'A 


oj 1 + y’2X(l - X) 


h o hA ) = 1 


h = 


_i_ r 

\ J 

o o 


MDas 


The solutions to this system are 


t 2en A 

dh _ In (1 - a*)x + a* 


^(hjBCh) a 


f 0 ( T l) 


= 1 





+ a f^o \ [■£- 


0\\ 


= (l - a*)x + a* 


(C-25) 

(C-24) 

(C-25) 

(C-2o) 

(C-27 ) 
(C-28) 

(C-29) 


(C-50) 


(C-51a) 


(C-31b) 
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1 + Y2X(1 - X) 

i + /a 0 x(i - x) 



where the def ini lions 



o 


x = 



have been added to help simplify the notation. 
The first-order system is 


f 0 ( 'l)h.[(T}) - 2e K p (rj)B o (rj) + Kp jh^ tj) 


L o 


(C-32) 


(C-55) 


(C-34) 

(C-55) 


= - f 1 ( i i) h o( T i) ■ 2e f i + r w ) r \ (i) r ° \ u)b. (e)dft 

^ '“'O O Jo o 0 


dX 




(c-56) 

f 0 (1)f{(1) - 

+ ^(o^Cn) + | a 2 h x = o 

(C-57) 


f x (o) = o 

(C-58) 

f 4\) 

, ^ 2 \ 

= -n. f '[ TK ] = 

*1 °l A oJ 1 + y 2 X(l - X) 

(C-59) 
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2(1 - 


f i '( v ) - ’ V ° l \ 


2 


(CJO) 



h if n. 


A ° J = " Vh\j 


= -26 T] 


'*1 


(c-M) 


The solutions to this system are 

r 


^(■l) * -2e*p (x)B 0 (xKtj^ + 



A 


a \ ^ 

1? 

1 

H 

1 + 1 

/2x(l - x)j 

(l - a*)x + a* 


1 (1 - a*)x + a* 

- — In 


+2e 1 + r 


w 


«p WB,(*>4 f 

D r\ v v 




O w X 


f 0 (fi)Kp (i)B o U) 


dl 


f 2 (x) = 2 — i x 


k + 


1 - 


yK ** 1 - 


X) 


\i + y2h o x(i : x ) J 


(C-42) 

(C-45) 




" 4 lT ~ \ 

o 


(C-JA) 


h l = 


- r 1 

h o J h i( x ) ds 


h 0 + T ( 1 ■ h o) 


(C-45) 
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The P-L-K solution. - It can be seen on careful Inspection of 


equation (C-46) that the first-order terra h. (x) displays a singular 
behavior in the vicinity of the wall (x = 0)? Consequently, the 
assumed expansion for h(x; kp) diverges as the origin is approached 
and the perturbation solution is not uniformly valid. However, if 
che coordinate x is perturbed, the solution can be made uniformly 
valid. Thus, according to the P-L-K method, 

x = y + kp xj(y) + ... (C -*£) 

where y is the transformed variable. The enthalpy when expanded 
in terns of k, 0 with coefficients as functions of y, not x, 
becomes 




-X , kp J 


b D (,-) 






and the nondimensional stream function 


(c-49) 


f(x; kp) = f*(y) + kp f*(y) + ... (C-50) 

According to Pritulo (ref. 1*8), the coefficients in the P-L-K- 
expansions can be related to the coefficients in regular expansions 
in the following manner 

h*(y) = h Q (y) (C-51) 

hj(y) - hjy) + xj(y)h£(y) (C-p2) 

* 

See chapter IV for a more complete discussion of this singularity* 
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f*(y) = f 0 (y) (c-55) 

f*(y) = f x (y) + xj(y)f ’(y) (c-5ii) 

The arbitrary quantity x*(y) should be chosen so as to eliminate 
the singularity in h*(y). An obvious choice is 

x x (y) = -^(yj/h^y) (C-55) 

The transformation of the independent variable by means of 
formula (C-43) removes the singularity from the domain of the problem, 
because y takes on some small positive value when x is zero. 

Hence, the first-order term h*(y) is nonsingular throughout the 
domain of the problem 0 < x < x^. 



APPENDIX D 


THE RADIATION DEPLETED SHOCK LAYER - 
MATHEMATICAL DEVELOPMENT 


The system of equations governing the flow in the stagnation 
region of a radiating shock layer is derived in the text using the 
substitute kernal approximation (see chapter VI). This system is 


f ( T )h * ( T ) 


- | eB”(r) - | k 2 f ( T )h ' ( t) = 0 (D-l) 


2f(Tj)f"(Tl) 


f’(n) 


+ a h( rj) = 0 


(D-2) 


f(0) = 0 


(D-3) 


(v 


= 1 


(D-L) 


f *(^) 


a = 2 f£ 
^2X(l - X) X\e^ 


(D-5) 


K Ta ) = 


This set of equations is subject to the additional condition 

fa P T & "2 k pl t “ T l 

f(r)h'(T) + 6/ 1 kp J B(t)e r| 1 dt 


- \ B(t) + l re 2 ^ 


(D-6) 


(D-7) 
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o 

When the parameter € is very much larger than unity and kp 
the asymptotic solution to the energy equation (D-l) is simply 
B(t) = + CgT. Substitution of this solution into the asymptotic 

form of the integral condition (D-7) gives Cp = Cg = 0. This 
solution obviously does not satisfy the boundary condition hjV^J = 1, 
which indicates that the asymptotic solution is not valid in the 
vicinity of the boundary (shock) at r = r^. This is not surprizing 
when it is recalled that the existance of a thermal boundary layer 
has been established on physical grounds in the text (chapter VI). 

In order to determine the form of the "boundary layer" equation 
valid near the shock, the stretched coordinate 

l = ( T A" T ] en (D-9) 

and the functions 


B b (6)=B(r) (D-10) 

f b U) = f(T) (D-ll) 


are introduced. The subscript b indicates that these functions are 
valid only in the boundary layer. 
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Rewritten in terms of the boundary layer variables, equation (D-l) 
becomes 


-€ 


3n-l 




+ € 




-e 2n |Bj({) = 0 


If n is set equal to 1 and € is allowed to grow without limit, 
the most highly differentiated term will be retained without losing 
the significant term of the unstretched problem. The resulting 
differential equation is 

f b U)h£U)J" + |b£U) = e‘ 2 £k 2 f b (|)h£(|) (D-12) 

The momentum equation, when written in terms of the stretched 
coordinate £ = f - rjj e becomes 


2f b (c)f;’(o 



-e” 2 a 2 hh(0 


(D-13) 


It has been shown that the boundary layer is characterized by 
the parameter € ^ and it would seem proper to expand both the 
boundary layer and asymptotic solutions as power series in 
However, f (tj) (where the subscript a indicates the asymptotic 

ci 

solution valid far from the shock) is not analytic in e" 1 near 
e ■*■ = 0, but is analytic in e”^^ 2 . Consequently, the solutions 
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will be expanded as power series in e r that is 


BJt) 


t a (n) 


L rW 2 V T > 


n=o 


- £ ‘ V2 *. (n) 


n=o 


B b ({) 


Z* Va \ ( » 

Jr-± n 

n=o 


f b (« - £ f b (£) 

l—* Y\ 


n=o 


n 


(D-14) 


(D-15) 


(D-16) 


(D-17) 


In addition, it will be assumed that the enthalpy h is an 
analytic function of B throughout the interval 0 £ B £ 1, and 
from physical considerations it will be assumed that h = 0 when 
B = 0. Then 


h(B) . 




(D-1S) 


where the dot ( * ) indicates differentiation with respect to the 
variable B q . Substitution of expansion (D-l 1 *), (D-15), and (D-13) 
into system (D-l) to (D-7) gives 


ro|o; 
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(2f a (n)f” (n) 

O O 


f a M 
o 


+ a h B 




+ € -1 ^ 2 / 


2f a (n)f" (ti) 

o 1 


2f' (n)f ' (r,) + 2f" (il)f (n) 
a o a l a o ®1 


+ a 2 h(B (t)) \ + ... = 0 

l a oJ a l 


-1/2 . /.x -1 

r ^u> + e - r + e - r ^u; + ... 
a o a l a 2 


= 0 



+ 


1 2 
2 ^ 



= 0 


(D-19) 


(D-20) 


(D-21) 


(D-22) 


26l 



V 2X (i - x ) 


K to 


l 81 



l 







^X 


' KH “ KwK-' lv ' 


A °B a (t)e 2 ^ dt - |B (r) + g 
o o * o 


+ £ ' l/a (l *p X A ° \w ** 2 k?|T " t, “ + ? ^ v%(>, 

-ly^iv'^VW^/N 


a 2 


/ ' p k pl T “ t 

(t)e 2 *1 


+ & ^ViW + 1 ^v%N ■ 1 V t) 


9 2 kp T 

+ V r^e r + 


S r 2 e + V (r)h Pa B X (t) 


07 O 


♦ § V 5 ^ 


T a -T 

<0 



6T, 


[V s) - \(\) 


d| ) + . . = 0 


where 


r = r 
o 


r T & - | let 

- r w *P J °\ < t)e dt 

° O O 


. . . s 0 

(D-23) 


dt 
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r, = 


Vp / A ° \ Me 2 *** dt 


+ **\\[\) 


(D-26) 


IV = r 


pV - I kpt 

'v 1 ^ J B a < t)e dt * r vV* B 

u o 2 


Ap - 2 


'’a') 
, o J 




„ ? T B /- A j. ,. v ~~ 2 Va„ r 


®b ^ 
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o\ o/. 


^ * 
^5 
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The shock layer optical thickness is determined from the condition 




' €T . 





d|\ + ... = o 


(D-2S) 
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The corresponding boundary layer equations are obtained by 
substituting the expansions for and into equations (D-12) 

and (D-13) with the result 



2f b (Of" (E) - 

L o o 


o 


y £ ' 1,2 ( 2fb 


■~.J5 t / f\ii! ! >\ •• / <r \ -r> ! y u 

"b vs/i b, vw T " A b ' i/A b 

o 1 o 1 J 


+ e 


Vs 0 c 


( £) + •••> + a: 0 


f b (0) + s” 1 2 f b (0) + € _J - f b (0) + ... = 0 


-1 


(D-30) 

(D-31) 


2 64 


U) - lU «' 1/2 / 


lim f ( 5 ) -lV+€— f lim f b 


!>"} 


+ e _1 f lim f, (5) + ^ 


L 5 *♦ » 2 


V2X(1 - X) 


* + • • • =0 


B ( 0 ) + e ' 1 / 2 B ( 0 ) + e " 1 R ( 0 ) + ... = 1 

o 1 2 


(D- 32 ) 

(D-33) 


lim B, (ft) -B a ^-. A 


T x A + € -1 / 2 


£ -* « 0 


oV. o, 


f lira B, (5) 

u - « i 


B 


i(\) - \ B % (\)} + £_1 {c V 5) 


B a.f T A^ - t A. B a, f T A - V E a f T A 'l 
o; ± ±\ oj d. o\^ o/ 


!\ ) 

ov o; 


K K (T A )i+--0 

1 o V oJ 


(D-34) 


f lim B^ (|A + e’ 1 / 2 f lim B^ ( |)\ 

1 00 o J 00 1 
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Systems (D-19) to (D-24) and (D-29) to (D-35) lead to a set of 

recursive systems for the solution of B , f , B, , and f. . 

a i a i b i °i 

Zero-order solutions .- The differential system which describes 
the asymptotic solutions to zero-order in the small parameter e“ ^ 
is 


2t. UK (n) 

o o 


b; (t) = o 

o 



f (0) = 0 
a o 





(t) 


(D-36) 

(D-57) 

(D-58) 

(D-59) 

(D-40) 


3 

2 


r e 
o 



= 0 


(D-41) 
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The solution is 


B a (r) . 0 



. a tefr — 1 x) , i ♦ \iix(T 


X) 


« - 1 (o) 


The zero -order boundary layer system is 


, (« + (i) - 4 0) * 4 a) i 

O V 0/ o o 

2f b (S)fJ (5) -[f£ (of « 0 

° ° I ° J 


f b (0) = 1 
o 


lim f (|) = l 
i -» » o 


Bfc (0) - 1 

o 


(D-42) 

(D-43) 

(D-4)|) 

(D-45) 

(D-46) 

(D-47) 

(D-48) 

(D-49) 


(D-50) 
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11m ^ (|) = 0 
| -* eo o 

lim B' (6) = 0 
l -» » o 

The solution to this system is 



o 


(D-51) 

(D-52) 


(D-53) 


f b (0-1 (D-54) 


First-order solutions . - The differential system which describes 

-1 2 

the asymptotic solutions to first order in e is 

B" (t) = 0 (D-55) 

a l 


2T) 2 f" (tj) - ki)f • (n) + bt (n) + 8x(i - x)h(o)B ( n ) = o 
®i ®i ®i a i 

f (0) = 0 
a l 



(D-56) 

(D-57 ) 
(D-58) 


(D-59) 
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B (t)e 
a l 


||t-r|k p 


dt - 2B ( t ) + i 
a l * 


r l 6 


’I V 


+ (D ' 60) 

The system for determining B (t) is identical to the system 

a n 

X 

for obtaining B (r). Therefore, 
a 

o 

B (t) = 0 (D-6l) 

1 

and equation (D-56) becomes 


2rj 2 f " (ti) - 4 i}f ' (n) + 4 f (rj) = 0 
a l a l ®1 

subject to the boundary conditions (D- 57 ) and (D-58) . Inspection of 

the preceeding equation indicates that condition (D- 57 ) is satisfied 

automatically, so that another independent condition or equation must 

be specified in order to obtain a nonarbitrary solution for f (q). 

a l 

This condition can be obtained from the differential system for terms 
of second order. 

The system which determines the first-order boundary layer 
solutions is 


*(\)V 5) + tOV 5) V 4) 


+ lv< 5) = < 4 1) + < 4 2> 6 


(D-62) 


269 


y o - o 

(D- 63 ) 

y o) = 0 

(D-6U) 

lim f b (0 = 0 

— > CO 1 

(d-65) 

B (0) = 0 
1 

(D-66) 

J 3 

H 

\m 

11 

o 

(b-67) 

lim (6) a 0 

(D-68) 


£ -* » 1 


The solution to this trivial system is, of course. 


y « - ° 


(D-69) 


fv (0 = 0 


(D-70) 


Second -order solutions . - The differential system which determines 
the second-order asymptotic solution is 


B" (t) = 0 
a 2 


2r] 2 f" (tj) - 1+T)f ' + bf (r,) = -tj 2 

a 2 a 2 a 2 o 


(B-71) 


2f (ri)f" (n) 
a l a l 


+ ** 


h(0)B (rj) 
a 2 


(D-72) 
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(D-75) 


2 . > 

2 *P J 


r 


, , - |kp|t-T| 

B q (t)e dt 

a 2 


2B a 2 < T > * I r 2 e 


" 2 k P T 



(D-76) 


The solution to this system is 



(t) 


2 p + K 1 - r w)*p\ 


(D-77) 


The second condition for the quantity f , (l> can be obtained 

a l 

by evaluating equation (D-72) at tj = 0 , which gives 


f ’ (0) = a \| h(0)B (0) 

“l 2 


(D-78) 
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With this condition, the solution for f (tj) can be completely 

®1 

specified with the result 


f (n) = a\[h(0)B (0) ( | V (D-79) 

i. 2 


(jl\ 

1 -fJL\ 

Pa 

Pa 

V o / 

V °J 


In addition 


and 


»V V h(0)E (0) 

A o a 2 




*1 ’ ^(O) 


(D-80) 


(D-81) 


As before, it is necessary to look to the next higher order 
system in order to obtain a second independent condition or equation 

for f (tj). This condition is 

a 2 


■ft fn) - n 

Q \ W / “ ~ 

a 2 


/ n _Ao \ 


In order to solve equation (D-72) for f (tj) it is necessary 

a 2 

to express the optical path length t as a function of tj. This 
can be done with aid of the definitions of t and tj. The result is 


t. 0(e 

A I tj 
o\ 'A / 


T = 


(D-83) 
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* 


» 


* 


Now, equation (D-72) can be written in the form 

2 n V ( n ) - V' (n) + (t,) 

& 2 a 2 a 2 



?(* * r ,K T A, 

1.1 + |(l - r v ) Vi - " 

o 


The solution to this equation is 
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(D-85) 


4 


273 


and 


** «^( 0 ) 


[^( 0 )] 


K ( 0 ) 


A _2 a 0 
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1/2 

€t a +6 t. +t. 
n A A A 

— J_ -A, 


1 


K p 1 f®h (1)1 

~ l "O j 

^(o) 
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d! 


(D-86) 


Radiative flux and standoff distance . - In chapter V, it was 
shown that with the substitute kemal approximation the radiant flux 
passing into the cold wall can be expressed in terms of the black 
body emissive power of the gas adjacent to the wall. Substituting 
the solution for B & ( 0 ) into this expression (5.5) gives 

R 



(D-87) 


The ratio of the shock standoff distance to the shock standoff 
distance for radiationless flow is given by the condition 




Expanding this expression in powers of 



gives 


(D-88) 


2Jk 


A a € " 1< \ % 


Ho) + h(0) 


•p (0) 


B a (0) 
a 2 
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X 


1/2 

€T a +€ i/£; T A +T 


*1 


h K u, )^(V s, ) as 
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APPENDIX E 


LIST OP SYMBOLS 


Symbol 

Definition 

Units 

a 

constant defined by equation (2.54) 

none 

OL 

n 

n th-order coefficient in the pertur- 



bation expansion of the constant a 

none 

a* 

constant defined by equation ( 4 . 15 ) 

none 

B 

black-body emissive power, B = — T*; 



chapter II only 

/ 2 

erg /cm -ster-sec 

B 

nondimensional black-body ^missive 



power; except for chapter II 

none 


Planck function, defined by equation 



(2.11); chapter II only 

erg /cm^-ster-sec 

*A 

nondimensional Planck function; except 



•fz-v-r* aVio>s *r*r 

•*> UCi, XX 

none 

5 

nondimensional black-body emissive 



power; chapter H only 

none 


nondimensional Planck function; 



chapter II only 

none 

B o 

B 0 (’i ) = B[b o (n)] 

none 

\ 

VD - 

none 
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B 2 

value of the nondimensional black- 



body emissive power in the interior 



of an optically thick shock layer 

none 

B v 

value of the nondimensional black- 



body emissive power in the gas 
adjacent to the wall in an optic silly 
thick shock layer 

none 

5 

• 

constant defined in chapter V 

none 

5 

constant defined in chapter V 

none 

B 

a 

nondimensional black-body emissive 



power in the interior of a radiation 



depleted shock layer 


B 

a 

n 

n th-order coefficient in the pertur- 


bation expansion of B & 

none 

®b 

nondimensional black-body emissive 



nAt.rav -? v> 4* Vi "V, 2 _ __ 

M.XX VUC 0UUUA. UUUliUtU ^ tJJT 



in a radiation depleted shock layer 

none 

=b 

n 

n th-order coefficient in the pertur- 


bation expansion of 

none 

B* 

constant in formula (6.3k) 

none 

b 

constant defined by equation (5.26) 

none 

C 1 

constant defined by equation ( 4 . 19 ) 

none 

C 1 

constant of integration 

none 

C 2 

constant of integration 

none 
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c 


E 


n 


F 


''1 



f 


f* 

n 


f 

a 



velocity of light 

exponential integral function of order 
n 

nondimensional stream function defined 
by equation (5. 51) 
function of h defined by equation 
(2.5*0 

function of h defined by equation 
( 2 . 54 ) 

function of h defined by equation 
( 2 . 54 ) 

nondimensional stream function defined 
by equation (2. 55) 
n th- order coefficient in the 
perturbation expansion of f 
n th- order coefficient in the P=!=K 
expansion of f 

nondimensional stream function in the 
interior of a radiation depleted 
shock 

n th- order coefficient in the pertur- 
bation expansion of f 

a 

nondimensional stream function in the 
shock boundary layer in a radiation 
depleted shock layer 


cm/sec 


none 


none 


none 


none 


none 


none 


none 


none 


none 


none 


none 
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g 

g 


fit, 


n 


H 


h 


h 


h 


XX 


n 


h* 


n 



n th-order coefficient in the pertur- 
bation expansion of f^ 
stream function} chapter II only 
nondimensional stream function in the 
viscous boundary layer; defined by 
equation (2.85) 


n th-order coefficient in the pertur- 
bation expansion of g 
nondimensional enthalpy defined by 
equation (3*55) 

static specific enthalpy; chapter II 
only 

nondimensional enthalpy; except for 
chapter II 
Planck's constant 

« -4“V> _ nx nrtp -P-P-t /> -t rvv-,4- -f — 4- , t*, r. ■«+■ 

J-J- UU- WbXAXUJLOUU JLX1 UUW OUA “ 

bation expansion of h, the non- 
dimensional enthalpy 
n th-order coefficient in the P-L-K 
expansion of h, nondimensional 
enthalpy 

n th-order coefficient in the expansion 
of h, the static specific enthalpy; 
chapter II only 


none 



none 

none 

none 

erg/g 

none 
erg /sec 

none 

none 

erg/g 
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h 


h 




XI 


n 


h* 



h 


w 


h 


w 


I 


nondimensional enthalpy; chapter 
II only 

average value of h, the non- 
dimensional enthalpy; defined by 
equation (4.5) 

n th-order coefficient in the pertur- 
bation expansion of h 
value of h for which 7 (the 
exponent in the correlation 
formula Kp(h) =■ Ch^) changes 
total enthalpy; chapter II only 
value of the nondimensional enthalpy 
in the interior of an optically thi ck 
shock layer 

static specific enthalpy in the gas at 
wall conditions; chapter lx only 
nondimensional enthalpy in the gas 
adjacent to the wall; except for 
chapter II 

approximation to the ratio of convective 
heat transfer to the stagnation point 
in a radiating gas to that in a non- 
radiating gas 

divergence of the radiant flux vector 
chapter II only 


none 


none 


none 


none 

erg/g 


none 


erg/g 


none 


none 
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I nondimensional divergence of the 

radiant flux vector; except in 

chapter II none 

I nondimensional divergence of the 

radiant flux vector; chapter II 

only none 

I n th-order coefficient in the pertur- 

bation expansion of the nondimensional 
divergence of the radiant flux vector, 

I none 

I nondimensional divergence of the 

radiant flux vector defined by 

equation (3.36) none 

i nondimensional enthalpy in the viscous 

boundary layer; defined by equation 

r-i Ah \ ^ 

\ c- J XJt vu« 

i n th-order coefficient in the pertur- 
bation expansion of i none 

J nondimensional divergence of the 

radiant flux vector in the viscous 

boundary layer; defined by 

equation ( 2 . 87 ) none 

J n n th-order coefficient in the pertur- 
bation expansion of J none 
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h 

k 


k 


eff 


*P 


Pe 

Pr 

P 



c 

<1 


specific intensity of radiation 
mass emission coefficient 
Boltzmann's constant 
effective coefficient of heat 
conduction including energy 
transport by molecular collisions 


and by diffusion of reacting 
species 

Bouguer number* kp « 

s 

direction cosine between the 
direction of a bean of intensity 
and the i th-direction 
Peclet number 
Prandtl number 
pressure 

stanaaxa pressure of air, l.ul^ x 10 


erg /ca?- ster- sec 
erg /g- cm- ster- sec 
erg/°K 


erg/cm-sec-°K 

none 


none 

none 

none 

dyne / cm^ 

/ 2 

dyne / cm 


n th-order coefficient in the expansion 
P 

i th- component of the combined radiant 
and conductive heat fluxes 
n th-order coefficient in the 
expansion of ^ 

component of the conductive heat flux 
vector in the tj- direction 


dyne / cm^ 


erg /cm -sec 


erg /cm -sec 


erg /cm -sec 
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% 



R 


% 



R 




R 



T 


component of the conductive heat flux 
vector in the i th-direction 
component of the radiant heat flux 
vector in the tj- direction 
component of the radiant heat flux 
vector in the i th-direction 
component of the monochromatic 
radiant heat flux vector in the 
i th-direction 

nondimensional rate of radiant heat 
transfer to the vail 
gas constant for air, 2.882 x 10^ 
body nose radius 

shock radius in the vicinity of the 
stagnation streamline 
Reynold's number 
position coordinate 
defined by equation (5*13) 


:/r- 2 


erg /cm -sec 


erg /cm -sec 


/ 2 

erg /cm -sec 


;/=- 2 


erg /cm -sec 


none 

cm 2 /sec 2 -°K 

cm 


cm 

none 

cm 

none 


defined by equation (3.l8) 


none 


reflectivity of the vail 

radiation source function 

position coordinate 

nondimensional variable of integration 

temperature 


erg 


none 

/cm^ -ster- sec 


cm 


none 


K 


standard temperature, 273.16 
temperature of the vail 
nondimensional variable of integration 
component of gas velocity in the 
r-direction 


VS 4* Vs as* Am a AP4 A 4 a A i. 4 u. I.U A 

UMA-UAVACA WCUiViCUO Ali UiC 


expansion of u 

component of gas velocity in the 
i th-direction 
volume 

free- stream velocity 
component of the gas velocity in the 
z- directi on 

n th-order coefficient in the 
expansion of v 


* r. 0 4“ a -5 w +Vi a 4> maa a aJ — ^ , 

Xll WUC u l nun 1 m inrru JJXC Ui C p 


chapter III 

normalized Dorodnitsyn coordinate, 
x = ty • chapter IV 

n th-order coefficient of the P-L-K 
expansion of x> chapter IV 
value of the transformed coordinate 
for which tj = 0; chapter III 


°K 

none 

cm/ sec 

cm/sec 

cm/sec 

._3 

cm 

cm/ sec 
cm/sec 


a/i 


cm/ sec 


none 


none 


none 


none 


coordinate in the transformed 
planej chapter IV 
position coordinate 
local angle of inclination of the 
bow shock from the stream direction 
mass extinction coefficient 
inverse Boltzmann number, 

r » 

Euler's constant, y ■ 0.577216 
exponents in the correlation formula 
for Kp (equation (4 .18) ) 
shock standoff distance 
shock standoff distance for non- 
radiating shock layer 
ratio of shock standoff distance 


for radiating and nonradiating 
shock layer, A = A^A^ 
displacement distance for the viscous 
boundary layer 

exponent in the correlation formula 


B 


5 


n th-order coefficient in the pertur- 


none 

cm 

none 

2 1 

cm*"/ g 

none 

none 

none 

cm 

cm 


none 


cm 


none 


bation expansion of the displacement 
distance, & 


cm 


none 


radiation cooling parameter, 

6 - 4oT s V°» W - 

transformed nondimensional Dorodnitsyn 
coordinate in radiation depleted shock 
layer 

transformed optical path length in 
optically thick shock layer 
Dorodnitsyn coordinate defined by 
equation (2.J4) j chapter II only 
nondimensional Dorodnitsyn 
coordinate 

nondimensional Dorodnitsyn coordinate} 
chapter II only 
variable of integration 
location of the shock in terms of the 
Dorodnitsyn coordinate} chapter II 
only 

nondimensional location of shock in 
terms of the Dorodnitsyn coordinate 
nondimensional location of shock in 
terms of the Dorodnitsyn coordinate} 
chapter II only 

n th- order coefficient in the pertur- 
bation expansion of n. 

A 


none 

none 

g /cm 2 

none 

none 

none 

I 2 

g /cm 
none 

none 

none 


3 * 
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n th-order coefficient in the P-L-K 
expansion of n 
9 n constants defined by equations 

(B.lla) and (B.llb) 
mass absorption coefficient} 
chapter II only 

k^ nondimensional mass absorption 

coefficient 

nondimensional mass absorption 
coefficient} chapter II o nly 
Kp Planck mean mass absorption 

coefficient} chapter II only 
Kp nondimensional Planck mean mass 

absorption coefficient 
Kp nondimensional Planck mew mass 

absorption coefficient} chapter 
II only 

Kp n th-order coefficient in the pertur- 

n 

bat ion expansion of K p 
k r Rosseland mean mass absorption 

coefficient 
X wavelength 

^ boundary layer parameter, A = Pe 

s 

4 coefficient of viscosity 

4 f second coefficient of viscosity 


none 

none 

2 / 

cm / g 

none 

none 

2 / 

cm / g 
none 


none 

none 

2 / 
cm / g 

cm 

none 

dyne- sec /cm 2 
dyne- sec /cm 2 


transformed nondimensional 


Dorodnitsyn coordinate in viscous 


boundary layer none 

thickness of the viscous boundary 


layer in terms of £ 

none 

density 

g j CuP 

standard density, 1.288 x 10"^ 

glen? 

density in the interior of an 


optically thick shock layer 

glen? 

n th-order coefficient in the 


expansion of p 

g/cn? 

Stef an- Boltzmann constant, 5.669 x 10 - ^ 

1 2 

erg /cm -sec- 

area 

2 

cm 

mass scattering coefficient 

2/ 

cm / g 

transformed monochromatic optical 


path length in viscous bo”™ :,D ' , 'y 


layer 

none 

thickness of the viscous boundary 


layer in terms of 

none 

optical path length in a gray gasj 


chapter II only 

none 

normalized optical path length in 


a gray gas 

none 

normalized optical path length in 


a gray gasj chapter II only 

none 


288 



It 


n 


ij 


4> 


0 


n 


X 


t 


n 


monochromatic optical path lengthj 
chapter II only- 

normalized monochromatic optical 
path length 

normalized monochromatic optical 
on - cnapter ijl oiuy 
monochromatic optical path length 
in the s- direction 
thickness of the viscous boundary 
layer in terms of 
shock location in terms of r 
n th-order coefficient in the pertur- 
bation expansion of r A 

A 

shock location in terms of t. 

A 

component of the viscous stress 
tensor 

functions of T defined by 

equations (B.9a) and (B.<7b) 

constants defined by equations 

(b .10a) and (B.lOb) 

density ratio across the normal 

shock, X * p ip 
m/'s 

constants defined by equations 
(B.12a) and (B.12b) 


none 

none 

none 

none 

none 

none 

none 

none 

dyne / cm^ 
none 
none 
none 


none 


none 


solid angle 

^ constant defined by equation (5.22) 

^2 constant defined by equation (5. 30) 

Subscripts 

s indicates value of dimensional quantity 

at normal shock equilibrium conditions 
M indicates value of dimensional quantity 


none 

none 


in the free stream 


